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PREFACE 


ERTS digital data have been compared with, ground-measured bi-directional 
reflectances in 13 separate test node areas. The basic philosophy direction 
of the studies has been 

1. How do ESTS Spectra for any given pixel vary with Season? 

2. How does seasonal vegetation (grass, etc.) respond with time of 
the year, when seen both on ERTS and in comparable ground , 
measurements? 

5 . How does vegetation growing on knovm mineralization, with measured 
anomalous metal content in its leaves, appear in ERTS type spectra 
during its growing cycle? 

1|-. Can ERTS spectra be used to differentiate rock and soil materials 
of economic interest in mineralized areas? 

To do these, experiments we have had to develop software to do the 
following tasks . 

5 Read, unpack and handle ERTS COT tapes on both IBM and PDP-10 
computers . 

6 Develop fully interactive programs for use by geologists (both in 
undergraduate and graduate classes, as well as by professional 
research workers, all of whom have had little or no programming 
experience. Data displays on CRT’s were developed to have a 

map -like format for these purposes. 

7 Develop programs to reduce field radiance measurements to bi- 
directional reflectances, in turn to be correlated with ERTS 


1.1 StJMMARY 

1.1.1 Software (STMSORT) has been developed to read and process ERTS CCT 
tapes to produce the following products, 

1.1. 1.1 Location prints ("shade prints") of approximate scale l:24a000, 
aspect ratio 0,9, on a conventional lineprint. Such prints 
can ..be made from original data, ratios, etc. 

1, 1,1*2 Algorithms written to do unsupervised clustering , edge detection , 
reflectance , and atmospheric corrections . Debanding and 
deconvolution functions were also programmed. 

1.1*1. 3 Images , correctly scaled, shewed and with a 1,0 aspect ratio, 
of either raw data, as Black/T^hite prints, or color composites 
of any three sets forming CIR or other hybrid imagery (eg. Ch 5 
+ R54 + (Ch. 7-Ch.5)). 

1.1.2 Applications to a variety of terrain classification problems were 
completed, with special emphasis directed to the effect of varying 
percentages of vegetated cover on the responses of soils and rock. 

1.1.2. 1 Biomass studies; relating the ground-^measured reflectances to 

2 

the ERTS reflectances, and both to the biomass (dead weight/m 
of vet:etation) - Detailed studies of soil types, vegetation 
species, etc. were made to find correlations. Soil types could 
be clearly differentiated either year in May by the relative 
vigor of grasses in the spring "dieback" period each year. 

Grassland species were too evenly mixed by the feeding cattle to 
give any differentiation of soil type. 

1.1. 2.2 A detailed study of the effect of heavy metal (Mo) poisoning on 
a mixed Pinon Pine arid Juniper "forest" in western Nevada, showed 
a positive correlation of plant growth x^ith Mo content in ashed 
leaves for Juniper* With Pinon Pine, the correlation was negative, 
and the trees showed morphological changes in proportion to Mo content 

Since the conclusion of tills contract effort, considerable work has been performed 
in this and other U,S. mining districts (Yerington, Goldfield, Nevada) and overseas 
(Almeria, Malaga and Rio Tinto, Spain; Karas j ok, Norway etc.) * These works are 
already appearing in publications. 






Hardware was procured ^ Integra tedj and experimental techniques 
perfected to produce meaningful data on reflectance of certain materials 
for comparison with ERTS . 

Significantly in all these projects the efforts have been directed 
to the spot-locality Aground truth" in an effort to derive 1;1 correla- 
tions with ERTS data* 
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A. CHANNEL 4 (-21 DN) AMP 3 
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2. 2.2. z. 19 - 


2 . 2 . 2 . 2.20 - 


C. CHANNEL 7 (-13 DK) AMP 3 

PINE NUT MTNS, NEVADA, 1289-18063, May 8, 1973 72 

A. RATIO 5/4 (-0 DN) AMP I. 

B. RATIO 6/4 (-0 DN) AMP 1 

C. RATIO 7/4 (-0 DN) AMP 1 

PINE NUT MTNS, NEVADA, 1289-18063, May 8, 1973 73 

A. RATIO 5/4 (-41 DN) AMP 1 


B. RATIO 6/4 (-63 DN) AMP 1 
G . RATIO 7/4 (-28 DN) AMP 1 
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C. RATIO 7/4 (-28 DN) AMP 2 
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CH 6 (-42 DN) AMP 3 (Green) 
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B. RATIO 7/4 (^-5 5 DN) AMP 3 (Green) 

RATIO 6/4 (-63 DN) AMP 3 (Red) 

RATIO 5/4 (-66 DN) AMP 3 (Blue) 
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1.0 INTROEUCTIOH . V 

This project was proposed In early 1971 at which time ERTS was only a 
possibility which might produce a few images of reasonable quality for analysis- 
The prime instrument on board clearly was the RB¥ camera, and only a few research 
groups were Interested (more than casually) in the MSS scanner and its concomitant 
digital tape processing. In no way could anyone have foretold the incredible 
success of ERTS-1, still running now in 1976 four years after its launch. 

ERTS^Z (now LMESAT-2) is performing even better than ERTS“1. 

To write an introduction now seems difficult- Almost all of What we 

hoped to achieve has been surpassed manyfold. In the original statement of 

work we hoped to locate and measure seasonal responses from a five square 

kilometer target of grasses and soil on the Stanford Grassland site. Today we 

argue about spectral response from individual 0-4 hectare pixels, and whether we 

have their location fixed to i 50 ra. We have established a 0.95 correlation 

2 

between ground-based ERTS-type raeasurements, and the biomass ratio (gm/m wet 
we ight/dead weight of cover) of the site. We have a 0.70 correlation with 
biomass ratio with several ERTS overpasses from comparable seasons. 

In our other area of primary research, that of exploration for mineral 
deposits, we have progressed from wondering if we could find a "mineral district" 
of several square kilometers, to again debating the spectral reflectance of 
soil and rock areas to the nearest 50 m in position- Now we define the 
mineralogy and chemical composition of the iron-oxide stainings using R54 
ratio — the ratio of the ERTS-derlved bidirectional reflectances of Channels 5 and 4 

Finally, concepts as stretched- and enhanced-images^ sampling, 
shade- and dot-prints (to either 1:62,500, li50,000 or even 1:24,000 scales) 
have become practical realities in everyday use — not only by "high-powered 
computer experts", but at Stanford, by graduate and undergraduate students 
in various academic majors- Therein lies the most significaTlt aspect of this 
study effort, the transfer of technology, from concept, to classroom with 
practical use in field work, performed by growing groups of students. 
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g?.0 BODY OF THE REPORT 

This report covers a funded period fi’om May 1913 through middle 197^5 
with access to LAKDSAT-1 (ERfS) data retroactive to July. 1972. The research 
efforts initiated under tiiis contract (HAS 5-21884) are still continuing 
with various levels of funding from a wide variety of sources. Today the 
project is called STAI^SEARCH — the application of TjAEDSAT data directly 
to mineral exploration, 

Several papers have heen puhlished, others are in press. Four of these 
are included here as Appendices . 

2,0.1 ORGAHIZATIOR 

The report is divided into three main, sections; 

2.1 Direct Image Interpretation 

2.2 Eaduction and Processing of LAMD3AT (CGT) Tapes 
2; . 3 Interpretation of Tapes 

By far the greater effort has heen directed into tape analysis, with excellent 
results well heyopd our original objectives. In fact today (in 1976), it 
is difficult to write about the primitive beginnings in 1972/73. 

Section 2.1 deals mainly with efforts at transmission densitometry . 
using the 4 spectral images, and met with only moderate success. 

Section 2.2 covers the following themes; 

2.2.1 Unpacking CCT tapes and data manipulation 

2.2.2 Spectral Image production - Enhanced B/W and Color 
images . 

Section 2.3 Interpretation of Taped Data 

2.3.1 Atmospheric Effects 

2 . 3.2 Stanford Grassland - Soils 

2 . 3.3 Stanford Grassland - Vegetation 

2 . 3.4 Stanford Grassland - Biomass ; 

2 . 3.5 Californian Grassland -Biomass/Soils 

2 . 3.6 Biogeochemistry ^ Molybdenum poisoning as guide to a 
Mineral Deposit 

2.3.7 Snow mapping 

References, Results and Kew Technology appear at the rear followed by reprints 
of the four published papers . 
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1 DIRECT IMAGE INTERPRETATION 

2.1*0 Feasibility of Using EP^TS MSS Imagery for Spot Transmission Measurements 
2.1*1 Image Densitometry and GCT-Generated Sba.de Prints 

The effort reported here concerns the density/ transmission measurements 
on selected ERTS frames, designed to be related to the MSS tape output (CCT) 
and shade prints* 

The Stanford Grassland test area was examined initially and conventional 
photo interpretation techniques used to obtain relative information over 
the total time period available in our tape collection* 

2.1* 1*1 Image Densitometry-Positive Imagery - 70 ram 

Two approaches have been used- Firstly conventional photo interpretational 
descriptions of features of local geological interest were made, in particular 
noting tone-contrasts and changes* Secondly the 70 mm positives were measured 
(a) using a HcBeth Quantalog TD102 densitometer with a specially prepared, 
smaller (0*7irmi) diameter aperture, and (b) with an ISCO SR spectrbradiometer, 
using a fiber optics probe used to observe the image-plane brightness of the 
transparency when projected with a lantern-slide projector. . These measurements 
can be made up to 6 meters away from the projector, by which time the fiber 
optics subtend only 0.5 mrad, equal to a circle containing above 15 ERTS 
resolution cells (magnification 35. 3x) . See Figure 1- A wavelength of 0*625 
micrameters was used for maximum sensitivity of the ISCO and all 4 ERTS MSS 
images were measured, using the grey scale wedges for .calibration. 

All these sets of data have been tabulated and the spectral radiance (at 

—2 . . —1 ' —1 

the center of each MSS filter) plotted (mwatts cm ster O.lym ) , in Figure 2- 


2.1,2 Discussion 

An experiment was conducted to evaluate the feasibility of using ERTS—1 
MSS imagery for spot transmission measurements on selected sites* Using the 
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15“-unit: linear gray scale on the 70 imn positive imagery as the base of 
comparison^ a correlation was made between the transmittance measured at each 
site on the imagery and the radiance values (coniits) recorded on CCT tape 
by the spacecraft over the same sites. 

The imagery (1075-18170-6) was enlarged (35-3 diameters) by a 750-watt 
projector and focused on a white cardboard screen 6 meters from the lens- 
The end of a 3 mm fiber optic, set flush with and normal to the face of the 
screen, was aimed directly (critical) at the projector lens- An ISCO SR 
spec troradiome ter recorded the light intensities at a wave length of 0*625 
micrometers (selected for maximum scale readings) - See Figure 1- 

Observations were made on 30 sub-sites (Wiltren test site 58, east of 
Travis A. F.B.) which were pre-selected for their differing image "brightness” 
and total areaj all are dried salt lakes* Observations were also made on "average* 
gray backgrounds and on small lakes. The site observations were preceded and 
followed by a series of similar readings on the 15-unit gray scale beside 
the image. Transmission values (ranging from 1 to 15) for each site were 
established from curves (essentially straight lines) of the gray scale values 
(ISCO readings vs 15-unit gray scale)- Unfortunately, during the 3-hour 
experiment, there was a significant decrease in light intensity of the lamp 
between the two series of gray— scale readings, thus producing a time-range 
of transmission values for each site. This range of values was compared xdlth 
the radiance values (counts) obtained from the tapes over the same sites 
(Figure 2) . The values from the tapes had to be averaged because at 6 meters 
the 3irati diameter fiber optic covers an area of about 15 pixels. 

There is very good correlation between the density values obtained from 
projected imagery and radiance values received by the detector on the space- 
craft. A shorter experiment time and closer integration of gray scale and 
site readings would avoid the variation in lamp intensities- See Figure 2. 
Detailed data outputs from the tapes are shown in Figures 4, 5 and 6. 
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2.JL3 CONVENTIQMAL PHOTOITSITERPIIETATIOHS OF EE.TS CHANREL 7 - STANFORD 

The Stanford Grassland test site is about 5000 acres of rolling grassland 
at about 450 feet elevation situated due east of NASA/ARC (Moffett Field) 
on the west shores of San Francisco Bay. The area is now bisected by 
Interstate 280 which iruns approxiinately M-SE along the "grain" of the country. 
Both the grassland and the freeway are clearly visible on channel 5 imagery 
contrasting strongly with the enclosing cities . The site includes 4 
Lakes-Felt, Lagunita, Searsville and Bear. 

(a) ERTS 1075-18173 (October 6, 1972) Sun El 4l°/azljnuth 146° 

Freeway slightly darker than the fields. Felt Lake, black with 
white rim (drying up). Menlo Country Club Golf Course (CCGC) white, 

Stanford University Golf Course (SUGG) white. Santa Clara gravels (Qsc) = 
Monterey shales (Tm) , about equally bright with Page Mill basalt (Tpb.) 

Butano sandstone (Tbu) darker- Stanford Linear Accelerator (SLAC) dark. 
Searsville Lake dark, For tola Road white, Jasper Ridge serpentine (Ksp) = 
sediments enclosing. Lake Lagunita dry. 

Qsc = Tm = Tpb > Tbu 

(b) 1111-18181 (November 11) 60% clouds no use 

(c) 1112-18181 CHovember 29) Sun El 27° /azimuth 155° 

All transparencies badly "cracked", need new prints 

General fields white, no contrast with Menlo CCGS, Stanford GC lighter 
tone. SLAG dark, freeway darker, Searsville Lake dark, Portola Road grey. 
POOR IMAGE. 

(d) 1147-18181 (December 17) 90% clouds, no use 

(e) 1165-18175 (January 4) Sun 24°/ azimuth 151° 

Good contrast. Qsc light-grey peaked on hills /valleys, contrasting with 
Tm, Freeway darker, SLAG darker. West end of SLAC patchy topography. Tra > Tpb. 
Where Tbu without trees approximately equal to Tm. Serpentine on Jasper Ridge 
(Ksp) lighter than sediments, both darker than Qsc, 

Tm > Tpb- : 

Tbu w/o trees = Tm 
.Ksp <■ Tbu, both > Qsc . 
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(f) 1183-18175 (January 22) Sun 26^/azimuth 148° 

Excellanl: contrast . Qsc light speclcled (but different to 1165) 
conti'asting with Tm. Freeway darker, Tni lightest, Tbu little lighter than Tpb, 
above equal to SUGC, SLAG dark, Serprentine (Ksp) lighter than seds, darker 
than Qsc . 


(g) 

1201-18181 

(February 9) 80% clouds, no use 

(h) 

1219-18182 

(February 27) 100% clouds a no use 

(i) 

1237-18183 

(March 17) 40% clouds, ground obscured 

(J) 

1255-18183 

(April 4) Sun 49°/azimuth 134^ 


Good grey scale, overall tone light, Searsville Lake now very dark 
(full of water). Freeway and SLA.C dark-medium grey, trees along creek same 
grey. Topographic relief effect now lost. Field at south end of Felt Lake 
darker tone than that around the lake. Tpb = Tm and only slightly brighter. 
Lake Lagunita now full (black). Qsc even toned (not speckled) , Menlo CCGC 
equals Stanford's UGC in light tone. Alpine Road Quarry Lake and Foothill 
Park Lake clear, black surroundings brighter. Golf courses now more clear 
x?ith fairTrays appearing relative to trees. 

(k) 1273-18183 (April 22) Sun 55° El/azimuth 129° 

Golf courses now the brightest areas with patterns of fairways showing 
clearly, Qsc - Tm “ Tpb.-t Tbu.. Campus trees a little more grey. Topography 
not clearly shown. Jasper Ridge about equal to Ladera. Tbu on ridges a little 
darker than Ksp. 

(l) 1291-18182 (May ID) Sun 60 El/Azimuth 123° 

Image transparency hassy, paper print better^ 

Golf courses clearly brightest areas with Menlo CCGC being brightest, 
?reeway^SLAC darker area, Lagunita Lake medium black. Field to south of Felt 
Lake not so dark in contrast with those around lake. General fields around 
SLAG more bright across creek towards Ladera (Webb Ranch) (Tus) boundary.. 

Natural open grass patches on skyline almost as bright as golf courses. 


2.1.4 MMSURFMEhTS ON SPECIFIC IMAGE SETS (ALL 4 CHANNELS) 

2.i.4.1 Densitometry using McBeth Quantalog (0.7 mm aperture) 
Conventional D log E curves for each channel of Frame 1075-18173 were 
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prepared using special graphs develc^ped by the Principal Investigator T^hile 
in Australia, These show the 15 step, linear-wedge values on a logarithmitic 
abscissa in radiance (xlO watts. Cm -ster see Figure 3. 

Ill the normal manner the density on each step o£ the wedge was measured 
to create the "caaracteristic curve" for that channel. Subsequently then 
the density of selected target areas was read^ plotted, and from the curve 
the radiance read out. 

Because it was not possible to make the aperture disc less than 0.7mm the 
spot size on a 7 Oftiiri formal transparency was very large (2.51an, or 70 pixels) 
and only major terrain elements could be measured. 

. The resultant spectral radiance through MSS bandpasses appear in Table I 

"•2 —1 —3 

(watts. Cm .star , .Olym x 10 ), 

2. 1,4. 2 Transmission Measurements using the Stanford 6-Meter System 

1. Stanford Grasslands Test Site 

Methods used have been explained above. Spectral radiance data 

—2 —3 

are listed in Table II. (Watt. Cm ..ster .Olym.xlO ) for the Stanford 
locality. 

2* . Site 58 a Area East of Travis AFB, California 

An area of small dry salt lakes was located and measured using the 
system,- The data from these measurements have been used to create Figure 2. 

In detail Figures 4, 5 and 6 show sequentially on a USGS Topographic map of 
the area, a shadeprint of the area from CCT tapes and the actual raw numerics 
from a tape-dump (NUMBER); indicating the method of data extraction. 

2.1.5 CONCLUSIONS 

1, A reasonable agreement was found between the two densitometer systems. 

2* The large area subtended even by the 6 m projection system makes the 

system not very useful for our work. In addition diffraction effects from 
using such a small aperture (0.7mm) worry us with the McBeth unit. 

3. We have decided to leave this optical method and concentrate on the taped 
data output, which gives us the ultimate in ERTS resolution (0,4 hectare 
or 1 acre) , 
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TABLE 2. 1.1.1 

SPECTRAL RADIANCE OF ERTS TRAMSPARENCIES USIMG McBETH QUAMTALQG 
(MEASUREMENT AREA = 70 PIXELS) 

DIAGE A IMAGE B 

IMAGE 1075-18173 (W.Cnr2.ster~1.0.01lini.xl0-^ 1183-18175 




OCT 6, 

1972 


JAN 22, 1973 

LOCALITY 

CH.4 

Cll.5 

CH.6 

CH. 7 

CH.7 

JASPER RIDGE (total) 

0.51 

0.34 

0.42 

0.65 

0.44 

BEAR GULCH RESERVOIR 
(W. MENLO PARK) 

0.64 

0.38 

0.50 

0.71 

0.57 

LAKE LAGUNITA 
(DRY) 

— 

0.44 

0.5R 

0.74 

0.57 

SAND HILL ROAD 
(TDS) 

0.62 

0.43 

0.51 

0.76 

0.62 

PELT LAKE AREA 
(QSC) 

0.61 

0.42 

0.44 

0.70 

0.63 

SALT PONDS ON BAY 
MARSH ROAD REDWOOD 

0.79 

0.75 

0.58 

0.51 

— 

CITY 

NOTE: 

AREAS 

COVER 

70 PIXELS 

(2.5 KM WIDTH) 




TABLE 2,1*142 

SPECTRAL RADIANCE OF ERTS TRANSPARENCIES USING STANFORD 6 METER PROJECTION SYSTEM 


(MEASin^MENT AREA = 15 PIXELS. ) 
LOCATION SEE ATTACHED . 
niAGE 1075 - 18173 (W. Cm~2^ster~^* 0 . Qlum xlO 


(OCTOBER 6/ 1972) 


LOCALITY 


San Francisco Bay Salt Ponds 


P ond. A . (Wes tern most) 
Pond B (Striped) 

B 

Pond C 


Stanford Test Area 


Grass cover over: - 
Monterey Shale (Tm) 

Santa Clara gravel (QSc) 

Unnamed sandstone Webb 
Ranch (Tus) / 

Serpentine Jasper Ridge (Total) 
Grass over Serpentine 
Trees over Franciscan 
Felt Lake area (Lake -1 Qsc) 


n.d** 

0.49 

0^42 

0*57 


Golf Courses 


Stanford Driving Range 
Palo Alto Hills Golf Course 
Sharon Heights Golf Course 
Menlo Country Club Golf Course 


not determined j usually not specifically locatable* 































Figure 2. 1.1. 4 Map of Site 58 - East of Travis ATE, Calif prnia 
Scale Approximately 0.5 mile = 1 inch. 
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Figure 2 . 1 . 1. 5 Shade Print of Site 58 Numbered Sub-Sites May Be Found in Numerical Data Sets, Figure 6 

















2,2 REDIiCTION AMD 3PROCESSIKS OF LAtTPSAT TAPES (CCT) 
2.2.1 STAN50RT; ‘Philosophy and Operational Use 



2 . 2 . 1 . 1 Intro duction 


With the increasing application of ERTS multispectral scanner data 
in the field of agriculture, forestry, geology and landuse manageraehtg urban 
studies and hydrology, a need has developed for a system "which "will allow any 
user, from the most experienced to the novice, -with any (or no) level of 
training in computing and programming to use ERTS tapes with minimal effort 
and j,ost. It is also desirable that he have at his disposal all of the 
techniques which have been developed for reduction, interpretation, and evaluation 
of such data, for example, density-eUcing of images, inter-band ratioing, 
clustering and classification techniques, plus the removal of instrumental and 
atmospheric effects. Eor speed, flexibility and ease of operation, an 
interactive program is the most appropriate form. 


For these reasons STMSORT, a fully interactive tape reading program 
was developed. The program has been written so that the operation is self- 
instructive, the user being queried prior to each step and being required 
to input short answers at the terminal, depending on his requirements . The 
user may examine output from earlier stages of the program to assist him. in 
his decision to the current query. Because of the ’tutorial’ nature of the 
program, it has been used also for the instruction of graduate students in 
the use of ERTS data (AES classes 294, 133, etc.). 

2. 2.1. 2 STANSORT Program Philosophy 

The program was developed on a PDP-10 computer at the Institute for 
Mathematical Studies in the Social Sciences at Stanford University. The 
source language is SAIL, an ALGOL- type language taking advantage of many of the 
PDP-10 features. The programs are broken into basic procedures, giving an open 
structure in which any new techniques may be simply inserted with a minimum 
of program change. The IMSSS PDP-10 operates under the TEHEX monitoring 
system, so some small changes would have to he made to the program to 
implement it on any other PDP-10 operating under the DEC operating system. 
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2.2. 1.2.1 Principal Aspect - Detailed Analysis of 20 km (7.7 mi areas) 

- ■ ■ ■ ■ - J i 

It should be stressed that STAHSORT is designed for the detailed analyses 
of a small area on a large scale (approximately 1:20,000) of sections of the 
ERTS images. Additional areas may be added simply to cover more extensive ; . 

scenes. It is assumed that the users will have located (at least to within several ( 

miles) the area in the ERTS image in which they are interested. Such a . , i 

decision may be arrived at by examination of the complete image, or from enlarged | 

L . 

portions (1:210,000 scale) of the image produced from the tape using the Stanford \ . 


image generation program IMAGE (SRSL 74-12, Honey, 1974, 75-4, Levine, 1975). 



2*2*1*3 ERTS System and GOT Tapes 

The ERTS system consists of satellites in a suii'^syiiclironous, almost 
polar orbit 3 at an altitude of 920 Icilometers. Either satellite crpssses 
the equator 3 traveling in a southerly dlxection at 0942 local time. Every 
18 days the satellite covers the same ground track, so that sequential 
coverage of any ground scene is available, with an 18-day period. The 
instrument package on the satellite consists of three ReturnBeam Vidicon 
(!RBV) cameras, a Multispectral Scanner (MSS), a. Data Collection. System 
(DCS), and two wide band Video Tape recorders • 

The RBV cameras, sensitive in the ranges 475-"575 nm (blue-green, 
channel 1) 580-680 nm (red, channel 2) and 690-830 nm (near infrared, 
channel 3) have not been utilized due to an early electrical malfunction* 

These instruments will not be described further* Such 3-band data would also 
require modification to the S TANS ORT (4-band) programs. 

; The MSS is a four-band scanner, with bandpasses of 500-600 nm (green, 
channel 4), 600-700 nm (red, channel 5), 700-800 nm (infrared, channel 6), 
and 800-1100 nm (infrared, channel 7) * These wavelengths are in the reflected 
region of the spectrum - no thermal radiation is detected, although ERTS-C 
(1977 launch) will carry this as a fifth channel* Each of the bands has a 
bank of six detectors to assist the scanning speeds. Across-track (E-W) 
scanning of 185 km is achieved by means of an oscillating mirror in the 
optical path ^ Six scan lines are imaged at once in each of the four bands. 

The. video outputs are multiplexed and either stored on the video tape recorders, 
or, if in range of a ground receiving station, encoded and transmitted. 

Recorded data is encoded and transmitted later. Both of the tape recorders 
in ERTS-1 have failed, hut ERTS-2 is making recording to be dumped* 
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Each image is made from 2340 scan lines each of 185 km length. A 
scan line consists, after processing, of 3240 digital units (picture 
elements, pixels) which are packed as S-kit "bytes on tape. In addition 
56 bytes of calibration information are included at the end of each record 
on btilk tape. Thus for each image there are about 3 x 10^ bytes of 
information. A set of computer compatible tapes for a scene consists of 
four tapes, each tape representing one quarter of an image, a strip 46*25 

■ ; 5 • 

km E-T-? and 185 km Each tape contains 7 ’71 x 10 bytes to date, with 

a 40 -byte scene identification and a 624-‘byte annotation record appearing 
as the first two records. 

Each data record on the tape consists of 2340 bytes of data, representing 
4 channels for one quarter of the image. The data is in an interleaved 
format, as indicated below; 


BYTE HO- 12 

3 4 

5 6 

7 . 8 . 


CHMHEL NO. 4 4 

PIXEL NO. 1 2 

5 5 

1 2 

6 6 
1 2 

7 7 

1 2 








2 

Each pixel in a scene represents an area of 57 x 79 m (approximately 0.4 
hectare, or one acre). 


2* 2,1.4 Operational Use of the STANS0E.T Program 
2.2-1.4.1 Initial Steps 

The user is required to have logged in, mounted his tape and carried 
out any preliminary measurements to determine the location of his area 
(Appendix 1), He is then queried by the mala program for basic information 
regarding the tape. The program reads the Header and Annotation records and 
displays identifying information on the screen -- Scene ID, date of acquisition 
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of data and which tape (of the set of four) is being read* This necessary 
Stage prevents mounting the wrong tape and subsequent (costly) confusion! 

If the tape is correct , then the ”Yes” answer causes relevant information 
from the header and annotation records to be stored in the housekeeping 
file. This sequence of actions is typical-data:, shown from the tape, a 
question, which requires decisive action by the user* 

The user than enters the co-ordinate of the position on the tape 
(image) and the size of the area^ The program checks the coordinates 
to he sure they exist on that tape^ then the tape advances to the area 
which is read in and repacked. The program then enters a series of 
procedures, querying the user and proceeding according to the response. 

The Flowsheet outlined in Figure 1 illustrates the cperations . and their 
normal sequence of execution* This sequence was derived during development 
of the program as the most common one a new user would follow. For the 
more ^perienced user, or for a re-examination of a scene as wxXh. another 
algorithm the '^JUMP! TO*^ operation may be utilized, With this feature 
the operator may skip forward or back to any procedure and break the 
normal sequence. 

2* 2*1, 4.2 Header 

The Housekeeping File stores the user-designated title for the area, 
the ERTS tape identifier and dates, the relevant information from the 
header and annotation records, such as frarae-center coordinates, sun 
elevation and azimuth and satellite heading- This information appears 
at the top of every output as. indicated in Figure 2, . 

The majority of the queries and the effect of the replies are 
evident from the flowsheet. We will proceed with the sequence and a brief 
description of each procedure will be given. More detailed description 
■are presented in Section V. 
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STANFORD REMOTE SE^'81NG LABORATORIES 

STANFORD UNIVERSITY 
CALIFORNIA 
U.S,A, 

TeU C^i5) i|97'"a7'^7 F'.R^HQNEV 

+++t ++ +++ +++++ + + + +++ + + + + + + + + + +++ + + ++++ + ++ + + Tt'{'t ++ + '/ ++■* + + + t+++ + + + + tt + 


ERTS tap® 1273*l8i8300 Tap® nginbep 2 of H 

Expbsyre date 22 APR 1973 

Coopd^natea of frame penter N37**37/W122«0H 

Sgn elevation 55 

Sun aail myth 129 

Satel H te heading 190 

Area begins at pow 1299,p^x«l 1293 of frame 
scale 1 HQRIZ (E-W) 1 r 22^40 f VERT CK-^SJ 1 *16670 


REDWOOD CITY TEST FOR AES 133 CLASS 

n «a wtvaentatMMoifHMiaBiMtstsM 


Shadepr i nt of 5 

ERTS raw data(t debandedf deconvoluted to 'sharpen image* 
Levels and thefr symbols* 


■ a? - 8 
l»9 - 10 
Ms'll - 12 
«=13 w U 
®=15 - 16 
Sal7 f 18 
7.=19 -.20 

^^a=21 22 ' . . 

+=23 9 ZH 
?=25 - 26 
«=27 w 26 
- \=29 30 

,=31 - 32 
=33 * 3^ 

NOTES Minimum value of area i* 7^ D.C, level bf 7 taken from y«lue»*. 


Figure 2.2.1. 2 Header showing Housekeeping Information 
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Some of the procedures interact with each other. For example^ 
histograms of ratios may be obtained by a ”juinp^^ thereby presenting 
the histogram procedure data in ratioed form. Such interactions are 
presented by dotted lines on the flowsheet*. 

For speed and efficiency all arithmetic operations are performed in 
integer arithmetic,, except where the roundoff would seriously affect the 
result of the procedure. A great deal of careful consideration was taken 
to minimisje any unnecessary time consuming operations within the procedures, 
for therein lie ^ the principal cost generators. 

' . After reading in the tape data representing his area, the user may deband 
the data. This function attempts to remove much of the ^striping” of the 
scanner data* 

2* 2. 1*4. 3 Deconvolute 

The data may then be **deconvoluted", a procedure included in an attempt 
to remove the effect of the overlapping fields of view of the scanner across 
a scan line. This appears to provide significant image improvement, 
particularly when searching far fine: details such as roads, rivers, airports, 
etc. 

2. 2*1. 4. 4 Histogram 

Histograms of the data may be obtained at this point* The histograms 
are. almost indispensible for a clear interpretation of the. data, and as an 
indicator for some of the subsequent procedures, such as shadeprinting, 
ratioing and clustering. We feel however that automatic use of histograms 
to control the grey scale of a shadeprint is not advisable , and prefer to 
leave this option available to the user* This is a particularly attractive 
aspect of an interactive program. 
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Shade Print 



2, 2. 1-4.5 

The next procedure* SHADE J PRINT^ preseifts the user with his first look 
at the data in a ^map^ form on the display screen* Typically the experienced 
user jumps" to SHADE initially to make sure of his location and to . quickly 
refine the co-ordinates if necessary. He is queried as to the correctness 
of the current location: if not correct j the program skips to the step of 

asking for a new area; if correct the user may print the shadepriiit with 
various levels of slicing (stretching or enhancement) for any of the channels. 

2. 2, 1.4*6 Edge 

Edge detection is included to enable a search for high-contrast changes 
as may occur at boundaries such as land-water interfaces or, possiblya of 
curvilinear geologic features such as faults or intrusives, or a vegetation 
changes. 


2.2. 1^4. 7 Ratio 

Interband-ratioing allows the user to select the numerator band, de- 
nominator band and the levels a,t which slicing is to be performed. The 
result also may be displayed in ^raap^ form with a shadeprint representing 
the ratio values printed. 

2. 2.1. 4.8 Numerics 

Numeric data from either raw data (or data converted to reflectance) 
may be printed for more detailed examination, by a call to NUMBER. 

2. 2. 1.4. 9 Clustering 

A cluster analysis may then be performed on the data in an attempt to 
separate population classes, or if some information about expected classes 
is available a supervised classification may be performed, . 

2.2.1.4.10 Special Function 

The special functions procedure (not implemented to date) will act as 
a run-time compiler, allowing the user to insert any new (special) function 
by which the data will be evaluated. This will provide an extremely 
versatile and powerful otpion for testing new ideas on manipulation of the data 
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2.2.1.4.11 Reflectance 

The data flag is then checked to see if the data is already in the 
form of reflectance. If not, the user may convert data to reflectance 
provided he has certain standard measurements to insert. After conversion of 
of the data to reflectance the program recycles through the entire sequence 
(D-I)v 

2.2.1.4.12 New Area, New Tape 

Having satisfied himself that no more information can be obtained for the 
current area, the user may then move to another area of the tap e* or to 
another tape and repeat the sequence of operations. 

2.2.1.4.13 Scale of Lineprinteir Output 

A note on the printout format from SHADE! PRINT, RATIO and CLUSTER. The 
N-S scale is approximately 1:19,000, the E-W scale approximately 1:22,000. 

The skewing (3.5 degrees) due to rotation of the earth has not been taken 
into account, as the resulting ’’map^ is slightly distorted, both in the 
magnitude and direction of its axes. This in no way interferes with the 
interpretation steps as all the output is similarly distorted in precisely 
the same way. 

2. 2. 1.5 Detailed Description of Each Procedure 

2. 2,1.5. 1 Debanding (Removal of Banding or Striping) 

Noise on the satellite data appears to take two forms. The first, of 
a fairly randopa nature, arises apparently from digitization on the satellite, 
followed by calibration (which necessitates f loatine point numbers) . For 
Channels 4, 5 and 6, noise is introduced during decompression of the data, by 
usitig the-look up table (Thomas 1973) which has missing values. An 
exaidination of the histograms for the three decompressed channels shows much 
more apparent noise than for Channel 7, thus it would appear ’that the 
majority of the random noise arises during the decompression process. The 
second form of noise, a less random, banding with a six-scan line period, 
probably arises from etror used in the calihration expression for each 
of the tour detectors in the bank of six per spectral band. 
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The first form, by far the least sign Lf leant of the two^ is the most 
difficult to remove. One crude approach is to convolute the data with a 
smoothing function, thereby reducing any noise spikes, but in the process 
effectively ’defbeusing^ the image by lowering the spatial rsaolution. This 
technique obviously also removes most of the banding. Because of its effect 
on the resolution, however, it is an undesirable approach when looking for 
fine details but has attraction in some regional geochemical studies 
(Lisaur, 1975) . For examination of targets with no high-frequency Information 
expected, such as areas of shallow (< 20 m) water to estimate depths^ this 
approach may be employed. This technique was used in our initial studies 
but was abondoned in favor of a "debanding** algorithm* True smoothing is 
being re--introduced^ after debanding^iu the new "STANSOkT 3", as a branching 
option relative to deconvolution. One thus has a three-way branch - raw 
data, smoothed data or deconvoluted (sharpened) data* 

The problem of * striping'* (banding) of the imagery has been considered 
by Algazi (1973), who presents an algorithm for removal of this noise. 

In this technique, the mean and variance for each channel across a whole 
tape record (810 pixels), and for a large number of scan lines of the 
image are calculated. The means and variances for each line are calculated 
as the data is being read in, and comparison with the global means and 
variances, provides an offset and gain factor to he applied to each scan 
line for each channel* This appears to effectively remove the banding, but 
does not affect the random noise discussed above* For intermediate and high 
radiance targets i.e* where this random noise is low relative to the signal, 
the image quality is improved substantially. 

Debanding is generally a necessary step if ratioing is to be performed 
on the data as the noise effects are increased noticably by the ratios. 
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It may be more appropriate to automatically smooth the data using 
AlgaziVs algorithm as the data is read in; at present this is a ‘later 
option^ the debanding factors being calculated, not for every line, but 
as a set using the first twelve lines of the area of the image under 
study j the factors being stored, and printed out for reference* 

If the data is debanded a flag is set, and the housekeeping file 
is modified to indicate debanded data on output* 

2, Z* 1*5- 2 Sharpen (Deconvolution) 

The area sampled by tne scanner across the ground track is nominally 

• . 2 ■ ■ ■■ ■■ ■ ■ 

79.3 X 79.3 meters , but has an overlap because of oversampling. This 

overlap is 21*8 meters, if one assumes constant mirror velocity^^ There 
is, ideally, no overlap between adjacent scan lines , although this may 
be disputed at the extremities of each line due to slight broadening of 
the ground footprint of the field of view* 

This overlap of pixels results ‘.ight degradation of the image* 

This may be reduced dramatically by xvolution^ of the data to remove 

sdme of the effects of the overlap* 

The un-normalized digital function performiTig this operation may be 
represented diagramatically by 



^"The mirror velocity is not constant but is a fixed^ continously varying 
function directly related to the E-W position of the pixel* We make a 
simplifying assumption that it averages 13*75%. 
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For any pixel (across the scan line) , the enb inced value Is given by 

' X,'= N(0.725 X. - 0.1375 (X._, + X )), i= 2...n~l, ( 1 ) 

where N is a normalization gain factor given by 



The normalizatloti factor may be evaluated by the program- For average 
radiance targets it has been found to have a value of 2-0- 

This enhancement procedure has been found to give greatly improved 
quality in recreated imagery^ though for use rath line-printer output 
the enhancement is not imniediately obvious* Image contrast is heightened 
as well as apparent spatial resolution- 

2- 2. 1-5- 3 Histogram 

This procedure provides output on the screen, (and if desired, on the 
lineprinter) the present type of data being processed — either raw or 
smoothed data, enhanced data ox ratioed* The histograms are scaled for the 
width of our line printer to have a maximum ordinate value of 68- The 
frequencies of counts for each level are printed down the right hand margin 
on the line printer. 

The procedure itself is straightforward, but provides a powerful tool 
in the interpretation and planning of processing of the data at later stages 
of the program. As mentioned above we differ in program concept from many 
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other groups and do not use the histogram areas to automatically control 
the next SHADEPRINT step. We reserve this for an operator-decision. 


2. 2. 1.5.4 shade: print (Map-like Output) 

The shadeprinting procedure enables the user to display on the 
screen a symbolic 'map' of his area, with approximate gray scale lines 
being indiated by appropriate alphanumeric characters (no graphic terminal 
is available at present , although we are reprogramming an IMLAC unit for 
this purpose). The shadeprint may be printed in either a replica of the 
screen image, or overprinted twice per raster line to give better gray 
scale representation.wlth the appearance of a over-enlarged newspaper 
illustration. (See below) . FIGURE 2.2. 1.5.1 
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DN"Offset , Two functions nicay be performed on the data at this 
stage to improve the presentation. The first j removal of a DN (digital 
number or counts) offset allows the user to have a lower cutoff level for 
the data. In the case of raw data, for the visible chaniiGls. (4 and5) 
there is a significant contribution due to radiation backscattered by 
the atmosphere. For Channel 4, this may reach 18 digital counts. This 
may be removed to allow better slicing intervals, and to look at another 
useful range of the data. Another user may be interested in only looking 
at very high radiance targets, for ejcample snow and clouds (Itten, 1975). 

All data below 120 may be disregarded, and with intervals of 1, the remaining 
screen counts will allow partial differentiation of wet and dry snow. 

Stretch . The second function involves variations in different 
levels of slicing, or changes in the grey scale increments. Although each 
of the levels are the same (i.e. "linear stretch") this allows the user to 
produce the optimum ^ map ^ of his scene on the screen and in the subsequent 
printout. Varying the DN offset and the slicing interval, allows the user 
to look at specific types of scene- water, forest, rangeland, or snow. 

No attempt has been made to implement the more automatic ^histogram 
corrected’ form of shadeprinting, wherein equal areas of the frequency 
curve are allotted^ to. equal grey scale intervals.. We . felt that some 
small but (possibly important) details may be lost using this technique, and 
each user must make the decision which best fits his data set. 

2- 1.5. 5 Linear 

In most applications of photogeology, linear and curvi-litiear features 
are selected based upon the concept that such boundaries are structurally- 
controlled features^ A search may be made for those boundaries, which show 
as contrast- edges between materials of different reflectivities, using the 
edge-detection algoritlim. In this algorithm an attempt is made to eliminate 
single- pixel noise, by requiring that any contrast change be present simultaneously ‘ 
in at least two of the channels. If a pixel-pair possesses this contrast 
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their position is marked "by an asterisk symbol^ chosen so as to appear 
non-directional to the eye of the observer, and hence will not influence his 
decision of linear joins* 

The user may decide the contrast change (or tolerance level) required 
to 'define a boundary for each of the channels. A* map* then appears on 
the screen, and the levels may be changed to produce greater or lesser 
detail. The product may then be printed. 

We have not used this spectral-lineation detector as much as we had 
hoped. The technique works satisfactorily for obvious delineations. It 
has the advantage of removing unwanted (presimably) information. One 
possible application being Investigated is the automatic "^contouring” of 
shallow (< 20 m) water data to prepare bathymetric maps, by using the 
known spectral aspects of water to define depths in clear sediment-free hays. 

2.2. 1,5.6 Ratio 


The application of interband^ratioing to enhance and accentuate features 

is a well established technique in several remote sensing groups. For 

enhar.cement of alteration zones in mineralized areas ratios of Channel 4 (green) 

over Channel 5 (red) yields some excellent results (Rowan, et. al, 1974) 

CH.7 CH7 

For vegetation studies ratios of or 7 ;;;^ give a good indication of vigor, 

CH5 

although some subtle changes may be brought out by (Lyon, et.al,, 1975). 


The user selects his numerator and denominator channels. Maximum and 
minimhni values for the ratios are displayed, and the user may request a 
histogram of the ratioed data. After examination of the histogram (and, if 
required, printing), a DM level and a slicing interval are selected. A 
shade-printed map is presented of the results which may be modified before: 
and after printing* Any number of ratios may be examined in this 'way, and 
hard-copy produced as desired. 



2, 2. 1-5*7 Ratio of Ratios (RATSiRATS) 

Ratios (or products) of ratios effectively provides a combination of 
the form 

(channel u) X(channel m) 

(channel ±) XCchaunel j) 

This technique was implemented on the request of Lizaurj in an attempt 
to increase discrimination of the data, in particular in alteration zones. 

The results of a study using this method are presented in Lizaur (1975) * 

The procedure is very similar to RATIO - the required channels are typed 
XTif histograms may be examined and a shadeprint prepared, examined and printed 
At present this proceedure is being replaced by a more general (-f, X, or 
algorithm to manipulate matrix-^airs* 

2, 2. 1,5. 8 Number 

This procedure outputs the data in its current form as it occurs in 
memory, i.e. raw data, debanded or deconvoluted, or reflectance data. The 
quantity of output is large, so that the use of NUMBER is only recommended 
when some definite study requiring the numeric information is to be attempted. 

2,2,1*5.9 Cluster 

Clustering of the image Inforniation into unique classes while using a 
fully^lnteractive mode is one of the most important aspects of the program. 

The algorithm used is extremely straight forward and fast, although it lacks 
the statistical rigor of the more usual clustering techniques . Some of these 
techniques will be discussed before describing our clustering procedure. 

Classification of data may be carried out in two modes. One is the 
so-called ‘’supervised^' classification, whereby previously defined populations 
with their representative spectra input to the program, are assigned to these 
known types according to some nearest-neighbor or least-distance criterion. 



This can either involve some preprocessing of the data to eliminate unncessary 
data (Andrews, 1972, Sebestyan 1962) followed by testing with the distance 
criterion, or the use of ’’table look-up” techniques in which a range of data 
is input for each pre-determined class and the data sorted by comparison 
(Eppler, 1974). The second of the two techniques is extremely rapid and is 
somewhat similar to our clustering algorithm described herein. 

The second mode of classification, ”non-supervised" classification 
(or clustering) uses natural differences in the data to classify them into 
arbitrary groups (Andenburg 1971). Several types of rigorous, but time - 
consuming statistical proceudres, are available for the clustering of data 
(frelat, 1974). 

The technique used in STANSOB.T is to divide four-dimensional space 
into hyper-rectangular ’’pigeon-holes” , whose locations in space are 
determined from the data. The significant difference with our approach 
is that the size of the rectangular "pigeon-holes” are determined interactively 
by the user. The data is then sorted into its appropriate position. 

Inspection of the shade-printed ’map’ allows the user to preserve the patterns 
or to redo them with variable pigeon-hole sizes (tolerances). 

In examining an initial area, the user selects the tolerance (or gate-width) 
he will allow around each band of the spectrum of any class. The program 
then begins to cycle through the data, the first pixel being arbitrarily 
assigned as the first class (”A") . The remainder of the data is then compared 
with the spectrum of the first class. If the spectra fit within the tolerance, 
they are assigned as the same class. After examining all member of the array, 
the program recycles, taking the first pixel, unassigned from a previous cycle 
and classifies it with the next symbol, comparing all remaining unassigned 
pixel spectra to the current ’standard’. The program continues to cycle in 
this manner until all data has been classified, or until no classes remain 
(maximum 26 with all others assigned to a blank symbol). A classification ’map’ 
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appears on the screen, and the user decides if the result is satisfactory: 
if yes, then he may print; if no, the tolerance may he varied- The final 
set of ^standard* patterns i-e- the mean of each class, are printed out, 
and stored for classification of subsequent areas. The printed output may 
be. used to evaluate each class- . Sometimes these spectra are found to be 
closely similar, although just outside the tolerance range. The technique 
appears to over’-classify some data T^hich could possibly be merged into one 
large class, perhaps during a subsequent hand-analysis stage (coloring, etc). 

The classification is sensitive to the starting spectrum. If it is 
started at a different position, the ^standards ^ and classes will be different. 
However, the gross clustering remains very similar - only the assigned symbol 
and partially, the 'standard^ spectra for an area will vary. 

It should be emphasized that this clustering procedure is used mainly for 
smaller areas, for example summing up to 4000 to 5000 pixels. It can be 
modified for examination of. complete tapes or sets of tapes, and would remain 
a very fast technique, but at present we are more interested in detailed 
analyses of small areas (^200 square miles) with mining interest. 

The clustering procedure may be used in four modes, un~normalized or 
normalized on raw data, and similarly on reflectance data. In addition, for 
reflectance data, the procedure may be used in a classification form to search 
for materials with known reflectance spectra which are typed into the program. 
Such spectra would, perhaps, result from field measurements in the area. 

It is found that for surfaces with undulating topography, that the spectra 
of similar materials are virtually identical, but the relative magnitudes of the 
spectra vary according to whether the pixel is on a sun-facing slope or on a 
partially-shaded slope. For moderate variations in slopes it . is found that 
normalizing (or ratioing) to one of the channels removes most of the effect 
due to topography. Without normalization, two similar spectra may be classified 



differently. For extensive areas of uniform slopes, or horizontal areas 
(coastal plains, lakes or ocean scenes), the necessity for normalization 
is removed, so computation times may be reduced by avoiding normalization. 

Both are present in STANSOBT as readily-available options. In addition, 

in the un-normalized mode, the user may choose to discard a channel, particularly 

if that channel is excessively noisy. 

In conclusion, it should be emphasized again that this algorithm was 
chosen primarily for its rapidity and comparatively low cost. Testing 
of the procedure against other techniques has been limited, but one case 
(Itten 1974) produced virtually identical clustering results to a procedure 
in LAESYSsat a fraction of the cost, in a much shorter time. 

2.2.1.5.10 Special Function 

Using this procedure the user could type in any form of manipulatory 
algorithm he wished. This procedure has not been implemented as yet, but 
its basic aim is to provide an extremely versatile facility by which to 
evaluate functional approaches to enhance the data, and for feature extraction, 
without the necessity for xnriting a separate program for each function. 

2.2.1.5.11 Reflectance 

The digital data present on the ERTS tapes represent the radiance signal 
received at the sensors. This signal is composed of atmospheric backscattered 
radiation, and radiation reflected from the ground target. The magnitudes 
of both of these signals are a function of sun angle, and their relative 
values may vary significantly. For inter-season comparisons therefore, it 
is preferrable to have the data in a more standard form, such as relative 
reflectance. To achieve this, several standard targets must be available 
within each scene. These targets should be extensive (over 9 pixels in size) 
with near-constant reflectivity throughout the year, and easily locatable. 

Once the standards are chosen their bidirectional reflectivities should 
be measured relative to some laboratory standards (barium sulfate or smoked 


magnesium oxide) • The user may then search for these targets on the 
tape, extract the digital readings for the areas and derive regression 
coefficients for the conversion of the satellite data to reflectances,. 

It is desirable to have several such areas of known reflectance within a 
scene so that the converted data may be checked, as atmospheric effects 
over a scene may vary (iDuggin, 1974). 

With the data now converted to reflectance the user may then perform 
all of the previously described functions on the data. A reflectance 
^flag* is set by the program for each area to avoid attempts at conversion 
of areas already in a reflectance form. 

2.2.1.5.12 Completion of a Study 

Having completed all of the required steps above, the user is then 
asked if he wishes to examine another area on the current tape. If so, 
then he inserts the relative position and sis e of the new area. The system 
checks the coordinates to be sure they axe valid for that tape and then 
the tape is repositioned. Data for the new area are read in. Should another 
area not be required on the current tape, the user may then look at another 
area on another tape, which must be mounted. For both of these steps, control 
is stored (such as clustering standards and ranges) along with the relevant 
header information. 

If the user does not wish to continue, he may exit from the program* 

Any printout the user may have requested exists as separate files on 
disc, and must be printed with a separate program. 
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2.2.2 m INTERACTIVE PROGBAM FOR PRQUUGING CQMFUTER-EMANCED ERTS IMGES 
(^^IMAGE*^) PART II - METHODOLOGY DEVELOPMENT 

2 . 2 . 2 ^ 1 Introduction 


This paper covers the operational or applications phase of 
the imagine software program during which the methodology was finalized. 

A series of black and white and color-composite images of several 
areas in California and Kevadaj generated during this phase^ are 
presented and discussed. 

It should be emphasized, that it is the operational’ 
simplicity and the interactive nature of the program, utilizing the 
PDP-10 computer and a keyboard controlled CRT display that represents 
the core and the strength of the program. The best capabilities of 
both man and machine are utilized* With IMAGE it is possible for the operator/ 
investigator to Initiate functions and react quickly to each step 
in the program. Parameter changes may be made rapidly to optimize 
the final enhancment, a tape generated and the required images made 
by use of the DICOMED Image Recorder*. Used in conjunction with the 
S.RSL STMSORT Program (Honey, Lyon, 1974) for preliminary investigative 
purposes, computer- enhanced ERTS images may be produced economically : 
and quickly* 

A general discussion of the initial phase in which the need, criteria 
and evolution of the Stanford Remote Sensing Laboratory interactive program 
CU-IAGE) for producing computer-enhanced ERTS images are covered by 
F. R. Honey in SRSL Report 74-12, 

*Made available to us by the courtesy of NASA Ames Research Center, 

Moffett Field, California. 
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.2. 2,2,2 Image Generating. Program 5?unctxo.ns 

The image generating program was designed to reformat the interleaved 
image data on tUe ERTS GGTs in such a fashion that it eould he manipulated, 
as desired, to enhance areas o£ interest and then generate digital tapes 
that are compatible \rLth a CRT recording instrument (siich as the DICOMED 
D47 Image Recorder). These tapes then contain the necessary intensity 
information to produce enhanced black and white, or, by using 2 or 3 such 
black and white images, to make color composite images of the areas of 
interest. The images obtained are recorded on standard 4X5 inch 
Polaroid or Sktacolor film, covering an area of 300 pixels 
horizontally and 250 vertically at a scale of 1:210,000 approximately. 

As indicated in SRSL Report 74^12 the image program, which is an 
interactive one, designed for use with a PDP-10 computer, locates the 
area of interest and then makes a series of mathematical functions and operations 
available which the user may apply to study the area before generating the 
final image* These functions will be discussed below in a sequential 
fashion as they ■'po* ar in the program^ 

2- 2. 2. 2,1 The am initially makes geometric corrections in the image 

format due to the lateral skewing of the progressive lines of the scanning 
device, caused by the orbital inclination, rotation of the earth , and the 
higher sampling rate in the scan direction vs, the flight direction, 

2, 2, 2* 2, 2 After locating the area of interest and specifying its extent, 
a debanding function is made available* This option selects 6 lines of 
data, a full 810 pixels wide and determines their 6 "average means”. 

The ” deb ending factors” then are the quantities necessary to make all 
the average means up to equal average intensity. The striped pattern 

40 






is due to the imperfection of the on-board sensor radiometric 
calibration which results in a banding pattern repeat at six line 
intervals in the recorded images* 

2- 2. 2* 2*3 A deconvolution function is also supplied which tends to eliminate 
the scan line overlap of pixels due to the over-high sampling rate, thus 
s trengthening their contras t * 

2. 2* 2, 2, 4 Since the scanning system is designed to cover a large dynamic 

range (due to wide variations of scene albedo and sun angle) the brightness 
range of any particular image may only cover a part of the dynamic range* 
Contrast enhancement of the image is made possible by stretching this 
brightness range. This program presently provides a linear stretch 
(uniform contrast increase) over the entire range of the image. 

2. 2. 2* 2. 5 To facilitate determination of the desired limits, a histogram 
of the scene, showing the frequency distribution of the image digital 
numbers (scene brightness) now stretched between 0 and 255 M may be 
obtained. Study of the histogram then enables setting the most desirable 
limits for optimum contrast or enhancement of a particular scene 
element or area. To assist In this determination maximum and minimum 
values, means, and standard deviations are also provided in a print out. 

A DN value subtraction may be made x-rhich effectively zeros the left 
hand extreme of the histogram. 

2. 2.2* 2* 6 Amplification of the remaining values is made linearly to 

coincide with the 255 DN value limit (or beyond), at the operators discretion 
should more white or more black seem desireable for a particular problem. 
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2.2. 2, 2. 7 Itidiyidual cbannels may be enhanced and imaged or ratioed 
pixel“by^pixel to show the variations in the slopes of the spectral 
reflectivity of the two bands. Ratio stretching tends to enhance 
the spectral reflectivity differences and also minimize radiance differences 
due to albedo and topography or slope. The histogram capability^ 
the DN subtraction and linear amplification are also provided to 
facilitate enhaneement in the ratio mode- 
2- 2- 2. 2. 8 After generation of the image tapes they are then utilized to 
control the output recording of the DICQMED P47 Image Recorder to 
produce enhanced black and white images. This unit may also be utilized 
to produce enhanced color composite images of two or more channels or 
ratios of channels in which the color variations represent differences 
in spectral reflectivity. Interchange of the channel or ratio with the 
colors selected is also possible by the user to maxiniize color contrasty 
and improve interpretability. 
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2. 2.2. 3 Methodology Study 

After the desired functions and flexibility had been designed 
into the image enhancement and image generating program, a series of 
images were produced which covered areas under investigation by the 
SRSL. These studies covered (s-) the ERTS-ground correlation site in 
the hills adjacent to Stanford University (Stanford Grassland Sitet), 

(b) serpentines adjacent to Crystal Springs Reservoir and the San 
Andreas Fault in central California, (c) the Yerington, Nevada open 
pit copper mine and (d) a geobotanical anomaly related to a molybdenum^ 
rich area in the Pine Nut Mountains, Nevada. 

The development of the methodology evolved, as images were generated 
relating to these studies. Detailed analysis of the enhanced images 
and the results of these studies are covered in Section IV of this report 
A discussion of these images as they relate to the methodology and the 
.application of the program in producing computer enhanced images follows. 

2. 2. 2. 4 IMAGE ANALYSIS 

Initially a preliminary study of the area of interest was made, 
viith STANSOET, to pin point the location and assess the problem. The 
first image generated on the full 4X5 inch polaroid format is geo- 
metrically corrected but unenhanced (see Figure 1, 6, and 11) . 
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With all our images the full 4X5 inch Polaroid format was 
utilized (300 by 250 pixels of ERTS data) which produces about a 
1:210,000 scale. Histograms are then generated for each of the 4 
channels and selected ratios that are believed would contain informa- 
tion of significance. By studying these liistograms one may decide 
upon the DN levels which it would be desireable to use, to stretch by 
a DN subtraction, with an amplification of the remaining levels. Figures 
1 thru 20 indicate various combinations that were used for enhancement 
and show the images generated. During this preliminary phase many 
more images were generated then req^uired for analysis in order to gain 
a more complete understanding of the capabilities and limitations of the 
system. In general our review of the images shox^s an overall improvement 

in contrast with a DN subtraction and amplification (stretching) ." A 

comparison of Figure 1 and 2, 6 and 8, 11 and 12, illustrate this 
improvement. However, it should also be noted that an increase in 
amplification often will eliminate information by boosting data beyond 
the 255 DN range. This can be seen iin Figures 9D, 13 A,B,C and 18 B,C. 

Of course, this is not universally true through the entire image or for 
every channel. In fact, some areas appear to be optimized; therefore, 
the specifics of a problem must be considered and often the aesthetics 
are sacrificed to the enhancement and specific information content 
obtained. 

Review of the ratioed images where by each channel is divided by 
Gh 4 and particularly Figures 10, 14 and 19 illustrate the minimizing of 


radiance differences xdiich appears to flatten the topography. With 
the DN subtraction and amplification the contrast improvement is 
again evident (compare figures 14 and 15; and 19 and 20) . The 
possible loss of information content with increased amplification 
(stretching) is show in figures 16C and 21G. It should be noted 
that the images presented here only show integer amplification 
values. The program was subsequently modified to include the greater 
f^xdxibility of floating point values. After generating the computer 
enhanced tapes for the area of interest the best combinations of 
stx'aight channel data or ratios were selected for the generation of 
the color composites which are then available for further study. 

The colours assigned to the various channels or ratios can be varied 
to suit the investigator and possibly improve image analysis. A 
cursory review of the color composites presented as figures 22 through 
23 indicate the improvement in interpretability by the superposition of the 
various bands or ratios and the addition of the color dimension. 

It is obvious from the above that continued optiiuiziation of 
the enhancement process will occur as additional type areas are studied 
and more understanding is gained of the relationshp betoeen terrain types 
and spectral behavior. 



2,2*2,4,1 Stanford-l^ASA/Ames, (1075/18154; October 6 1972 

image) 

In this series of images, Tigures 1“5 a progression of DIT subtraction 
and amplifications were applied to individual EE.TS channels as well as 
ratioed channels in order to optimize or enhance the information content „ 

In Figure 1, it can ‘he seen that the image generated from the raw EFTS 
data is generally of low contras t(^particularly channels 4 and 5, with 5 
somewhat better than 4) # Features such as the turbidity of the bay, the 
wooded areas as well as watered (bright) grassy areas (golf courses, etc), 
runways, main highways are evident hut require careful scrutiny to distinguish 
then from the background. In channel 7 the contrast is better with the 
reservoirs, hay, bayland ponds and road networks quite evident; however, the 
topogi‘aphy of the hills is poorly defined* The turbidity of the shallower 
creeks and the inshore areas of the bay can also he seen by contrast with the 
darker adjacent salt ponds* 

The DN subtraction and amplification demonstrated in Figures 2 and 3 
improve the overall contrast with the turbidity features of the hay and its 
shores more sharply defined* Felt Lake in the lower central area of the image 
is now very evident in channel 5; however, close examination of Figure 3 indicates 
that the clipping of the DH level has effected the outline (reduced) by 
eliminating the damp transitional areas of the shoreline* The wooded areas on 
the left and bottom side of the images are now easier to see in channel 5 as 
well as the road network and watered grassy areas of channel 7. 

Figures 4 and 5 contain images obtained by ratioing channels. Study of 
these images indicate both a gain and a loss of information content as a result 
of the ratioing and enhancement procedures. For instance, in Figure 5, at the 
higher DN clipping and amplification levels a peculiar salt and pepper effect 
tends to blend features in the 5/4 ratio. Over saturation (over amplification) 
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of the signal in the 6/4 ratio ‘while emphasiziTig the ' *anford Linear 
Accelerator because of the background brightness of the lower left quadrant 
of the image tends to obliterate other detail in that area. In general, 
the topographic relief is subdued in all images of Figure 5, which tends 
to emphasise the cultural features such as the developed areas and 
road networks* The DN and amplification levels of Figure 4 seem to he 
more optimum, as the salt and pepper effect is not apparent, the topographic 
detail is more subdued hut the highlights are more apparent without appearing 
washed out. The road networks and cultural features also seem optimized 
as well as the bay and bay land features. 
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2. 2. 2. 4. 2 Crystal Springs Reservoir, California 
1075rl8154, Oiiober 6 , 1972 images 

Tn this study, an investigation was made relative to the possibility 
of differentiating, serpentine outcrops and soils to the east of Crystal 
Springs Reservoir, the surface manifestation of the San Andreas fault in 
this area. In this series of images Figures 6-10, the effect of increasing 
amplification with a fixed (optimum) DN subtraction on specific information 
content can be seen. In Figure 6, the raw data was utilized to generate 
the images. As before, the contrast level in generally poor in channel 4 and 
5. It is in channel 6 and 7 that the serpentine area, east of Crystal Springs 
and south of the Crystal Springs Golf Course (bright area in center of image) 
can be detected. The initial DN subtraction shown in Figure 7 has the overall 
effect of darkening the overall aspect of the images with some detectable 
contrast increase in channels 4 and 5. In Figure 8 the amplification increase 
to 2 improves the contrast and tends to isolate the serpentine area somewhat, 
particularly in channel 7. An increase in amplification to 3 in Figure 9, 
further isolates the serpentine in channel 7 but the over saturation reduces 
the Information content elsewhere by the general increase in brightness level. 
With progressive amplification it is possible to emphasize the difference in 
character of the topography and ground cover east and west of the fault zone. 

This is seen best in channel 5 of Figure 8. Also evident with this increase 
in amplification is the Increase in contrast of the watered grassy (bright) areas 
(golf courses, cemeteries, parks etc) in channel 6 and 7, the bay turbidit^^ 
in channel 5, and the road networks in channel 7. The shore outlines of the bay 
and Crystal Springs are sharpened appreciably as the contrast is improved. 

Ratioed Images are presented in Figure 10. It is evident that while the 
display of the serpentine area is optimized in the 7/4 ratio further DN 
subtraction and amplification would only tend to eliminate information content 
elsewhere. 
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2. 2.2. 4. 3 Yerington Pit, Nevada, 1397-18154, August 24, 1973 
2, 2. 2. 4. 3.1 Bands 4, 5, 6 and 7 

Figures 11-16 represent images made from ERTS tape 1397-18051, from 
August 24, 1973. This tape was selected for study because it was made by 
ERTS only 13 days after a RB57 underflight (for SKYLAB SL3) obtained excellent 
photographic coverage with B/W and color films. In addition the high sun 
elevation (67*^) of this date aids the spectral content. 

All four bands have been processed on each figure. In sequence Figure 11 
represents the 4 raw data images. Figure 12 shows them with suitable digital 
values (DN) removed, but still with amplifications of 1.0. In Figure 13 
the amplification has been raised to 3.0, stretching the data linearly. 

(By this point Channel 5 (-45DN; amp 3) is starting to show a ragged appearance, 
due to noise. 

At first glance there is liftle difference between any of the four images, 
except in the obvious vegetation at the right hand margin, of the irrigated 
crops (alfalfa etc.) of the flat, agriculturally-rich Mason Valley, in 
which the town of Yerington lies off the image to the lower right. A few 
small areas of vegetation may be seen similarly north of the shorelines of the 
mostly dry, Artesia Lake playa. Elsewhere in the Singatse Range, centrally 
running N-S through the image, no obvious spectral differences exist. The 
banana-shaped patch at the central right best seen on Channel 4 and 5, represents 
the gardens and housing of the mine. 

The black and darker grey triangles in the upper right are noticeably 
different in Channels 4, 5 a.*d 6. The area darkest in Channel 4 is acidified 
ferric sulfate leach liquor, used to extract *^oxide** copper from the adjacent 
piles of mineralized rock. It has a deep rust-red color to the eye and here 
is black even in Channel 4, regardless of its depth. The cigar shaped patch, 
best seen in Channel 6 and 7 is a shallow pond of water at the lowr end (northern) 
of the tailings pond. This pond becomes essentially transparent in Channel 4. 

Clipping suitable digital values (DN) off the data sets seems to bring 
the values into a better portion of the grey scale of the (Polaroid Type 107) 
film. It does somewhat increase the information content, but not their spectral 
differences. Further amplification to 3.0, with the same DN offsets, serves to 
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brighten the whole scene, although now the lightest areas are hopelessly 
overexposed . 

A false color (CIR) image was made using the following settings of the 
3 best Images, and colored using the filters as indicated; 

Ch 7 (-17DN), amp 2 RED 

Ch 6 (-42DN), amp 3 GREEN 

Ch 5 (-45DN), amp 3 BLUE 

A black and white representation of the original color print is shown 
in Figure 22A 

2. 2. 2. 4. 3. 2 Ratio of Bands (R5/4, R6/4, R7/4) 

Comparisons between ratio images (Figures 14, 15 and 16), with 
the raw data sets (Figures 11-13) show immediately that most of the topographic 
effects have been removed, providing a ”flat image**. Figures 14, 15 and 16 
represent successive steps in DN subtraction and increasing **stretch’* 
(amplification) . 

The ratio images now differ slightly from each other (which makes the 
color print show differences) . Of these one now sees a darker patch in the 
lower right center now appearing which correlates with the granite (Kg) 
on the geological map and is terminated by a roughly E-W line (fault). The 
best simple ratio for lithological purposes is Ch 6/4. 

Figure 22 (even in its B and W form in this report) shows these contrasts 
clearly, particularly those which result from the beneficial removal of 
topographic effects and the heightened degree of the rock type discrimination. 
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2. 2. 2.4, 4 Pine Nut Mountains, Nevada (Molybdenum-vegetation anomaly) 
1289-18063, Ray 8, 1973 


A similar series of images for Channels 4, 6 and 7 appear in Figures 
17 and 18 with DN subtraction and increasing amplification (stretch). A 
large area of snow, taking up most of the right hand (E) edge of the image 
dominates the scene. Because the ERTS system clips all data values above 
127DN the snow brightness variability (now reduced evenly to 127DN counts) 
causes problems in subsequently ratioing steps. (In addition while making 
these particular images the NASA/ARC DICOMED unit was suffering ‘’electronic- 
overshoot** problems resulting in blacks anomalously appearing in the snow) . 

Notice also that the areal extent of the ’’white** snow is larger on Ch 4 > Ch 5 > Ch 7, 
in accord with the known spectral pattern for late spring snow melting. 


Several localities are keyed to symbols on the figures. 


Figure 17 (a) — Mo/vegetation lies (1 cm) to left (w) of this 
letter (Figure 17c). 

(b) — Sugar Loaf Hill, andesitic plug intrusive 

(c) — Suspected caldera (outlined by arrow points) 

of Double Springs Flat traversed by Highway 395. 


Figure 18 (a) — Small farming area E-W patch (Figure 18b) 

Figure 19 (b) — Trace of 3-5 year old fire, which burned off the pines 
and junipers, now replaced bv sagebrush. 


Figure 23 shows the false color (CIR) composite print (here is 
black and white) of these images. 


Rationing (in areas outside the snow problem) again removes the topography 
allowing the observer to concentrate on spectral differences. Figure 19 has 
zero DN offset. Figure 20 has the best DN offset, and Figure 21 increases the 
stretch (amplification) from 1.0 to 2.0. Figure 23 (B/W copy) shows the effect 
of composite-color-ratios. 

The best simple ratio is again Ch 6/4, in which the (sparse) broad-leaved 
vegetation (near Spring of Wales) shows white patches. The spotty white/black 
pixel in the snow are due to the ’’over shoot” problem. 
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Wati^r in Double Springs Flat shov%^s black spectrally in the 7/4 and 6/4 
ratios (lower left corner). The firo trace (b) does not appear in Ch 5/4 
but is emphasized in Ch 6/4 and Ch 7/4. 

Channel 5/4 ratio however does show patterns in the geology in the 
upper right corner correlatable with the published geological maps. 

The Mo-vegetation anomaly is not easily seen but may be discerned about 
1 cm to the left of (b) in Figure 21B, as a grey rounded patch. 
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2 . 2 . 2 - 5 Conclusions 


The rapid interactive and iterative nature of the program 
leads to the continuing development of applications technique as 
experience is , gained in an expeditious fashion* An integration 
of the SRSL STANSOE.T program with the image enhancemer ‘ program 
(IMAGE) would improve the efficiency of the process* Although, at 
the expense of lowering the overall speed drastically. Also recommended 
is the use of an on lines real time TV -viewing system so that the 
computer enhancement may be viewed, modified as desired, and finalised 
images produced more rapidly* 


2,3 INTERPRETATION OF LANDSAT DIGITAL TAPES 


2,3,1 ATMOSPHERIC EFFECTS 

With multiple coverages available from LANDSAT-1^ it has been 
possible to effect a calibration procedurej using an approximately 45 acre 
target, in the northern part of San Francisco Bay,. This target is composed 
of carbon black, a waste product of oil refining nearby, and was found by 
searching Channel 5 and 7 iamges for very black targets — this target was 
very black inboth. 

The fallowing is an abstract of a published paper (appearing in 
Appendix A), detailing the use of black targets for atmospheric correction, 
The simultaneous use of a "bright” target of known bidirectional reflectance 
(using the 4 LANDSAT filter bandpasses) can transform the radiance data 
sets to reflectance data sets. This has obvious use in comparing ground- 
measured reflectances, and in searching for matching spectra elsewhere 
in a LANDSAT tape. 
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2/3. 1.1 Abstract of Published Paper 
(in Appendix A) 

A COMPARISON OF OBSERVED AM) MODEL-PREDICTED ATMOSPHERIC PERTURBATIONS 
ON TARGET RADIANCES MEASURED EY ERTS 


R.J.P. Lyon 
F. R. Honey 

Stanford Remote Sensing Lab 
Stanford University 
Stanford, California 94305 

SUMMARY 

In order to be able to compare results from ERTS MSS data over a 
series of tapes, the perturbing effects of a variable contribution due 
to radiation scattered by the atmosphere into the detector field of 
view, and of the variation in the irradiance on a target with solar 
zenith angle, must be eliminated. These two effects may be compensated 
for, or entirely removed, by studying selected targets in a scene, 
one (or more) of low (zero) reflectance, one (or more) or high, knoim 
reflectance. In some scenes, however, suitable reflectance targets 
may not be obtained. T-Jhen this occurs, atmospheric modelling must he 
employed to arrive at some values for the atmospheric scattering 
contribution, and for the irradiance on the scene. 

Two targets of measured, constant reflectivity in the area of San 
Francisco, California are studied. The first standard, a waste 
products treatment pond at an oil refinery near Suisan Bay, having an 
area of approximately 0,3 square miles, and bandpass reflectances of 
<0.5% in all Tour bands, is assumed to have a zero contribution to the 
radiance recorded by ERTS. The radiance observed then arises entirely 
from atmospheric scattering. The variation in these radiance values as 
a function of solar zenith angle Is compared x«7ith models for atmospheric 
scattering* 

A second target, a concrete parking apvxn for aircraft at Moffett 
Field, California, assuming that it remains dry during the period of study 
has constant reflectances of 27.8, 31.0, 30.0, and 32.3 percent bandpass 
reflectances in four MSS equivalent channels* Using these values, the 
radiance observed by ERTS may be corrected for the atmospheric contribution, 
and thus values for the irradiance on the target may be calculated. Thes 
values may be studied as a function of solar zenith angle and compared ^nlth 
results from models , 

The technique of using standard targets within a scene is applied to 
a specific scene which contains an area oi measured reflectivity. 
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.3.2 GEOLOGY AND SOILS OF THE STANFORD G[IASSLANDS SITE 


n 


2, 3 » 2.1 Introduction 

The foothills behind the Stanford campus are exemplary of the rolling 
topography encountered the length of tbe San Francisco Peninsula* The 
rocks are predominately of marine origin and have undergone considerable 
deformation since their lithification. Local relief is substantial with 
hills rising to elevations above 5 OO feet. Topography is rolling ac a 
result of structural control and a well developed soil profile* 

2,3*2. 2 The Project Area 

The area is characterized by a rolling hill and swale topography 
indicative of the climate and geological structure of the region, Vege- 
tation is predoTainantly annual grasses; a few large oaks dominate the 
grassy knolls , Slopewash deposits and soils are thick enough to conceal 
most of the bedrock. Outcrops occur on steep slopes where erosion is able 
to remove surficial material faster than Its production. Faults express 
themselves in minor control of topography. The 6il*“year rainfall average 
is 15 inches on the east (lower) side and 2Q inches on the western (higher) 
side. Over 80^ of this falls between December and March. Streams occupy 
deeply incised valleys* Rejuvenated headward erosion in these valleys 
indicates a change in runoff conditions. This is interpreted as the result 
of overgrazing by cattle. 
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2* 3. 2.3 Geology, Rock Units 


The rock types encountered in the project are marine in origin. During 
part of the early Tertiary^ these rocks must have accumulated in an off shore 
environment. A schematic statigraphic section is included in the appendices. 

1 Butano Formation (Eocene^ Tbu) . 

This unit is composed of medium grained sandstone beds averaging 60 
feet in thickness and dark grey-green shales averaging 60 feet in section. 

The sandstone is massive^ fragments angular, poorly sorted, and arkosic. 
Relatively thin beds (1-2 ft.) of grey silt and mudstone are included in . 
these larger beds. Outcrops weather to a yellow buff color and are much 
harder than unweathered rock. 

The shale beds of this unit are greenish to brown grey and contain a 
predominate proportion of clay minerals. Atchley and Grose (I 960 ) state 
that core specimens contained as much as me ^tmorillinite. Those which 
appeared dry would shrink upon exposure to the atmosphere. The resultant 
cracks represented a minimum of of the original volume* 

Edwards (I96I) examined microfossils from the mudstones of the unit 
and placed it in the Eocene, 

The unit is probably more than 2000 feet in section, Alon^ 

Francisquito Creek at Searsville Lake there are more than hOOO ft 
vertical section exposed* 
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Detailed Discussion: 


The unit is of marine sedimentary origin and is dated by fossil con- 
tent as Eocene age. It presently occupies the core of a faulted anticlinal 
fold and is overlain with angular unconformity by younger Miocene rocks. 

The unit has undergone at least two major episodes of burial ^ deformation^ 
uplift j and erosion. The component sediments j originally sand and clay mud 
were lithified to sandstone and clay shale and subsequently deformed and 
distorted. 

The lithology of the unit is characterized by irregular alternation of 
sandstone beds and clay shale beds. The individual beds vary in thickness 
from 1 foot to more than 100 feet^ and within any given section the ratio 
of sandstone to shale may vary from 70:50 50:70* 

The Eocene sandstones typically are fine- to -medium grained^ poorly 
sorted^ and variably cemented, Beds of mappable continuity are commonly 
50 to 90 feet thick 'and often contain 1 to 2 foot thick interbeds of shale. 

t 

The contacts are Irregular and the attitudes of the beds may vary abruptly. 
Individual beds of sandstone or shale may pinch or swell in thickness, 
terminate abruptly, or coalesce with adjacent beds . 

The sandstones consistently are light grey in color when fresh; they 
generally weather to white or light brown with, conspicuous iron stains on 
fracture and bedding planes. Outcrops generally are "case-hardened" and 
are often harder tban fresh rock. Road cuts along Junipero Serra Blvd. 
show criss-cross joints which, yield 2 to 3 foot blocks in the weathered, 
case-hardened rock. The Eocene sandstones tend to be somewhat harder than 
the Miocene sandstones . 
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1 ,\ 

Tha Eocene shales consist perdominantly of clay rocks ^ but include 
minor sandy clays ^ silty clays and marls. Most of the clay rocks have 
the appearance of mudstone^ but thinly laminated shale Is also found- 
The degree of distortion in the shales varies from severe shearing to 
blocky fracturing. The common distortion of the clay beds appears to 
be primarily the result of. wide- spread, intr a -bed movement during deforma- 
tion and folding and only in small part the result of shearing associated 
^Ith fault movement * 

Mappable shale beds range in thickness from 20 feet to rarely over 
loo feet. Individual beds usually contain a mixture of several varieties 
of shale and commonly 2 to 3 foot interbeds of sandstone. In several core 
intervals, tectonic mixing of sandstone and clay was noted. 

Several varieties of clay shale were recognized in the drill cores but 
none of these could be correlated with the weathered varieties exposed in 
the trenches. The dominant variety observed in the drill cores is a blocky, 
hard, dry., medium-dark olive-grey or brownish-grey clay with occasional thin 
horizons of sheared, light green-grey clay. This light grey clay was analyzed 
by X-ray diffraction and found to contain approximately 50^ montmorillonite. 
Some^^hat less abundant is a moderately stiff, hard, dark, greenish-grey 
chioritic clay shale. 

Most of the shales appear "dry*^ when first cored, but on exposure to 
air, they further dry out and develop extensive shrinkage cracks amounting 
to as much as 5^ of the original volume. Many of the clay horizons contain 
variable amounts of authigehic pyrite which, in the weathered zone, oxidizes 
to form abundant secondary gypsum. Practically all of the clay shales are 
abundantly fossilif erous . 
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2.3. 2.3. 2. Page Mill Basalt (Miocene. Tpb) 

'<■' This unit is a series of separate volcanic flows. Detailed study by 

Atchley and Grose (195^0 sequence of flows and mapped their 

locations. The basalt of the 9 flows they identify have three distinctly 
'C' different textures. (See Map II and Appendix II.) 

The massive basalt is fine grained and very hard. However, it is 
extensively fractured and columnar jointed. The vesicular flow is tuff- 
& aceous, massive, and well statified by horizons of gritty pumiceous debris. 

The breccia has fragments of various sizes incorporated into it. Outcrops 
weather to a distinctive bro^ra. 

d Cummings (1956V concluded that this unit was contemporaneous with the 

Mendigo diabase and tuff exposed in Woodside. 

Poland (1959) recognized the submarine environment of deposition of 
d the basalt without finding chilled pillows. 

The unit varies in thickness. It is as much as 400 feet thick in the 
Quarry off of Page Mill Road. 

Detailed Discussion; 

The Miocene volcanics include three principal rock types: fine- 

grained flow basalt, blocky volcanic agglomerate, and local volcanic 
tuff and tuffaceous sandstone. Previous mapping distinguished the basalt 
flows from the agglomerate. The present map further refines the distribu- 
tion of the basalt flows; the tuffaceous rocks were mapped as part of the 
agglomerate unit. The volcanic agglomerate occurs sporadically intermixed 
with basalt flows in the Page Mill quarry, and as a thin capping on the 
slopes and hilltops immediately west of Page Mill Road. 
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The agglomerate is a composite mixture of fragmental volcanic debris 
imbedded in a matrix composed principally of hardened ash and mud^ It 
appears to be a mixture of volcanic fragment mud flow sediment^ and 
landslide debris ^ such as is often found on the flanks, of active vol- 
canoes , It contains angular to subanguiar fragments of volcanic rock^ 
usually vesicular basalt^, which range in size from one to several inches 
across. The agglomerate is a soft weak rock and its dominant occurrence 
in the. Page Mill quarry gives an erroneous impression of conditions in the 
volcanic rock sections. 

The basalt flows j .four of which are recognized, are well exposed in 
the Page Mill quarry. Along the west limb of the anticline^ they lens 
into a single flow which thins and swells in thickness along a narrow band 
extending beyond Alpine Road, T^here penetrated by tunnels^ the basalt unit 
consists, of a single flow_, complete with vesicular flow surface and baked 
lower contact. The rock is an extremely hard^ dense, fine-grained basalt 
but is highly fractured, jointed, and transected by hair-line cracks which 
result in easy breakage. The basalt flows typically contain abundant pyrite 
disseminated \7i thin the rock or filling cross-cutting fractures. There is 
considerable chloritic-type alteration along the fractures associated with 
the pyrite mineralization. 

The tuffaceous rocks consist of massive, well- stratified horizons of 
gritty pumiceous debris mixed locally with coarse sandy materials , Much 
of this rock appears to be water-deposited, but some of it is clearly of 
volcanic origin. The rock Itself is variable in hardness, density, color 
and composition, and generally forms conspicuous brown- stained outcrops. 
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2, 3. 2* 3*3. Lower and Middle Miocene Sandstone (Un-naTned; Ts ) 

The lower part of this imit is. richly, fossiliferous and rests uncon- 
formably on the volcanic unit. Fragments of marina invertebrate shells 
(such as barnacles, and pelecypods ) are .abundant . Calcareous cement mal^ies 
this course grained j well sorted^feldspathic sandstone more resistant to 
erosion than the. upper part of the unit * Hogg (I 9 S 5 ) 

presence of megafossils indicated a shallow^ envixonment of deposition. 
Weathered, rock is case hardened. 

Higher in the section the sands are finer and moderately .well sorted* 

. . The rocks, are yellow-grey to olive, green and quite friable^ lacking the 
calcareous cement seen lower in section. Calcareous concretions are large. 
Shell fragments form resistant beds* . This rock VJeathers to a light, grey. 
Outcrops are rare and slop ewash cover, thick, Antonnen (I966) concluded 
that the sorting and feld spathic nature of the sandstone indicated a. rela- 
tively immature sand being deposited near the source. 

Los Trancos Formation has been proposed as a name^fdr this unit but 
remains unadopted. 

Detailed Discussion; 


The unit is composed predominately of massive^ thick-bedded^ fine- 
grained to silty sandstone which is poorly to moderately cemented. Some 
of the rock Is very poorly cemented^ almost a: firmly compacted sand. 
Inter stratified are local 1 to beds and lenses of coarse-grained 

sandstone which also varies in hardness and cementation. Clay -.or. shale., 
strata are virtually absent. The moderately cemented sandstones consti- 
tute .the bulk, of the rock^ but even these are: relatively soft and cores 
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can Le broken by hand or with light hammer blows. Much of the softer rock 
can be crushed by finger pressure. Only a small percentage of the sand- 
stones can be classed as hard to very hard^ such that a 2“inch core would 
require several haimner blows to fragment* The harder rocks include thin 
beds of calcite-cemented fossiliferous sandstone scattered through the 
lower portion of the unit^ and occasional boulder-size calcareous . con- 
cretions - 

The fine-grained sandstones are. medium grey in color^ weathering to 
light buff or light brown^ and are quite friable. Most of . the sands are 
fairly well sorted^ with generally subangular grains^ and are only 
slightly permeable. Many of the sandstones are ^*arkosic*^ or .'"feldspathic” 
as they contain a variable^ though high^ percentage of feldspar minerals* 
Distinctive features of the Miocene sandstones are their, massive character^ 
the general scarcity of shale horizons^ and the. presence of . thin fossil- ■ 
fragment beds* Near the western part of the area the sandstones are inter- 
bedded with the Miocene volcanic sequence and are found beneath the basalt 
flows * 


2 * 3 - 2 . 3 . 4 . 


Upper Miocene Silicious Siltstone 


This unit is represented at the' southern boundary of the project area. 
The lower sandstones grade into it conformably. The silicious cement indi- 
cates a much deeper, calcium deficient environment of deposition. The chert 
is very fractured/ 'Weathering then produces fragmented debris* 



‘fl 

2 . 3 . 2 . 3 . 5 . 

The Santa Clara formation^ of widespread occiirfence along the inargins 
of the Santa Clara Valley, is coramonly described as a terrestrial deposit 
of Plio-Pleistocene age. In the Stanford Foothills, this formation is 
partly of marine origin and probably per dominantly of early Pleistocene 
age. Along Arastadero Road south of Felt Lake^ nearly flat lying beds of 
; sands and gravels of this formation are well exposed. These beds are near 
the base of the formation and have fossil horizons containing a marine 
pele.cppod. fauna * Similar beds are. found in the hills south of the inter- 
section of Page Mill and Arastadero Roads* This fauna has been Identified 
by Dr. Keen (personal communication^, July^ 1959) being no older than 
Pleistocene* Thus^ in the central part of the mapped area^ the Santa Clara 
rocks would appear to be essentially Pleistocene in age though rocks of 
: pliocene age may occur furtherto the southeast* 

The Santa Clara formation consists perdominantly of poorly bedded 
unconsolidated sandstone with layers and lenses of interbedded sandy gravels , 

The gravel is mostly , of small size^ from l/k*^ to hut larger sizes: are 

i-.-r 

not uncommon* The rocks are generally similar to local deposits of recent 
alluvium although there are some obvious differences in the source areas of 
the gravels. The Santa Clara conmionly shows evidence of accumulation from 
local sources such as the Eocene sandstone and the Miocene basalt. 

The Santa Clara rocks rest unconformably on all of the older formations 
^ and at one time may have covered most of the area. Remnants of gravel are 

found almost everywhere even though the bulk of the formation has been 
eroded from most of the hill areas. The Santa Clara rocks are gently folded 
in the central part of the area^ but they have been greatly .affected by 
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faulting along the San Andreas sone and have been extensively folded and . 
faulted along with the underlying Nlocena rocks ^ near the border of the 
alluvial plain^ northeast of Junipero Serra Boulevard* 

Some late Pleistocene and Recent terrace gravels ware mapped with the 
Santa Clara rocks since^ in the absence of good exposures^ the two are 
practically indistinguishable- The thickness of the Santa Clara rocks in 
the mapped area is not known; however the formation is several thousand 
feet thick in the region iinmediately southwest of the San Andreas fault. 

2* 3. 2 -3. 6 Recent Alluvium 

Recent alluvx^ is present in the foothill area only along the stream 
valleys; it forms no deep valley fillings except along the. San Andreas rift 
zone* Extensive thick deposits of alluvium are present along the northeast 
margin of the foothills • 


2. 3. 2. 4 Structure 

The rocks described above are all involved in local distortion result- 
ing from folding during the Late Miocene orogeny which produced the coast 
ranges . 

The project area is dominated by an anticlinal fold which has been 
breached by erosion. 

Atchley and Dobbs (I960) point out that the uomplexity of distortions 
obscures the structure. However^ one can recognize the dip of certain beds 
in the Butano Formation and the dip of the Lower Miocene Sandstones . 

Reverse faults displace parts of the volcanic unit and involve the 
Lower Middle Miocene Sandstones as well. This structural irregularity is 
expressed in valleys cutting across the inclined beds of the Butano 
Formation. 



2* 3* 2^5 Geometric Considerations of ERTS Imagery 

The common approach for dealing with areally distributed (two- 
diitiensional). data is to prepare an appropriately scaled base (usually 
orthographic topography) *'^ou which different data, types can be plotted. 

The ERTS. scanning ,system presents certain technical problems to this 
approach. The actual ground area for which spectral reflectance is 
sampled is a rectangle approximately I 87 by 259 feet East-^Test and North- 
South respectively* Because the orbit of the satellite is not exactly 
polar j these rectangles are oeiented approximately lO^E of North* In . 
addition/ the MSS collects six such lines of data simultaneously - Each 
successive scan takes place after the earth has rotated to the east 
slightly (60 ft). Hence^ everj^ six lines are offset to the west almost 
1/2 one such element. The images produced by Goddard are rectified to. 
compensate for these geometric, characteristics and hence j are almost 
completely orthogonal, finally^ during a subsequent overflight the 
picture element may be collected as much as ^Q^ hoth North and East of 
the pixel from the previous overflight . 

The Stanford ERTS -tape reading system presents the data on a simple 
line printer (10 characters, per inch' Horizontally and 6 characters verti- 
cally). The results is a skewed and stretched ”image*‘ with scales of 
1:18482.64 and 1 :22171 - 5 ^ order to locate particular sites 

within such ^'images" it is necessary to find recognizable features such 
as bodies of water or highways (best seen in band 7) SiXid interpolate 
between them. This difficulty is confounded when one wishes to collect 


the actual nutneric value for one pixel as numeric printouts, are stretched 
horitorLtally X5 order to legibly present double digits- And^ o£ course^ 
when images from different overflights are to be compared this. is. the final 
problem' of determiniTig just what piece of ground was samp led - 

Costly and time-consuming efforts would be necessary to develop a 
system to present, various shade-prints geometrically corrected and at 
appropriate scale* It would be even more difficult to present numeric 
data in a form compatible with an orthogonal base map* 

One alternative to converting the "imaged* to orthogonality is to 
convert the orthogonal base to the geometry of the printouts * Such 
skewed and stretched maps could be layed directly over printouts* 

•Another approach is to construct templates for collecting discrete 
data points from either type of data format* 


2. 3. 2, 6 Soil Sf?..'pling 

Once actual locations of sites were established in the field their 
positions were plotted on a topographic map with a scle of l:2i^-00. 
Distance measurements made to features appearing on the map (such as road 
intersections^ structures^ hill-tops^ etc,) made it possible to locate 
sites to within a few feet. 

Samples of soil were taken at each site. The percent of weight com- 
prised by water was determined by oven drying. Color of soil was deter- 
mined by comparison with the Geological Society of America koclc-Color 
Chart, 

The soil type and bed-rock for each site was determined on the basis 
of field reconnaissance hnd mapping^ consultation of literature on the 
area^ and physical examination of samples. The actual procedures used 
are described iti the Appendices. The individual sites x-rere plotted on a 
template for numeric printout. When the dam of Felt Lake and several 
clumps of trees are located on both printouts (band 7 shade prints are 
darker and numerics have lowest numbers at these locations). Statistical 
analysis of the data should reveal the relationships if any ^ which exist 
between geologic characteristics and spectral reflectance as measured by 
ERTS , 
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Detailed Techniques Used 


2, 3. 2. 6.1 Soil Samples 

A black of top soil was e:K:cavaE:ed with pick-haimner and shovel. 
Following triinniiTig a fist size sample was placed in a plastic bag with 
a label and sealed. 

2. 3-2. 6.2 Soil Moisture 

E, ' 

A small amount of each sample was cut from each blocks placed: in 
crucible, weighed, and dried in an oven at 120^0 for 8-10 hours, and 
weighed again while still hot. The difference was used to calculate 
the percent of original weight which the moisture represented. 

2 . 3 . 2 . 6 . 3 Color 

The dried samples were set next to the color chips in the Rock- 
Color Chart and the Munsell number of that which, was most like the .soil 
was recorded. 

2.3.2. 6. 4 Soil Type * 

Local “ Local variations of soil types result . from down slope move- . 
ment of surficial materials and the resulting differential distribution 
of clay size particles. The classification in this category was made on 
the basis of field examination as well a> textural study of samples. 

Soil Conservation .Service ~ The SCS Soils Report for Santa Clara. . 
County describe different soil types and shows their areal distribution 
on photomaps of scale 1:12000* , A portion of this mapping is shown in 
Appendix VIII. The descriptions of the classifications are given in 
Table I. 
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2. 3. 2. 6. 5 Bed Eoclc Units 

The geologic unit for each, site was determined by fieldwork as well 
as using existing maps. Descriptions are included in the section on 
Rock Units . 

2. 3. 2. 6. 6 Inclination and Slope Azimuth 

With site locati'ons plotted on the small-scale topographic map it 
was a simple matter to determine gross slope characteristics. A line 
was constructed perpendicular to adjacent contour lines. Its azimuth 
was noted. The distance between contour lines was then used to trig- 
onometrically calculate dip of surface below horizontal. These measure- 
ments appear in Table I. 

2, 3. 2. 6. 7 Elevation and Lambert Goordinates 

These were determined by simple interpolation on the grided topo- 
graphic map. These are noted in Table I as x^ell. 
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TABLE 2. 3.2.6. L 

STANFORD FOOTHILLS 
ERTS-A TEST SITES 

SOIL DRY % SOIL topographic LAMBERT CO-ORO 


y 

0 1 ... 

STAT 

descriptor 

COLOR 

MOISTURE 

SLOPE 

elev 

(ZONE 

111) 

f i • ■ ■ 

U 

NO, 

1 2 3 

HUE 

4 5 

11/15/74 

INCL AZMTH 

FT 

North 

EAST 

i 

0-029 

slr-gme/tbu 

* ■ 

/ 


14.0/I75SSE 

380 

1516262 

335158 


B-030 

FIL-(5ME/T8U . 

5YR 

5/1 

5.4 

4.7/243SWW 

360 

1516349 

335034 

■ J • ' 

8-0 31 

CSW-GME/TBU ' 

5YR 

4/1 

11. B 

I1.3/338NNW 

374 

1516425 

334905 


0-032 

CSW-GME/TeU 

5YR 

3/1 

10.1 

9.5/003NNE 

390 

1516451 

334801 


B-033 

CLR-LTD/TPB 

bYR 

4/1 

11.9 

5.7/0 X2NNE 

444 

1516454 

334502 


B-035 

CLR-LTD/TP8 

iOYR 

4/2 

10.7 

6.3/335NNW 

4b7 

1516476 

33429T 


C-041 

CLR-LTD/TPB 


/ 

« 

4.0/200SSW 

476 

1516596 

334078 

]J ■ 

C-042 

CSW-LTD/TPB 

lOYR 

2/2 

9,9 

3.2/201SSW 

460 

15X6569 

333951 


C-Q45 

CSW-LT0/TP8 

5Y 

4/1 

16.6 

2,3/X4SSEE 

466 

1516572 

333886 


C-046 

CSW-LTD/TP0 

SY 

3/1 

20.1 

2.8/089 E 

468 

1516507 

333847 

; ( 

C”047 

clr-ltd/tpb 

SYR 

3/2 

11.2 

4.1/095 Ek, 

.480 

1516321 

333739 


C-049 

CLR-LTO/TP8 

lOYR 

2/2 

17.2 

5,2/094 E 

493 

1516128 

333642 


D-OSl 

clr-ltd/tpb 

lOYR 

3/2 

10.9 

3.6/076 E 

512 

1515877 

333482 

. { 

0-053 

CLR-LTD/TP8 


/ 

■ '♦ 

4.B/074 E 

523 

1515690 

33329S 


D-054 

clr-ltd/tpb 


/ 

• # 

2.3/135SE 

516 

1515557 

333097 


D-0S5 

clr-ltd/tpb 

lOYR 

2/2 

13.9 

3.6/080 E 

SI9 

1515464 

332941 

: ( 

0—0 S6 

clr-ltd/tpb 

IOYR 

4/2 

9,0 

3.2/083 E 

528 

1515412 

332753 


E-060 

CSW-DAE/TaS 

5YR 

4/1 

17,0 

0. /320NW 

457 

1514286 

333305 


F“070 

slh-dae/tus 

IOYR 

4/2 

7.2 

7.1/028NNE 

4T5 

1514202 

333057 


F-072 

SLR-DAE/TUS 

IOYR 

3/2 

10,6 

14.Q/03QNNE 

496 

1514311 

332906 

! ( 

F-073 

SLR-DAE/TUS 

IOYR 

4/2 

7,2 

9,5/lliSWE 

486 

1514450 

332666 


G—080 

SLR-DAE/TUS 

5Y 

5/1 

8,3 

3.2/2B8SW 

483 

1614382 

332563 

1 

5-082 

slr-dae/tus 

5Y 

4/1 

15.1 

11.3/307NWW 

496 

1514296 

332406 

’ . r 

G-084 

slr-dae/tus 

IOYR 

4/2 

5,B 

9.5/237SWW 

430 

1514242 

332134 

1 ■ 

G-086 

SLR-DAE/TUS 

IOYR 

5/2 

6.6 

e,l/2S3SWW 

454 

1514197 

331815 


H-091 

SLR-OAE/tUS 

IOYR 

4/2 

7.9 

4.7/104SEE 

452 

1514096 

331612 

t ' ' 

H-092 

CSW-OAE/TUS 

IOYR 

3/2 

9.3 

. 8. /05S E 

445 

1514058 

331451 


H-094 

SLR-DAE/TUS 

5Y 

5/1 

10,1 

4, 1/061 NEE 

467 

1513964 

331246 


H— 096 

clr-dae/th 

5Y 

5/1 

11.2 

7a/0T4NEE 

477 

1513839 

331103 


I~005 

CLR-DAE/TM 

7YR 

4/1 

10.4 

0.5/190 S 

473 

1S13772 

330962 

1 j 

I-Q06 

clr-dae/tm 

5Y 

5/1 

10.4 

6.5/205SSW 

466 

1313719 

330 783 


1-008 

clr-dae/tm 

IOYR 

5/2 

8,7 . 

Q«S/HO*fr^SSW 

435 

1513610 

330562 


1-00 9 

CLR-AVO/QSC 

IOYR 

6/2 

8.3 

13, /207SSW 

405 

1513456 

330400 


1-010 

csw-avd/usc 

iOYP 

4/2 

9.3 

14, /235SWW 

36o 

1513318 

330322 
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ELEV 
FT 



\ 

STAT 

NO. 

SOIL ' DRY 

DESCRIPTOR i ■ COLOR 
1 2 3 hue 4 5 

.% SOIL topograph I 

MOISTURE slope 
11/15/74 INCL AZMTH 

LAMOEftT ( 
{ZONE 
NORTH 

CO-ORD 

III) 

EAST 



J-020 

CLR-AVD/QSC 

lOYR 

5/4 

7.9 

10. /130SE 

340 

1513606 

329426 

, - 

f 

0-0 1^? 

clr-avd/qsc 

lOYR 

5/2 

6.9 

S.5/200SSW 

358 

1513400 

329510 



0-016 

ClR“AVD/OSC 

lOYR 

4/2 

6.9 

9.5/141SE 

3/0 

1513081 

329558 



J-0 14. 

clr-avd/usc 

lOYH 

5/4 

6.9 

4,5/194 S 

368 

1512779 

329482 


( 

K-030 

CSv<~AVO/.jSC 

lOYR 

5/2 

9.2 

5,7/123SE 

3bb 

1312606 

329331 



K-031 

CLR-AVO/USC 

lOYR 

5/4 

7. a 

3.8/190 S 

397 

1512327 

329278 



K-032 

clr-avd/osc 

lOYR 

4/2 

8,(S 

5. /044NE 

393 

1512066 

329231 

■. 

{ 

K-034 

ClR-AVD/OSC 

iOYR 

4/2 

0 *0 

1. /285NWW 

400 

1511816 

329245 



K-035 

CLR-AVD/QSC 

lOYK 

4/2 

7.2 

10, /284NWW 

372 

1511557 

329295 



M-052 

CLR-AVD/OSC 

5Y 

5/1 

7.1 

4.1/312NW 

404 

1510937 

328267 


f 

H-054 

clr-avd/qsc 

iOVR 

4/2 

8,2 

2.5/316NW 

423 

1511102 

328403 



N-OOi 

SCL-PRD/QAL 

ioyr 

5/4 

6 fli 0 

2.8/068NEE ' 

381 

1510230 

330275 ■ 

A 

( 

N-002 

ScL-PRD/QAL 


/ 

n 

1.0/267 W 

383 

1510136 

330160 



N-003 

SCL-PRD/QAU 

: IOYR 

7/4 

5.5 

5, /275 W 

382 

1510072 

329943 



N-004- 

SCL-PRD/OAL 

i 

/ 


1, /064NEE 

392 

1510173 

329624 


( 

N-006 

scl-prd/qal 

1 

{ 

/ 


6.3/132SEE 

386 

1510211 

329378 



N-OOfl 

SCL-PRO/QAL 

1 

/ 

.*■ 

3,3/095 E 

382 

1510180 

329083 



N-OlO 

Scl-prd/qal 

1 . . 
1 

/ 

• * 

2.0/085 E 

381 

1510177 

328760 


f 

N-012 

SCL-PRD/QAL 


/ 

« 

4,5/l05SEE 

381 

1SIQ150 

328515 

rv 

'— ( 

N-014 

SCL-PRO/QAL 

. f 

/ 


6.3/090 E 

383 

1510124 

328216 


1 

. M-016 

SCL-PRD/QAL 


/ 


0,5/323nNW 

384 

1510000 

32/914 
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TABLE 2, 3. 2. 6.1 (ConC’d) 


/O'] 


1 LOCAL SOIL TYPE - 

FIL- MAN PLACED FILL OF MIXED TEXTURES 
SLH= sandy LOAM residual 
CLR= CLAY LOAM RESIDUAL 
s£Wi= Sandy slgpewash 
CSW= CLAYEY SLOPEWASH 
; SCL= SANDY CLAY LOAM 

2 SOIL conservation SERVICE SOIL NAMES “ 

GM£= GAVIOTA-LOS GATOS COMPLEX 
GCE= GAVIOTA loam 
LTD= LOS TRANCOS STONY CLaY 
DAE“ DIABLO CLAY 
AVD=» AZULE SILTY CLAY 
PRD= POSITAS-SARATOGA LOAM 

3 GEOLOGIC BEDROCK NAME, - 

T8U= DUTANO formation (SANDSTONE ANO SHaLEJ 
TPB= PAGE MILL BASALT (MASSIVE® BRECCIA* AND TUFF FLOWS 
TQSss BARNICLE BED SANDSTONE 
TUS= UNNAMED SANDSTONE 
TM = MONTEREY shale (SILICIOUS) 

QSC= SANTA CLARA gravels 

qal= old and recent alluviUal deposits 


4 value (LIGHTNESS) 

5 CHROMA (SATURATION) 

MUNSELL number : NAME 


5Y 3/i = between »QLIVE BLACK» AND »0LIVE GRAY * 

5Y 4/1 = I OLIVE gray * 

5Y 5/1 BETWEEN ’LIGHT OLIVE GRAY’ AND ’OLIVE GRAY’ 
5YR 4/i = ’BROWNISH GRAY’ 

5YR S/1 = BETWEEN ’LIGHT BROWNISH GRAY’ 

AND ’BROWNISH GRAY* 

5YR 3/2 = ’GRAYISH BROWN* ' 

7YR 4/1 = BETWEEN ’LIGHT BROWN* AND ’YELLOWISH BROWN’ 
lOYR 2/2 = ’DUSKY YELLOWISH BROWN* 
lOYR 3/2 = ’GRAYISH YELLOWISH BROWN’ 
lOYR 4/2 - ’DARK YELLOWISH BROWN* 
lOYR S/2 = ’YELLOWISH BROWN* 
lOYR 6/2 = ’PALE YELLOWISH BROWN’ 
lOYR 5/4 = ’MODERATE YELLOWISH BROWN’ 
lOYR 7/4 « ’GRAYISH orange* 
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TABLE 2 . 3 . 2 . 6. 2 



GENERALIZED STRATIGRAPHIC COLUMN 


ROCK UNITS ENCOUNTERED IN THE DISH HILL AREA 


(no vertical scale) 


Uoper Miocene Silicious Siltstone 


fine grained, well statified, 
clay and mudstones with silicious 
cementation . 

chert concretions numerous 
dark grey weathers to light buff 


Middle Miocene Sandstone 


olive-grey arkos 
scarcely interbedded with light 
grey clay and shale, 
very friable 
few fossil beds 


Lower Miocene Sandstone 


fine grained, moderately well sorted, 
well cemented, fossi 1 i f erous , 
feldspathic v/ith basalt pebbles, 
weathers light brown and case hardened 


Lower Miocene Volcanics (Page Mill Basalt) 


series of flows ranging 


mas s i ve , hard, fractured, 
columnar jointing 

^) tuffaceous, massive, well sorte 
^ statified horizons of gritty 
pumiceous debris 

3). volcanic breccia , fragments of 
^ variable size 


Eocene Sandstone and Shal e (Butano Formation) 


fine to medium grained, poorly sorted, 
variably cemented sandstone, beds 30-90 ft 
interbedded with dark green shale, beds 
20-100 ft. 



composit afterj Atchley 
* a proposed name s 


TABLE 2. 3. 2.6. 3 


STATXGRAPHIC COT.UMN..OF ROCK UNITS 

, DISH HILL, STANFORD, 
(no vertical scale) 


Unnagerl Sandstone (Los Trances Formation*) 

700 feet 

fine fp^ained, moderately well sorted, 
wen cemented, f ossiliferops, feldspathic 
sandstone 

grading up into calcivim deficient, friable, 
grey to olive-grey arHos scarcely interbedded 
with light grey clay and shale near center 
of secti on 


Page Mill Basalt 300i"teet 

series of volcanic flovis ranging 

massive, extremely hard, extensively 
fractured and pointed 
tuffaceous, massive, well stratified 
horizons of gritty pumiceous debris 
■ 1 , breccia vjith fragsnentr of variable size 


Buteno Forme tl on 2000-fr feet 

fine to medium grained, poorly sorted, 
varinblly cemented sandstone beds 
(3p-Q0 ft) with daric grey clays tone seams 
(2-3 ft) composed of 50?5 montmorillinite 

interbedded with medium darlc olive-grey to 
brownish shale seamed with course to 
fine sand (1-2 ft) 

Sandstone v/ea there to case-hardened, 
huff outer pus 

Shale -renders a dar'k blaclc, clayrich soil 
Grose, (195'»') f.Antonan (1966) .Edwards (1961) 5^3 

.11 unadopted 99 
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2.3.3 THE VEGETATION OF THE STANFORD GRASSLAND SITE: A BIOMASS STUDY 

2. 3. 3.1 Description of the Study 

, A study of the vegetation at selected sites in the Stanford grassland 
has been undertaken to aid in the interpretation of reflectance data from 

those sites. 

The Stanford grassland is a typical representative of the California 
Valley Grassland plant community (Munz and Keck* 1965). It has been subjected 
to grazing by cattle for decades which has changed the species composition 
entirely. Few of the original native species remain. Most of the species 
of grasses and broad-leaved plants found in the grassland today have been 
introduced from the Mediterranean region (Thomas, 1961; McNaughton, 1968). 

A preliminary study was done to determine the species of the grasses 
and broad-leaved plants growing at the study sites. (Figs. 2. 3. 3.1 - 3). 
Plants were collected in early May in various stages of flower and seed 
formation. The plants were identified to the. level of genus or species 
using local floras. The nomenclature is that of Munz and Keck. Specimens 
of each 3 species were dried and pressed for a permanent reference collection. 
The major plant species found in the Stanford grassland study sites are: 

Botanical Name Common Name 


Grasses: Bromus mollis 

Avena barbata 
Loliura mult if lorum 
Bromus rigidus 
Hordeum lepor inum 
H ordeum hystrix 

A description and drawing of each 


Soft chess 5, 10 
Slender Wild Oats 4, 11 
Ryegrass 9, 12 
Ripgut grass 6, 13 
Foxtail 7 , 14 


Mediterranean barley 8, 15 
! species is attached (Table I) . 


Broad-leaved 
plants : 


Erodium sp . 

Geranium sp . 

Medicago sp . 
Convolvulus arVensis 
Bellardia trixago 


Filaree, needle plant 

Geranium 

Bur clover 

Morning glory, bindweed 
Bellardia 


Eschscholzia californica California poppy 


Rumex sp. 


Sorrel 


16 

17 

18 

19 

20 
21 
22 


Initial observations also revealed that the vegeta,tion at the study 
sites was variable in species composition, plant size., percent cover and 
time of onset of senescence and drying. 
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TABLE 2. 3,3.1 

Lescrlptiou of the ma'jor grass species 

Grass, species are. Identified by the characteristics of their flowers, 
and, to a lesser extent, their leaves. The flower head is borne on a culm , 
or stalk, which raises it above the leaves. The rachis , or maiti axis, of 
the head may he branched in a variety of patterns and hears the splkelets > 
the basic unit of grass flower structure. At the base of eAch Apikelet is 
a pair of glumes , modified leaf structures whose shape .and texture are 
important in species identification. The glumes subtend one to many florets , 
small modified flowers* The floret consists of an outer lemma and palea 
enclosing a small inner flower. Glumes, lemmas and paleas may terminate in 
slender bristles called awns * 

The grass leaf consists of a lower sheath , which surrounds the stem, 
and an upper blade which diverges from the stem. At the junction of the 
sheath and blade is the ligule , a hairy or membranous extension of the 
sheath. The margin of the leaf at the junction may be extended laterally 
auricles , or lobes. See Fig. 2, 3*3,1 - 2*3. 3.3 

Avena barbata Brot* Slender Wild Oat C^'ig. 2. 3*3. 4) 

Annual; culm 30-60 cia, tall: flower head open, loosely branched; 
spikelets large and drooping at maturity; florets 2; lemma hairy, bear- 
ing a prominant bent axra 3 cm* long. 

Bromus mollis L* Soft Chess; Soft Brome (Fig. 2. 3. 3. 5) 

Annual; culm 10-80 cm* tall; softly hairy all over plant; flower head 
compact, 4-10 cm. long, fex7 short branches; spikelets compact, slightly 
flattened laterally, 15-20 cm, long; florets 6-12 per spikelet; lemma 
. rounded on back, with soft awn 5-9 mm., long. . 

Bromus rigidis Roth. Ripgut Grass; Ripgut Brome (Fig. 2,3.3 *6) 

Annual; culm 30-70 cm. tall; flower head branched, open, 6-18 cm* 
long; spikelets 2-5 cia- long; florets 5-7 per spikelet; lemmas 2,5-3 
cm. long, tipped with a stiff awn 3.5-5 cm. long; both lemmas and awns 
covered with stiff short hairs pointing toward the tip. 

Hordeum hys trix Ro th . . Mediterranean Barley (Fig. 2,3-3*8) 

Annual; culms sometimes bent, 12-35 cm. tall, fcliage hairy; leaf blade 
lacks auricle; flower head 1.5-3 cm. long, unh ranched, with bottle brush 
appearance due to radiating awns; glumes rigid and divergent, laclcing 
hairs on margin; spikelets in 3* s; florets ! per spikela.t. 

Fig. 2. 3,3,1 from Cranipfcon (1974) 

Fig, 2/3. 3. 2 from Munz and Keck (1965) 

Fig. 2. 3-3^4 from Abrans (1940) 


103 


Hordeura leporlmim Link. ■ Foxtail (Fig. . 2. 3. .3. 7) 

• Annual; culm 15^60 cm. tall; leaf with auricle; flower head 5-9 cm, 
long, unbranched, with bottle brush appearance due to radiating 
awns; spikelets in 3’s; florets 1 per spikelet; glumes with hairs on . . 
margin. . . . 

Lolium multiflorum Lam. Ryegrass; Italian Ryegx-ass; . Australian Ryegrass (Fig. 2.3.3 

Annual, culm 25-lOQ cm- tall; blades with auricles at base; flower head 
10-20 cm. long , unbranched narrow spike; rypikelets flattened, close to 
and alternating on rachis; 10-20 florets per spikelet; ^.lumes shorter . 
tiian spikelet; lemmas with awns, 

2.3. 3,2 BIOMASS AND REFLECTANCES STUDIES 

From May 15 to May 22. a detailed study of ..the vegetation at 44 sites 
was made L‘o determine the species composition, biomass and stage of plant 
growth. At each site the vegetation was treated in the following manner: 

1. A square, 0*5 m, on a side, was marked off at a randomly 
selected site. 

2. Reflectance relative to BaSO. was recorded. 

..4 

3. The plant species were detenuined, 

4. The percent contribution of each species to the total 
biomass was estimated by eye. 

5. All vegetation within the square was cut off at ground 
level and put in an airtight plastic bag and taken to 
the lab. 

6 / Ref lectance after cutting was recorded, 

: 7, .Total fresh weight of the vegetation at each site was measured , . . 

8. All plant material was dried ill ovens at 100° - 5°C, for 48 hours, 

9. Total dry weight was measured. 

The reflectance at each site was measured before and after the removal, 
of the vegetation cover, using the. ERTS radiometer, 15°FGV. bidirectional, 
geometry (Tab. ' 2, 3, 3 *2) The sites vary considerably with regard to species 
composition, fresh vreight, dry weight (biomass) and the ratio, dry weight/ 

: fresh. weight,(Tab. 2, 3, 3,2 and 2,3.3*3) ^ The dry weight/fresh weight measure- ■ 
ment indicates the degree to which the plants have dried out. As the 
vegetation dries , the green color is lost and the leaves turn to yellow- ; 
green then tan. The ratio of dry. weight to fresh weight is therefore* an 
indirect measure of the ”greeness^' of the vegetation. These data will be 
used to interpret -the reflectance data taken at the- same sites, - 
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A statistical study of tlie correlations between species, biomass. 


biomass ratio and reflectance, appears, elsewhere as SRSL Technical Report 
No. 74-7. 

• . TABLE 2. 3.3. 2 


site Number 

STAMFORD GRASSLANDS 
Total wet weight (g) 

Total, dry weight Cg) 

Dry weight 
Wet weight 

946 , 

158.3 

78.9 

, 495 

942 

117.3 

52. S 

.450 

941 

105.0 

62.8 

.598 

947 (green) 

213.2 

120.0 

.563 

947 (dry) 

79.7 

57.3 

.719 

949 

119.3 

71.0 

.595 

951 

V ISO.Q 

. 81.9' 

.546 

953 

131.1 

82.7 

.631 

954 

^ : S3..9- 

55.5 

.662 

955 

70.1 

54.6 

,779 

980 

342.4 

138.4 

.404 

982 

421.8 

173. 0 

.410 

914 

282.5 

121.0 : 

.428 

916 

306.4 

128,2 

.418 

917 

184.6 

103.2 

.559 

920 

193.9 

104.6 

.539 

930 

333.7 

160.0 

.479 

931 

271.7 

155.0 

.570 

932 

406.5 

180.0 

.443 

934 

307.3 

130.1 

.423 

940 

344.8 

122,3 

.355 

942 

430.8 

x49.1 

.346 

943 

313.9 

138.5 

,441 

944 

290.0 

132.3 

.456 

905 

354.5 

168.2 

.474 

906, 

168,0 

50.1 : 

.476 

908 

209.5 

95.8 

.457 

909 

192.5 

113.6 

.590 

991 

220.4 

98.9 

.449 

992 

463.5 

. 175.7 

.379 

994 

404.3 

165.4 

.407 

996 

423,2 

197.2 

.466 

986 

411.0 

152.4 

.371 

984 

317.7 

126.4 

.398 

950 

325.5 ^ 

127.1 

.390 

952 

403.5 

145.3 

,360 

954 

367.1 . ■ • ■ / .. 

■ 164.2 : 

.447 

970 

86.5 

45.1 

.521 

S72 

152.0 

74.1 

.488 

973 

354.9 

171,4 

.483 


191.7 

81.8 • 

.427 

931 

294.6 

131.4 

.446 

933 

86.2 

48.2 

.559 

936 

60.3 

42.2 

.700 
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5 
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70 
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954 

.5 


95 





955 

40 
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980 

95 
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30 
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30 
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■■■5 . 
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994 . 1 
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10 





996 / 1 

1 10 

30 


60 
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30 


. 40 
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TABLE 2, 3. 3. 3 

Stanford Grassland - Species Composition 

















Fig. 2.3.3.15 Hordeum hvstrix Mediterranean barle 
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Fig, 2,3.3.21 Rumex sp 


2.3;4 statistical COEItELATION of BIOMASS DATA VERSUS BI-DIRECTIONAL 
, REFLECTANCE ; 

V- ; 2.3.4.I nSTRODPCTIOH - 

Whei^ the biomass data, from the 42 sites of the Stanford Grassland 
study became available, a statistical study was started by Prelat, 
using both the "scatter diagram" (2 variable plot) program (BMD02D) , 
and the stepwise discriminant program (BMD07M) . 

Selected examples of the output have been included in this section. 

2. 3.4.2 COMPUTATIOHAL METHOD (BMD02D) 

As stated by Sears in the preceding; Section 0.5m x 0.5m quadrants 
were selected at random near the measutement stations on the Grassland 
traverse over the period May 15-18, 1974. These stations were at in- 
tervals of about 80m along a 5 Km N-S traverse across the 4 main rock/ 
soil types in the Grassland site. 

For the biomass portion of the study, the grass was clipped off 
the quadrant, and weighed (wet and dried) in the lab, where species 
counts were also prepared. 

■: Immediately prior to the clipping, however, the bi-directional 

reflectance measuremehts were taken, relative to BaSO^, using an ERTS- 
bandpass t'diometer (EXOTECH Model 100). The measurements were repeated 
after cll,' ting , to give an appreciation of what "stubble" reflectance 
(grass roots + stem bases + soil) might lock like as seen from EKIS. 
2.3.4. 2.1 111 e following ground measurement variables were available for 

statistical analysis (both before and after clipping): 

1. Ch 4, 5, 6, 7. brightness (radiance) 

2. Ch 4, 5, 6, 7. reflectance (bi-directional) 

(or R4, R5, R6, R7) 

3. Ratios : R7/R6, R7/R5, R7/R4, R6/R5, R6/R4 and R5/R4. (or R76. 
R75, R74, R65, R64, R54) 

4. Ratios : R74, R64, R54 (used to diminish the effect of sunlit- 
and shaded-sides of hills) . . 

The following biomass variables were available for each of the 
same 42 stations : 

1. Wet weight grass - . 
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2* Dried weight grass (Biomass, weight/m ) 

3- "Deadness" (Biomass ratio = Dried/wet weight 
4* Species (10) determination 

The following soil variables were available: 

1, Soil type, color, etc. 

2- Soil moisture (although not taken until November 1974, at 
the end of a dry summer, the relative moisture-holding capabilities 
should be consistent) - 

3. Azimuth and slope of hill at the station* 

From careful location of the ERTS CCT outputs we were able to give an 
approximate ERTS-radiance figure to each site, for 7 overpasses of 
the satellite. However, in common with many other experiments, we did 
not have tapes for the year of 1974, but had to make do with those of 
1973, unfortunately a more wet year. Obviously this is not a suitable 
situation but it was all we had with which to work. 

The analyses proceeded using four (4) main data sets. (Prelat 
xi/orked only with Set 1, the present author used all 4). 

2. 3.4. 2.1.1 Data Set 1: Eleven variables - 42 stations 

Biomass data, plus ground reflectance, before and after clipping. 

See Table 2. 3.4.1. 

Data Set 2: Nine variables - 42 stations 

Same data set, now normalized to Channel 4, before and after clipping. 
See Table 2. 3. 4. 2. 

2. 3. 4. 2. 1.2 Data Set 3: Thirty- three variables - 42 stations 

Grass color, biomass data, ground radiance reflectance and ratios, 
soil parameters and species data. 

See Table 2. 3. 4. 3. 

Data Set 4; Seventy-six variables - 42 stations 

ERTS-radiance, ERTS-reflectance (calculated) and ground (truck) 
reflectance, from 8 ERTS overpasses and 3 truck traverses. 

Cross correlation matricies were prepared for each data set, and those 
showing meaningful correlations selected for scatter diagram plotting. 
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Typical matricies x^ere: 
Groxiud Data 


R4 

R5 

R7 



1 

-.77 

.71 

Wet wt. 


- 

-.62 

• 

Dry wt. 


.70 

.90 

. 

Biomass 

ratio 

Ground Data 


i 



Biomass 

Dry rat. R54 

R64 

R74 



.95 -.79 

.77 

.80 

Wet wt. 


-.61 -.50 

.65 

.67 

Dry wt. 


.50 

-.81 

-.82 

Biomass 

ratio 


. 

- 

R54 




.96 

R64 



Data Set 1 

r 

Biomass ratio n 

to R5 (ground) +0.90 

Biomass (dry) to 
R5 - .62 


Data Set 2 

^ 

Biomass ratio to “ 

R74 -- . 82 

R64 - .81 

Biomass (dry) to 

R74 .67 


Ground Data 


Dry wt. BPS* R 5 r ?5 R74 R65 r64 

.95 . -.74 .82 .76 

.87 .82 -.83 -.82 . .77 

. 61 . 

.35 . 

.45 


Data Set 3 

Wet wt. 

Biomass rat* 

Soil moisture 
Broadleaf 
Morning Glory sp. 


* ” Channel 5^ bandpass radiance 
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CORRELATION MATRIX WITH ERTS SATELLITE DATA 
(NOTICE THE DATES REPRESENTED IN THE MATRIX) 


1973 - Wet 

ERTS 1309 
May 28, 1973 

BP7 


1973 - Drier 

ERTS 1669 
May 23, 1974 


BPS 


R5 


R6 


.71** 

.70 

- 


-.63 



.70 


• 


.88 


Data Set 4 


"Deadness" 

Biomass ratio - ground 
May 15-18, 1974 

Truck R75 - ground (3 dates) 


.80 


R5: ERTS 1309 May 28, 1973 

R6 


BP7: ERTS 1165 Jan. 4, 1973 

(wet year) 


linear 

0.71; see figure 2,3.4.19 for actual plot shape 





TABLE 2. 3. 4.1 Five Soil (+Grass) Types in Discriminant Analysis 
Stepwise Choice Sequence 


Single Pair Trio Quartet Quintet 


A 1 

5 steps 



Station 

Altitudes 

BP4 

ERTS-74 - 
5/23/74 

BP4 

ERTS-73 

12/30/73 

s 

Plant 
Species 6 
(Ripgut) 

— • 

. BP7 
ERTSI 74 
1/4/74 

Success 

52% 

69% . ' - 

.93% 

93% 

^00% C" 

A 2 

10 steps 

Station 

Altitude 

BP4 

ERTS-74v'. 

5/23/74 

BP4 . 

ERTS-73 N 
12/30/73 

Ground 

ratio 

R76 

5/18/74 


Success 

72% 

. 92% - 

100% '' 

100% 

(100%) 

5 steps 

BP4 

ERTS-74 

5/23/74 

Biomass 
ratio (D/W' 
5/15/74 

BP4 

ERTS-73 

12/30/73 

Station 

Slope 

BP7 

ERTS-73 

5/28/73 

Success 

60% 

69% 

90% 

86% 

88% 

B 2 

34 vai'. 

BP4 

ERTS-73 

5/28/73 

Biomass 
(Dry trt . ) 
5/15/74 

Station 

s3ope^ 

S7 

ERTS-73 

5/28/73 

R5 

ERTS-74 

5/23/74 

Success 

45% 

69% 

79% 

. 86% 

88% 

C 

34 var. 
f l.?c tnncc 

BP4 

ERTS-73 

5/28/73 

R5 

ERTS-74 

5/23/74 

R7 

ERTS-73 

5/28/73 

R6 

ERTS-73 

1/4/73 

R4 V, 

ERTS-73 

12/30/73 

Slier css 

45% 


88% 

88% 

98% 


Fig. 2. 3. 4. 2- 3 


Fig. 2.3.4.2./I 


Fig. 2. 


Fig. 2. 3. 4, 2.0 


Fig. 2. 3. 4. 2. 7 


NOTE: Crosshatching indicates selection of ERTS data 


Table 2.3. 4.2 


Groupings used In Discrimin a nt: Analysis 


Groups 


■j Steps 


Variables 

Used 


Deletions 


3-t- (Two tests) 


Altitude + Plants 


Altitude + Plants 
and all band pass 
brightnesses 


Table 2. 3.4. 3 Stepwise Choices: F-Value Results at End 



Single 

Pair 1 

Trio 

A 1 

40 

11 ' 

14 

A 2 

67* 

21 

16 

B 1 

15 

10 i 

9 

B 2 

14 

i 

10 

1 

7 

c 

1 14 

8 i 

1 

I 10 


Quartet Quintet 


= I 

■li ! 


* Highest Discrimlnability 


/37 


2, 3. 4. 2, 2 Data Set 1 (Reflectance Variables) 

The following date were available from the biomass and ground reflectance 
measurements at the 44 stations in the Grassland Survey* (See Table 2.3*4«4) 

TABLE 2. 3- 4. 4 BI-DIRECTIONAL REFLECTANCES 

2 

Variable 1 - total wet weight grass in 0*25 m 

2 

Variable 2 = total dry weight grass in 0.25 m 
Variable 3 = ratio dry weight/wet weight 
Variable 4 = reflectance channel 4 before cutting 

Variable 5 - reflectance channel 5 before cutting 

Variable 6 ^ reflectance channel 6 before cutting 

Variable 7 = reflectance channel 7 before cutting 

Variable 8 = reflectance channel 4 after cutting 

Variable 9 = reflectance channel 5 after cutting 

Variable 10= reflectance channel 6 after cutting 
Variable 11= reflectance channel 7 after cutting 
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TABLE 2, 3,4. 5 Station Data for eleven V&riables 


ID 

1 

2 

3 


variables 

j 

4 ^ 


^10 


942 

117.3 

52.8 

0.450 

0.057 

0.9R9 

n.241 

9.338/ 

' 0.192 

9.1ti4 

9.224 

0,299 

941 

105.0 

62.8 

0,590 

0.064 

0.111 

9,209 

0,281 

C.098 

0.149 

9.224 

0.797 

947 

213.2 

120,0 

0.5G3 

0.065 

n.iii 

0.239 

0.537 

0.128 

9.187 

0.251 

0.321 

947 

70.7 

57.3 

0.710 

0.07S 

0.130 

'^,206 

0.27S 

9.110 

0.173 

0.237 

0,^01 

949 

110.3 

71.0 

9.595 

0.072 

n,i05 

0.213 

0,2f4 

9.105 

0.152 

0,214 

n,270 

951 

150.0 

81, n 

0,5 46 

0.057 

1,n93 

0,230 

0,335 

0.123 

0.175 

9.251 

9.3lf 

953 

131.1 

82. 7 

0.G31 

0.067 

n,ii4 

0,206 

0.2P7 

9.nno 

9.151 

9.213 

n 9 p r 

■954 

83.0 

55.5 

0.062 

0.087 

0.143 

9.. 239 

0.590 

0.147 

0.20? 

n.279 


955 

70.1 

54.6 

0.779 

0,091 

0.149 

9,213 

0.29G 

0,129 

9.184 

0.245 

0,331 

980 

542.4 

13S.4 

0.404 

0.05 4 

9.965 

0.243 

0.414 

9, ''15 

0.13P 

n 0 no 

q 

982 

421,8 

173.0 

0.410 

n.066 

n,n7l 

0.203 

0,474 

0.103 

9.142 

0,1^0 


M4 

282,5 

121.0 

0.428 

0,055 

0,070 

0.253 

0.300 

0,19? 

n.iG 2 

n,22f> 

q opo 

916 

50G.4 

1.28.2 

0,418 

0,057 

O.nqi? 

0.246 

0,33? 

n.nriT 

9,142 

0.214 

0.768 

917 

184.6 

103.2 

0.559 

0.051 

0,997 

9.199 

0.314 


0.14G 

'’.214 

0,783 

920 

103.9 

10 4. r 

0,530 

0.048 

n.ns5 

9.187 

9 , ? R 1 

0.905 

9.144 

9.123 


930 

333.7 

160,0 

0,479 

0,053 

0.<'85 

o,2ir 

9,302 

0 , 

0.124 

1 . 1 4 

n,2T5 
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Fig. 2.3*4.! Biomass ratio (Dry/Wet) versus R5 (r-0.90) 
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Fio*. 2*3.4. 2 Wet weight versus CH 5 (r — ^0.77) 
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Fig. 2. 3.4. 3 Wet weight versus CH 7 (r-0.71) 
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Fig* 2. 3*4, 4 "Deadness*' (Biomass ratio, or Dry /Wet) versus R 4 (r=0*70) 
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Fig, 2*3<4.5 Dead weight (Biomass) 


versus R 5 (r=-'0.62) 
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2. 3.4* 2* 3 Data Set 2 (9 Variables) 

The data set after normalization to Channel 4 (for reduction of lighting 
effects and preliminary atmospheric correction. (See Table 2. 3. 4, 6) 

TABLE 2.3, 4. 6 BI-DIRECTIONAL REFLECTANCE DATA NORMALIZED TO CHANNEL 4 

2 

Variable 1 = total wet weight grass in 0.25 m 

2 

Variable 2 »" total dry weight grass in 0.25 m 
Variable 3 = ratio dry weight/wet weight 
Variable 4 = reflectance channel 5/4 before cutting 

Variable 5 = reflectance channel 6/4 before cutting 

Variable 6 - reflectance channel 7/4 before cutting 

Variable 7 = reflectance channel 5/4 after cutting 

Variable 8 = reflectance cbannel 6/4 after cutting 

Variable 9 - reflectance cbannel 7/4 after cutting 


134 



TABLE 2.3.4. 6 cont. Preceding Set Normalized to Channel 4 
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3.024 

3.223 

4.681 

1,543 

2!311 

3*290 

3,650 

5.1G1 

1-567 

2.146 

2.606 

4.135 

5.075 

1.456 

2,112 

2.75F 

3.628 

5.100 

1,549 

2.104 

2.P43 

4,54 3 

6.467 

1.517 

2,0^4 

2.723 

2.308 

3,205 

1,‘^15 

2,0P3 

2*758 

2.462 

3.542 

1.531 

2. Jin 

2.QF8 
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lUABLB 


5 

C, 

6*160 

+ 

6. cog 

. 

6.000 


5.920 

* 

'v^840 

. 

5-760 

#• 

5*680 

* 

5.600 

* 

5.520 

- 

5-440 


5, 360 


5*280 

* 

5.200 

. 

5* 120 

. 

5-040 

« 

4-960 


4* SDO 

* 

4-800 

* 

4- 720 

. 

4*640 

4 

4.560 

+ 

4.400 

•* 

4.400 

* 

4.320 

— 

4 , 240 


4. 16Q 


4. 080 

.* 

4» 000 


3 * 920 

• 

3* 840 

- 

3. 760 


3.680 

- 

3. £00 

. 

3*520 

’ * 

3-440 

- 

3*360 

y 

3*280 

m 

3.* 200 


3.120 

, 

3.040 

• 

2 » 960 

+ 

2. 880 

. 

2. 800 

. 

2*720 

* 

2*640 

. 

2.560 


2.480 

. 

2-400 

. 

2.320 

« 

2.240 

• 

Z.160 



VARIABLE 3 


U-360 


O.AIO 


0.510 


0 . ^60 


0*560 


0.610 


0.710 


0*660 




0-760 


O.dlQ 




1 1 


!• 


1 1 1 
1 1 

L 1 


11 


1 

I 1 


Z 1 

1 


1 

i I 


I I 

I 


,,, ,, 1-^. ^ * .4** e . * *9'. **.4'.. 

0*310 O.AIO OoSlO 0*610 

0-360 Oa'iBO 0.560 


0*710 




0.660 


■ *♦**.» 4 " 

0.810 


2.3.h.6 ^’Deadness" -versus R 61)- (r=".8l) 




LiL.. 


6. 160 
6-080 
6-000 
5.92 0 
5.640 
.5*760 
5.680 
5.600 
5.520 

5.440 

5.360 
5.260 
5.20 0 
5.120 
5*040 
4*960 
4.660 
4.G00 
4* T20 
4*64 0 
4*560 
4. 400 

4.400 

4.320 

4.240 
4*160 
4-OBO 
4.000 
3.920 
3.840 
3.T60 
3.650 
3.600 
3-520 

3.440 

3.360 
3.280 
3.200 
3*120 
3*040 
2.960 
2.860 
2.800 
2.720 
2.640 
2.560 
2.460 

2.400 

2.320 

2.240 
2*160 


0.760 


.rtIA8L£ 

16 


VARIAQL6 ^ 





18.000 

50.000 

98.000 ■ 138.000 

178.000 


218-000 




36.000 70.000 

110*000 

158-000 

198-000 
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10.2G0 







10-200 

10.050 






* 

10.050 

9.900 







9-900 

9.750 



1 




9-750 

9.600 






- 

9.600 

9*^50 

+ 





^ ■ 

9-450 

9-300 






. m 

9-300 

9-150 






* 

9.150 

9«OOD 






m 

9- 000 

b.650 






m 

a-050 

0. 700 

4, 





> 

8.700 

S-550 






. 

8-550 

S.40C 



1 



« 

8.400 

8.250 




I 


** 

8-250 

B. 100 






- 

0.100 

7. 960 

1- 



1 


. *■ 

7*950 

7.800 

* 






7-BOO 

7-650 



1 


1 

» 

7-650 

7-500 

* 






7-500 

7.350 






- m 

7-350 

7-200 




1 1 


*■ 

7.200 

7,050 

m 


u 



. ^ 

7-050 

6.900 






* 

6.90Q 

6-750 

• 





. 

6-750 

6*600 

m 


1 i 



— 

6-600 

6.450 

h 


1 

1 



6.450 

6.300 






« 

6.300 

6-150 

m 


11 



- 

6-150 

6-QOO 

m 


11 

1 


. 

6.000 

5.050 


1 1 

1 

i 


- 

S.B50 

5-700 

+ 


1 . 1 

I 



5.700 

5- 550 






. 

5*550 

5-400 

• 



1 


— 

5-400 

5-250 

* 

I 

1 



■ 

5-250 

5-100 

■ 

L I 




— 

5-100 

4*950 






> 

4-950 

4.800 



1 



• 

4.800 

4.650 

m 

1 




— 

4-650 

4.500 




1 


— 

4.500 

4-350 

* 

1 1 




• 

4-350 

4.200 

h 





*• 

4.200 

4.050 






— 

4-050 

3. 900 

* 





. 

3-900 

3-750 

* 

1 




— 

3.750 

3.600 


1 1 




■ 

3^600 

3-450 


1 





3-450 

3.300 


1 1 




* 

3-300 

3- 150 

m ' 




- • 

. 

3.150 

3*000 

m ' 





« 

3-000 

2-850 

* 





m 

2-350 

2-700 






<• 

2 -TOO 









16-000 

58.000 

90.000 136-000 

17S-DOO 


218-000 




36.000 78.000 

118-000 

158-000 

198-000 




2. 3-^. 7 Dead -ffeight (Biomass) versus H jh (r=.68) 




2. 3, 4.2. 4 

Data Set 

3; Original set; plus biomass, plus species and soil 
types (see Table 2.3.4. 7) 

TABLE 2.3. 

4.7 BI-DIREGTIONAL 

REFLECTANCE, CHANNEL BRIGHTNESS, BIOMASS DATA, ETC. 

VARIABLE 

LINE NO. ■ 

• DESCRIPTION 

1. 

2 

1 

Pseudo-C.I-E. coordinates **EASNX” 

2. 

2 

2 

Pseudo“C-I.E. coordinates "RASNY^^ 

3. 

2 

3 

2 

Net weight of grass/ 0.25 m 


2 

4 

Dry weight of grass/ 0.25 m 

_5. 

2 

5 

Biomass ration "deadness^^ = D/W (3/4) 
... 2 

6. 

3 

2 

CH 4 radiance (w/ cm /bandpass ERTS filter) 

7. 

3 

3 

CH 5 ” " 

8. 

3 

4 

GH 6 " " 

9. 

3 

5 

CH 7 " " 

10. 

3 

6 

Reflectance R 4 

X J* « 

3 

7 

R 5 

12. 

3 

8 

R 6 

13. 

3 

9 

R 7 

14. 

4 

1 

Ratio R7/R6 

15. 

4 

2 

R7/R5 

16. 

4 

3 

R7/R4 

17. 

4 

4 

R6/R5 

18. 

4 

5 

R6/R4 

19. 

4 

6 

R5/R4 

20. 


6 

Moisture percentage, Nov 1973 

21. 

5 

7 

Dip of slope, degrees 

22. 

5 

8 

Azimuth of slope , degrees from north 

23. 

5 

9 

Altitude, feet /sea level 

24T~~ 

6 

2 

Species count, percent, of soft chess(Bromus mollis) 

25. 

6 

3 

Ryegrass (Lo lium multiflorum) 

26. 

6 

4 

Wall barley (Hordeura hystrix) 

27. 

6 

5 

Foxtail barleyCH. leporinura) 

28. 

6 

6 

Slender wild oat (Avena barbata) 

29. 

6 

7 

Ripgut grass (Bromus rigidis) 

30. 

6 

8 

Morning Glory (Convolvulus; sp.) 

31. 

6 

9 

Needle (Erodium sp.) • 

32. 

6 

10 

Sheep sorrel (Rumex sp.) 

33. 

6 

11 

All other broad leaved plants. 
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TABLE 2 . 3 . 4 . 7 co n t . 


<■.9^25 92700 l.C 0.19C 0,290 l.O O.B60 1.Z85 0-515 

2.. ; 0.180 0.226 U7.3 52-8 C.650; 942»:RASNY 

5.5 542 0.412 0.675 1.647 3.034 0.059 0,078 0.254 0,337 
y.. ; 1,227 4.344 5,716 3.272 4.307 1.316 

C. C-942;CSW-LTB/TPB 10YR 2/2 9,9 3.2/201SSW 468 

942 — 10 ec — 10 


941; 94600 

; 0-200 


0.24C 


0.820 


1.167 0,515 


15. GIASS UP HILL TO NE HVi 


1516569 333951 


15. CRASS NR MW TRAILER HIL 


; 0-200 0,317 105. 0 62, B C.590; 54l,:PASNY 

; 541 0,514 1 .057 1 ,566 2.762J0.065 0.108 0.220 0.280 


; 1,274 2.599 ^.232) 2.041 3.322 1-628 
C-941;CLR-LT8/tPB fiOYR 2/21 

941 80 10 — 10 


4.0/200SSH 476 1516596 334078 


947; 100000 1,0 0.25C 0,480 1,0 0.985 1.505 0,515 

; 0.180 0. 231 213.2 12 C.O C. 563; 947, :PASNY 

; 947 0 .534 1,078 1.900 3.541 0.067 0.108 0-252 0.336 
; 1-335 3. 109 5,022 2.32 8 3.761 1.615 

C-947;CL9.-LT3/TPB SYR 3/2 11.2 4-1/095 E 400 

947 10 20 — 60 — 10 

9477 101000 I-O 0,350 0-665 1.0 0.937 1.337 0.515 

; 0.224 0.345 79.7 57,3 0.713; 947,;RASNY 

; 947 0,736 1-471 1,803 3.154 0.080 0.127 0.217 0-277 
r 1.277 2-180 3.443 1.7C7 2-697 1.579 

C-947;CLR-LTB/TPB 5YR 3/2 11.2 4.1/095 E 480 

947 10 50 20 — 20 

949; 102200 1.0 0-30C 0.500 1.0 0.940 1.220 0.515 

0.227 0.309 119,3 71,0 C.595; 949,:PASNY 

; 949 0,635 1.120 1.009 2.8B5 0.C74 0.102 0.229 0-263 
; 1.150 2-574 3.578 2.239 3.112 1,390 

C-949;CLR-LTfi/TPB lOYR 2/2 17.2 5. 2/C94 E 493 

949 70 20 -- 10 


IS. GREEN GRASS SITE 4 


1516321 333739 


15. VERY DRY AOJ PATCH 40RY 


951; 103400 1.0 0.26C 0.475 l-O 1.060 1-650 0.515 

; 0.178 0.264 150. C 81-9 0,546; 951,:RASNY 

; 951 0.554 1-067 2,052 3,875 0.059 0.090 0.243 0.334 
; 1.37/ 3.703 5.672 2-69C 4.120 1.532 

0-S51;CLR-LTe/TPB lOYR 3/2 10.9 3,6/076 E 512 

95 1 5 5 9C — — 

953; 104900 l.C 0.310 0.600 1.0 0.960 1-428 0.515 

; 0,205 0.322 131.1 82-7 C.631; 953, sRASNY 

; 953 0.655 1.333 1.850 3.364 0.069 0.112 0.217 0.286 
; 1.319 2-562 4.151 1,943 3.147 1.620 

0-953;CLR-LTe/TPfl tt0VR2/2) (9-81 4-B/074 E 523 

953 70 20 — 10 — 


1516321 333739 


15* site 6 NR OAK TREESNW 


1516128 333642 


15* SHORT GRASS NR OCT FNCt 


1515877 333482 


IS* GRASS SITE 8 SHRT BLKTI 


1515690 333295 


954; IIOOOO 1.0 0-4U 0.747 1.0 1.090 1-510 0.515 IS* GS SITE 9 CRX RD GRY F' 

; 0.23G 0.341 83.9 55.5 C.662; 954,:RASMY 
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; S54 0.857 l*6^i5 2.115 2,553 C.0?0 0.139 0,242 0-299 
; 1.234 2-154 3.342 1.745 2.707 1.551 

0-554 ;CLR-UTF/TP0 (lOYP 2/2> (9.8} 2.3/135SF 516 1515557 333097 

954 5 — 95 . — . 


955; ILIOOO l.O 0-45C 0.630 1.0 1.070 1,537 0. 515 15. SITE 10 HLF TO TOP DEAD 

; 0.236 0.354 7C^l 54.6 0.7793 955,:RASMY 

; 555 0.938 1,021 2.073c: Tr6l H^ 0.C54 C.147 0.207 0.293 
3 1.413 L.95C 3.119 1.400 2-207 1.568 

D-955 ;CU-LTB/TPB 10 YB 2/2 13.5 3-6/080 E 519 1515464 332941 

555 40 5 40 — 15 


98 0 5 U2500 1.0 0.26C 0.340 1,0 1. 195 2.105 0.515 15- SITE U ACROSS Y FORK 

; 0.164 0. 179 342.4 138.4 C.404; 980. sPASNY 

; 980 0,554 0.701 2.326 4.924 0.056 0.064 0.235 0,409 
; 1.733 6.421 7.327 3 .704 4.227 1,141 

G-980; SLR-DAE/TUS 5Y 5/1 0,3 3, 2/208SW W03 1514382 332563 

980 95 5 


9 8 25 U4000 l.O 0.320 0.375 1.0 1.450 2.510 0.515 IS. SITE 12 TO W OF 11 GREEf 

; 0.170 0.166 421,8 173.0 0, 410; 982t:RASKY 

; 982 0.675 0.855 2.043 5.857 0.068 0.070 0,285 0.468 
; 1.644 6. 70S 6.046 4.079 4,165 1,021 

G-98 2;S LR-OAE/TUS 5Y 4/1 15. 1 11.3/307NWW 498 1514296 332406 

982 20 15 60 — 05 


914; 93000 l.C 0,190 0.270 1.0 C.980 1.540 0.S16 15. GRASS PH«2 

; 0.163 0.193 282.5 121.0 0.428; 914,:RASNY 

; 914 0.412 0.632 1.890 2.622 0.C52 C. C68 0.267 0.359 
; 1.346 5.280 6.873 3 .930 5.100 1.300 

J-914;CLR-AV0/QSC 10 YR 5/4 6.9 4.5/194 S 368 1512779 329482 

914 30 20 40 5 05 


916; 95000 l.O 0.21C 0,275 1.0 1.000 1.495 0.516 15. GRASS PH #4 

; 0.17 9 0. 196 30 6 .4 128.2 0.4 1 6 ; 9l6tiPASNY 

; 916 0.453 0.643 1.931 3.518 0.059 0.066 0.259 0.337 
; 1.301 5,091 5.752 3 .912 4.420 1.130 

J-9l6;CLR-AVO/0SC 10VR4/2 6.9 9.5/141SE 3701513081329558 

916 30 10 30 30 


917; 101000 1,0 3,190 0.415 l.O 0.820 1.370 0.516 15. GRASS PH 46 

; 0.159 C. 2 80 184.6 103. 2 0. 559; 917t:RASNY 

; 917 0.412 0.940 1.566 3.230 0.052 0.095 0-210 0.313 
; 1.494 3.30C 6.008 2.209 4.023 1.821 

J-917{CLR-AV0/CSC lOYR 4/2 6.9 9.5/14ISF 370 1513081 329558 

917 60 20 — 10 5 05 


920; 101350 1,0 0-210 0.420 1,0 0.880 1.400 0.516 15- GRASS PH<f 8 

; 0.170 0-275 193-9 104-6 0, 539; 920f:PASHY 

; 9ZC 0.453 0.951 1.687 3,299 0.049 0.083 0.197 0,280 
; 1.^21 3.382 5.676 2.379 3.993 1.678 

J-920;CLR-AVC/OSC lOYR 5/4 7.9 10. /130SE 340 1513506 329426 

92C 10 TO 10 10 


930: 102800 1-0 0.22C 0-380 1.0 1.395 1.850 0.516 151 GRASS P H»1 0 A T T FElT L 

; 0.153 0.216 333.7 160.0 0,479; 930,;RASNV 

; 930 0.473 Q.866 2.123 4,336 0.052 0.082 0-249 0.382 
; 1.533 4.642 7-303 3.028 4.764 1.573 

K-930;CSH-AV0/(JSC lOYR 5/2 9,2 S.7/123SE 368 1512606 329331 

■iBO 5 30 60 — 5 
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921; lO'iaCO 1.0 0.25C 0.470 1.0 1.000 1.670 0-516 15. GRASS PH# 12 TO LIKE 

; 0.171 0.262 271-7 155-0 C-570; S3l,:RASNV 

; 531 0.534 1.057 1.931 3.921 0.060 0-095 0.233 0.338 
; 1.554 3.564 5.53C 2-452 3.872 1.579 

K-931; CLR-AVD/QSC lOYP 5/4 7,8 3.8/190 S 397 1512327 329278 

931 10 30 40 10 10 


922; 1057CO l.C 0.225 0-370 1. C C.960 1.637 0.516 15- GRASS PHfl 14 

; 0.169 0,227 406.5 100.0 0.443; 932,:RASNY 

; 9320.483 0,845 1.850 2. 845 C.C52 0.074 0.214 0.326 
; 1-519 4.419 6.243 2.909 4.109 1-413 

K-532;CLR-AVD/QSC lOYR 4/2 8.6 5. /044NE 393 1512066 329231 

932 25 25 — 40 5 5 


934; 111400 1-0 0.26C 0-370 X- C 0.980 1-605 0.516 15. GRASS PH# 16 

: 0.191 0.224 307.3 130.1 0,423; 934,;PASNY 

; 934 0.554 0-845 1-890 3,772 0.060 0-074 0.223 0.334 
; 1.502 4-509 5,550 3.002 3-696 1.231 

K-934;CLR-AVD/QSC lOYF 4/2 8,6 1. /285NWH 400 1511816 329245 

9 34 40 40 — 20 __ 


940; 113100 1-C O.iac 0.205 l.O I-OIT 1.680 0.516 15'. GRASS AT LAKE PH#18-TRK5 

; 0.154 0.150 344.8 122.3 0.355; 940. : RASNY 

; S40 0,392 0-494 1.965 3-944 0.C40 0.041 0.221 0.324 
; 1-467 7-988 8.118 5-443 5.532 1,016 

M-940;CLR-AVD/Q5C 5Y 5/1 (7.1) 4.1/312NW 404 1510937 328267 

940 — 40 30 20 10 


942; 114000 1.0 0.24C 0.284 l.O 1.320 2-16C 0-516 15. GRASS PH# 1 NEW ROLE 

; 0-156 0-155 43C.8 149.1 0.346; 942,:PASNY 

; 942 0.514 0.662 2.580 5.050 0.C53 0.054 C.262 C.404 
; 1-544 7-456 7-634 4-828 4-944 1-024 

H-942;CL».-AVD/QSC 5Y 5/1 (7-11 4.1/312NH 404 1510937 328267 

942 30 30 20 20 


943; 120400 1.0 0.340 0.513 1.0 1.460 2.46C 0.516 15- GS VRY HNDY PH 4 

; 0.171 0-212 313.9 138-5 0.441; 943,;RASNY 

; 943 0.716 1-14E 2.864 5.742 0.071 0.091 0,2 60 C.453 
; 1-618 4.951 6,405 3 .059 3 .958 1.294 

H-943;CLR-AVD/QSC lOYB 4/2 (8-2) 2.5/316NW 423 1511102 328403 

943 30 30 30 — 10 


944; 121700 L.O 0.305 0.467 l.O 1-270 1.960 0.516 15- GRASS PH# 6 HINOV 

; 0.183 0-230 290. C 132.3 C.456; 944,;RASAY 

; 944 0.645 1,050 2.478 4.590 0.C61 0,081 0.235 0.352 
; 1.494 4.331 5.720 2.90C 3-029 1-321 

H-944; CLR-AVD/QSC 10YB4/2 (8.2) 2.5/316NW 423 1511102 320403 

944 20 40 — 40 


9C5; 93500 l.O 0.28C 0-390 1.0 0.900 1.340 0.521 15. PH# 7 NO WIND DAY 3 

; 0,220 0,253 354.5 168-2 C-474; 905,:RASNY 

; 905 0.594 0.087 1-728 3.161 0.085 0.102 0.271 0.368 
; 1.360 3-613 4.337 2,657 3-190 1-200 

I-905;CLR-DAE/TH 7YH 4/1 10.4 0.5/190 S 473 1513772 330962 

9C5 50 40 — 5 5 — 


9C6; 94800 l-O 0.23C 0-480 1,0 0.080 1-340 0,521 15, SPARSE GS PH# 9 + SOIL 

; 0.180 C.303 16E-C 50-2 0,476; 906, :RASNY 

; 906 0.493 1.070 1.6:37 3.161 0.064 0.111 0.238 0.328 


TABLE 2. 3. 4. 7 cont 


t 1.277 2.966 5.120 Z,15A 3.72^ 1.729 

1-906; CLR-Oar/TM 5Y 5/L 10.4 6.5/H055SW 466 1513719 330703 

9C6 50 40 — 10 


9C0i 100200 l.O 0.17C 0.310 1,0 0.660 1,080 0,521 15. PH « 12 HILL SLDPNG H 

; 0.179 0.266 209.5 95.0 C.457; 9ca,;RASNY 

; 908 0.372 0.717 1.241 2.562 0.C4T 0.072 0. 170 0.257 

; 1.510 3.553 5.467 2.352 3.619 1.539 

r-90B:CLP-OAE/TM lOYP 5/2 0.7 0-5/204S9W 435 1513610 330562 

900 20 30 — SO 


909! 101600 1.0 0.25C 0.450 1.0 C.950 1.315 0.521 15- PH# 14 LAST BFR FRWV 

; 0.197 0.200 192,5 113-6 0.590; 909t:P.ASNY 

; 909 0.534 1.014 1,029 3.103 0.070 0.105 0.265 0.326 
; 1.231 3.099 4.641 2,517 3.770 1.490 

1-909 !CLP.-AVD/QSC lOYR 6/2 8.3 13. /Z07SSW 405 1513456 33Q400 

909 45 45 — — JO 


99 1; 103300 l.O 0.24C 0.360 1. 0 0.990 1-590 0.521 15, ACRX DIRT RD PH 16 

; 0.182 0.224 22C.4 98.9 0.449; 991t:PASNY 

; 99 1 0.5L4 C- 823 L.9LQ 3,737 0. 057 0.074 0.209 0.336 
; 1.609 4.563 5.062 2.836 3.642 1.284 

H-991;SLR-DAE/TUS lOYP 4/2 7.9 4.7/104SFE 452 1514096 331612 

991 30 20 35 5 5 5 


992; 104800 1,0 0.14C 0,160 L.O 0.840 1.330 0.52L 15. VRY GRN VALLEY PH# IB 

; 0.153 0-151 463.5 175.7 0.379; 992t:PASNY 

: 992 0.311 0.399 1.606 3.130 C.C40 0.045 0.197 0.317 
; 1,608 7,113 7.905 4,424 4.917 l-llL 

H-992;CSW-UAE/TU5 lOYR 3/2 9,3 8. /056 E 445 151405B 331451 

992 25 50 — 25 


994 ; 110400 l.O 0.3CC 0.410 1,0 1,270 1,910 0-521 15- PH 20 CL3 SKY UP HILL 

; 0,189 0-213 404.3 165,4 0-407; 994t!RASNY , 

; 994 0.635 0,929 2,478 4-474 0.070 0.083 0.266 0.426 
; 1.601 5. 145 6.C82 3.214 3 .000 1.182 

F-994;SLR-DAE/TUS 5Y 5/1 10.1 4.1/061MFE 4671513964331246 

994 30 60 10 

996; 111700 1.0 0.22C 0-420 1.0 1.230 1-865 0-521 15- PH I CLR SKY SHRT OF DR 

; 0,149 0.231 423.2 197.2 C.466; 996,;RASNY 

; 996 0.473 0.951 2.39T 4.3T1 0.049 0.079 0-244 C.361 
; 1.478 4.574 7.357 3.095 4.978 1.608 

H-9g6;CLR-DAE/TH 5Y 5/1 11.2 T. 1/074NFE 477 1513839 331103 

996 10 30 — 50 


9E4; 115500 1,0 O.ZC: 0.330 1.0 0.000 1.410 0-521 15i PH 4 

; 0.172 C-233 317.7 126.4 0.398; 984,;RASNY 

T 984 0.433 0.760 1.687 3.322 0,044 0 ,062 0.1 73 0 .2 95 
; 1.709 4.68 1 6.716 2,73 9 3 .930 1.435 . 

G-984;SLD-DAE/TUS lOYR 4/2 5.8 9.5/23TSWW 480 1514242 332134 

9 64 30 60 — — 0 5 5 


950; 93200 1-0 0.160 0.200 1.0 0.870 1-300 0.522 15. GREEN AT LAKE TRK 6 PH6 

; 0. 167 0.177 325.5 127. I C. 390! 9S0,:RASAY 

; 950 0,352 0,484 1,667 3 .069 0.056 0.061 0-284 0.381 

; I.34Q 6.276 6.861 4.685 5,119 1.093 

M-950; CLR-AVU/QSC 5Y 5/ 1 7.1 4. 1/312NW 404 1510937 328267 

950 5 5 5 EO 5 




TABLE 2. 3.4. 7 cont. 


95J; g^BCO l.f 0.14C o.iao 1.0 0.930 1.620 0-522 

; 0-132 0.144 403.5 145.3 C.360S S52,;PA5NY 

; 952 0.311 C.44L 1-789 '.006 0,042 0.048 0.262 0.397 
; 1-517 8.354 9.392 5.505 6. 189 1.124 

M-952;CLR-AVC/CtSC 5Y 5/1 7.1 4.1/312NW 404 

952 10 20 20 30 20 

9545 100100 1.0 0.25C 0.330 1.0 1.040 1.53C 0.522 

5 0.195 0.214 367.1 16.4.2 0,447; 954, : RASi^Y 

5 954 0.534 C.76C 2.012 3.599 G.070 0,079 0.27fl 0.364 
5 1,310 4,603 5.U5 3.513 3.949 1.124 

K-9545CLR-AVD/QSC lOYR 4/2 6.2 2.5/316NW 423 

954 20 10 — 40 10 10 10 

9705 102600 1.0 0 .24C 0.450 L.O 0.8 20 I,22C 0.522 

5 0.198 0.302 86,5 45.1 0.521; 970,JPASNY 

; 970 0.514 1.014 1.566 2.885 0.061 0 .096 0.205 C.277 
; 1,349 2.883 4.539 2 .137 3 -3 64 1.574 

F--g70iSLR-DAE/TUS lOYP 4/2 7-2 7. 1/028NNE 475 

970 50 30 2G 

972; 103600. 1.0 0,.Z7C 0.500 1.0 1,020 1.515 0,522 

; 0,189 0.2B4 152. Q 74,1 0.488; 972, : RASNY 

; 972 0.574 1*120 1.971 3.564 0.C£3 0.097 0.237 0.316 
; 1.335 3^.263 5.004 2-445 3 .750 1.534 

r-9T2;SLR“DAE/TUS lOYP 3/2 10.6 14.0/030NNE 496 

97 2 40 40 — 10 10 


973; 105200 1.0 0.26C 0.450 1-0 1. 1 10 U6B0 0,522 

; 0.178 0-251 354-9 171.4 C. 4 83; 97?, ; RASNY 

; 973 0.554 1.014 2.154 2.944 C.C61 0.C87 0.238 C.34S 
; 1.464 3.992 5.731 2.727 3.914 1.436 

F-973;SLR-DAE/TUS 10YR 4/2 7.2 9. 5/111 SHE 4B6 

973 5 20 — 70 5 

929; 1113C0 1.0 0.275 0.470 1.0 1.045 1,540 0-522 

; 0.193 0.269 191-7 81.8 0.427; 929,SRASNY 

; 929 0.584 1.057 2.022 2.622 0.062 0.088 0.214 0.306 
; 1.426 3,459 4.094 2.426 3.432 1.415 

B-929;SLR-GME/TBU / (9.8J 14.C/175SSB 380 

929 20 20 10 10 30 10 __ 

93i; U2500 1.0 0.25C 0.420 1.0 1-170 1.750 0,522 

; 0.171 0.235 294.6 131.4 C.446; 93li:pASNY 

; 931 0.534 0i951 2.275 4.106 0.056 0.079 0.239 C.345 
; 1,444 4.372 6.2C6 3.027 4.298 1.420 

B-931;CSW-GME/TBU 5YP 4/1 11,8 11.3/338NNW 374 

931 — 50 50 — — * 

933; 113800 1.0 0.365 0.685 1.0 0,910 1.320 0.522 

; 0.239 C.346 86.2 46.2 C. 559; <33, ; RA SN Y 

; 933 0.807 1.5L3 1,748 3.115 0,085 0.128 0,187 0,270 
; 1.44B 2.117 3. 161 1.462 2.184 1.493 

B-933;CLR-LTE/TPB SYR 4/1 11.9 5.7/C12NNE 444 

933 20 20 — 50 , 10 

936; 115200 1.0 0-360 0.640 1.0 0,890 1.315 0.522 

; 0.233 0^337 6C.3 42.2 C. 700; 926,tRASNY 

; 936 0.756 1,417 1.700 3,103 0.078 0.116 O.UZ 0.266 
; 1.459 2.288 3.398 1.56E 2.329 1.485 
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2.3. 4, 8 Wet weight versus dead weight (r=0.93) 
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243.4*9 '^Deadness” versus CH 5 (BP 5)(r=0*87) 
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2-3-4-10 ^Deadness^* versus R J5 ( 2 ^™ -63) 
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2.3.4.11 ''Deadness*' versus R 5 (r=.82) 
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2- 3 *^.12 *^Deadness*' versus E 7^ (r™, 82 ) 
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2.3.4.13 Wet veight versus R 5 (r=-.74) 
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2-3.ilflU Soil moisture versus E 5 (t^.6i) 
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2*3*4.15 Morning Glory (sp.) versus R 65 (r=,k3) 
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2. 3, 4, 2. 5 Data Set 4 (76 Variables); Original set of ground data plus ERTS brightness 

and reflectance. 

TABLE 2. 3. 4,8 ERTS BRIGHTNESS, ERTS RERLEGTANGES AMI GROUND BERLECTANCES 

( TOTAL QE 76 VARIABLES) . 

LETS BRIGtlTNESS ERTS REFLECTNCE. GROUND REFLEGTNCE 
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2,3. 4.3 Discriminant Analysis using Vegetation on Main Soil Groups 
(BMD07M) . 

Stepwise linear discriminant f^nalysis to calculate canonical transforms 
of the data (65 variables) was carried out » with results summarized as 
Tables 2v3-4.1, 2,3*4. 2, and 2. 3*4. 3 above.^ . The programs used were 
adaptations of the UCLA Biomedical set> genex-ally available on all TJ.S. 
computer systems (BMD series) . 

One pictoral aspect of the output is a two-dimensional plot of 
the first and second canonical variables^ which enable the viewer to see 
the best-fit-decision plane through the N-dimensional space of the 
calculation* Such output plots follow here as figures 2*3.4.23 through 
2.3.4.26* 


2. 3 *4.4 Computational Procedure Discriminant Analysis (BMD07M) 


Figure 2.3*4*23 
Figure 2.3*4*24 
Figure 2*3*4*25 
Figure 2,3*4-26 


5 groups, 65 variables, 5 steps (Al) no deletions: 

100% success in separating the 5 (soil and grass) groups* 

3 groups 2 test groups (also same as Al), no deletions: 
100% success 

5 groups, 54 variables (Bl) deleted altitude (1) and 
plant species (10): 88% success, ERTS data included* 

5 groups, 51 variables (C) , deleted bandpasses (4), 
altitude (1), and plant species (10); 98% success, 

ERTS data included - 


The GOWDEL (control and delete) option was used for Runs Bl, B2 (not sho^-m) , 
and G. Table 2. 3*4*1 should be read to show the results together with the 
above figures. 
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2.3 -h, 23 Canonical Plot (BMDOjM) j 5 group separation 63 varialiles 

Run. Al, no deletions. (Success 100^) 






i.k.2k CK.oaioal Plot (BMD07H) ; 3 groups + 2 test groups; 65 varisilos, 
10 step. Run A2, no deletions (Success - 100^) 
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2.3. ^.2^ Canonical Plot (BMDOTli) ; 5 groups « 51 variables, 5 step. Run C 
deietefl, bandpasses ( It-) , altitude (l) and plant species (lO) 
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2.3/5 CORRELATION OE ERTS SPECTRA WITH ROGK/SOIL TYPES IN 
CALIFORNIA GRASSLAND AREAS 

2, 3,5.1 ABSTRACT 

A seasonal study of ERTS-CGT data^ accomplished by means of four 
band spectra plots of normalized reflectance, indicates that in the San 
Francisco Bay and adjacent Coast Range grassland areas, soils mapping 
or classification by computer techniques is possible at the end of the 
dry or grass dieback season* Excellent correlation is sbom betx^een 
ground reflectance measurements and GCT data at three test sites and two 
different soil types: serpentina and sedimentary. The uniqueness of 

their spectra is then demonstrated by the successful application of 
STANSORT, a computerized classification technique developed by the 
Stanford Remote Sensing Laboratory. 


2 - 3 . 5. 2 


MTRODUCTION 


The primary purpose of this investigation ws to deteinniTie if the 
serpentine exposures and soils on the San Erancisco Peninsula could he 
detected uniquely by means of EETS imagery and/or the related CCT. data. 

In doing so it was also hoped to evolve a methodology which would be useful 
in conducting similar studies in the future. As envisioned, the imagery, 
(individual bands and color composites) were to be studied first to determine 
if these serpentine areas could be detected visually and then a study made to 
determine if any uniqueness existed in their four band spectra- This property 
if existent, could then be utilized as a basis for the development of a 
computerized classification program to automize the detection and mapping 
procedure* 

In the course of this investigation, it became evident that the seasonal 
response of the vegetative cover could be most important in obliterating 
or enhancing the information relating to the serpentine soilm Therefore, 
a careful systematic, spectral study of the CCT data for the yearly cycle was 
undertaken by means of four band radiance and reflectance plots of the test 
areas- Off season correlations of serpentine soil spectra vs- serpentine 
soil/ grass spectra were also made possible by means of a fortuitous grass 
fire in the study area which had exposed a large area of bare soil- These 
correla;tions ultimately led to the conclusion that the soil/ grass spectra 
were in fact essentially soil spectra at the end of the dry or dieback season. 

After an extensive ground measurement program had substantiated the 
unique character of the serpentine soil spectra the study was expanded to 
include a sedimentary area on the east side of the Coast Range upon which 



a yearly controlled burn occurred. The same seasonal trends were evident 
and a strong correlation between EIH!S reflectance spectra and ground measured 
..spectra, was again found. In addition ^ the spectra of the sedimentary soil 
was found to be distinguishable from the background as well as the serpentine 
soils studied on the San Francisco Peninsula.. 

The clustering program STANSORT developed by the 
.Stanford Remote Sensing Laboratory was then applied to the study areas \‘7ith 
significant success. 

2. 3. 5.3 AREAS STUBIED 

Two major exposures of serpentine rocks and soils mapped by the 
USGS, on the San Francisco Peninsula a were selected for study and are sho™ 

In Figure 2.151 Area 1 exposures consist of higlily V7eathered blue-gray serpentine, 
only a small percentage of which is outcrop , the remainder decomposed fragments, 
grading to a serpentine soil. These exposures are to the east of, and adjacent 
to, the Crystals Springs Reservoir segment of the San Andreas Fault Zone. 

They are surrounded by and occasionally penetrate the various rocks of the 
Franciscan assemblage. To a large degree the northern section of Area I is 
obliterated by housing developments and roads. Therefore, the study was focused 
on the southern section which is largely within the Crystal Springs watershed 
and is public land. 

The vegetation of this site, composed of annual broad leaved herbs and 
annual grasses, is readily distinguishable from the surrounding grassland 
on nonserpentine soil. It is marked by a different species composition, smaller 
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sise (height less than one foot), siiai’ser cover and earlier onset of 
senescence and di’ying. The dominant broad leaved herbes are Layia platyglossa 
(tidy tips), Qrthocarpus sp . (owl's clover) and Plantago erecta (California ; 
plantain) and the dominant grasses are Bromus mollis (soft chess) and Lolium 
mult if lor um. (ryegrass). 

Area II is approximately miles south of the Ciystal Springs Reservoir, 
again on the east side and adjacent to the San Andreas Fault Zone. Because 
of the housing developments and roads an open field area of roughly 60 acres 
at the south end was selected for study. The serpentine is heavily weathered 
and blue gray in color with only a small percentage of outcrop; the remainder 
decomposed fragments and serpentine soil. Interstate P80 transverses the south 
end of the area eiqiosing large amounts of fresh serpentine in the roadcuts. 

The serpentine vegetation of the Farm Hill Road site is clearly differentiated 
from the surrounding nonserpentine vegetation hy the same features that distinguish 
the Crystal Springs Road serpentin^> vegetation, i.e. , a different species 
composition, s mall er size, sparser cover and earlier onset of senescence. It 
is made up of annual broad leaved herbs and grasses and shares several species 
in common with the Crystal Springs site • The dominant plants are a grass , • 

Festuca sp . (fescue), and the broad-leaved herhs, Layia platyglossa (tidy tips) 
and Hemizonia sp . (tarweed). 

The third area studied consists of a sedimentary area located in the 
southeast quadrant of the Midway 7-1' topographic quadrangle (See Figure 2. 3. 5. 2) 
ThiS; area lies within the Lawrence /Livermore Radiation Laboratories Field Test 
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Site 300 and is largely a yearly controlled 'burn area. The study site is typical 
of the rolling grassy eastern foothills of ..he Coast Range. It is roughly 
600 acres in extent, crossing elevations varying from 1000 to l600.feet. The 
pediments are serai-consolidated sandstones with the outcrops again a minor 
pei*centage compared to the soils dei'iven from the sandstone. The vegetation 
is that typically described as a California valley grassland community, dominated 
by annual species of the grasses Broraus (bromegrass) , Festuea (fescue), Avena 
(oat) and others. 

2. 3. 5. 4 VISUAL STUDY 

A visual study of available ERTS imagery, both the individual bands and 
color composites covering the San Francisco Peninsula was accomplished. Also 
included was U-2 imagery taken during the ERTS Simulation Program. It was noted, 
in the ERTS frame date 6 October 1972, that a distinct dark gray pattern existed 
vdaich seemed to coincide generally with the Area I serpentine east of Crystal 
Springs Reservoir. (Fig. 2. 3.5.1) Study of the imagery before and after this 
date indicated that the pattern persisted with diminishing intensity back to 
26 July 1972 after which it could not be seen. The pattern was not evident again 
lontil 26 August 1973, at which time it was faintly discernible. The appearance 
and disappearance of the observed pattern seemed to correlate with tbe die-back 
and growth cycle of the grass in this area. 

Review of the ERTS color composites substantiated the above, with the 
pattern readily discernible at the dates noted, as a dark purplish tone. In 
addition, similar tones were also evident within Area II, south of Farm Hill Road, 
coinciding with the open field mentioned previously. 
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2 . 3 . 5, 5 


RADIANCE SPECTRA 


To study the possible uniqueness of the tones associated with the 
serpentine areas, the radiance values of ERTS-CCT pixels traversing these 
and adjacent areas where obtained and their spectra plotted. These pixel 
traverses, across the Crystal Springs Reservoir and Farm Hill Road areas, are 
indicated in Figure 2, 3. 5.1, 

Thble 2353- lists the mean radiance values, standard deviations and coefficients 
of variation relative to terrain types, across these traverses. Typical spectra 
are plotted in Fig. 2 Radiance throughout this report is presented as 

digital or word count levels. Should absolute value of radiance be desired 
conversion factors must be applied. . At this point no atmospheric corrections 
were made. Fig* 235.5 contains radiance spectra plots of the traverses indicated 
in Fig* 235. 1 The location of specific features was accomplished by means of a 
skewed pixel overlay’ of the proper scale and an ortho-photomap (l:2it000), as 
well as aerial photographs of the areas. 

It can be seen from examination of the spectra plots that the serpentine 
and grass areas as well as Interstate 280, water and the forested areas appear 
to have distinctive spectra. It is interesting to note that while the pixel 
spectra across the forested area in traverse AA are generally the same shape, 
peak values are evident at four points. The aerial photographs indicate that 
these coincide with the hilly terrain across which the traverse was made. 
Apparently, this effect is caused by the variation in sun angle due to hill slope 
The repetitiveness of the individual water spectra is also very striking. 
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' TABLE 2. 3. 5.1 - ^ REFLECTANCE STATISTICS PLOTTED FROM 
SELECTED ERTS PIXEL TRAVERSES AA, BB CC 

TREES WATER SOIL/GRASS I 280 



4 

■■ 1 

6 

1 

4 

5 

6 

1 

4 

1 

6 

1 

4 

5 

6 

1 

Mean 

19.33 

11.20 

20.88 

12.25 

19.14 

8.89 

5.67 

1.82 

26.53 

21.67 

20.88 

10.12 

38.86 

36.14 

32.4 

13.79 

Std Dev 

1.48 

2.09 

3.42 

2.72 

0.90 

0.82 

1.01 

0.57 

1.51 

1.48 

2.13 

0.86 

3.03 

2.88 

3.08 

1,53 

Coef of 

Var 0,08 

0.19 

0.16 

0.22 

0.05 

0.09 

0.18 

0,43 

0.06 

0.09 

0.10 

0.08 

0.08 

0.08 

0.10 

0.11 



Figure 2. 3. 5. 3 Radiance Spectra - ERTS Frame ID 1075-18173, 6 October 1972 
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iVlean radiance values from table 1 
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3*5.4 Radiance Spectra Traverses (By Pixel) ERTS Frame ID 1075-18173 







In traverse BB, the constancy of the pixel spectra of the serpentine 
soil/grass area on the east side of Interstate 280 as contrasted to the west 
can be correlated to the tones evident in the ERtTS imagery. The slight 
variability evident on the west side is attributed to variability in the soil, 
grass cover or both. It seems apparent that the spectra obtained is a function 
of the interaction of the soil and degree and type of grass cover which in turn 
is a function of the season of the year. 


2 . 3 . 5. 6 


SEASONAL REFLECTANCE SPECTRAL STUDY 


Based on the results of the visual and radiance spectra study it was 
evident that the serpentine soil/grass signature was unique at the 6 October 
1972 date. It also seemed possible that because the grass die hack in this 
part of Californ.Ia, was complete by this date, that the recorded spectra was 
essentially that of the soil. To substantiate the above, a systematic study 
of the soil plus grass interaction at the Crystal Springs area through the yearly 
cycle was instituted. Due to a fortuitous 15 acre grass fire, within Area I 
which occurred 1 July 1973 and easily seen in the ERTS imagery, it was also 
possible to include spectra of the devegetated burn area in the study for 
comparative purposes. This study was accomplished by means of four band spectra 
plots of the mean ground radiance values of the selected test areas and then the 
values normalized to band 4. In so far as was possible identical ten pixel 
areas within the following ERTS frames were utilized, covering an 11 month time 
interval : 
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ID 1075-18183 

6 Octaljei^ 1972 

ID 1165-18175 

4 January 1973 

ID 1185-18175 

22 January 1973 

ID 1291-18182 

10 May 1973 

ID 1309-18181 

28 May 1973 

ID 1345-18180 

3 July 1973 

ID 1363-18173 

21 July 1973 

ID 1399-18170 

26 August 1973 


In order to be able to compare results from the ERTS niultispdctral 
scanner data over these series of tapes, corrections were made for the . 
perturbing effects of radiation scattered by the atmosphere and. the variation 
in irradiance on the scene with solar zenith angle. These effects were 
removed by studying selected targets of low (zero) reflectance and high 
■ known reflectance (Honey and Lyon, 1974)* In the scene studied, a waste 
products treatment pond at an oil refinery near Suisun Bay, with bandpass 
reflectance of <0.5% in all four bands was utilized as the zero reflectance 
standard, A concrete parking aproii for aircraft at Moffet Field NAS 
California with reflectances of 27,8, 31*0,30*0 and 32,3 percent bandpass in the 
four ERTS channels was used for the high reflectance standard. The factors 
derived were applied as follows: 

p ^ Target Radiance-Dump Reflectance (Keas ) x Concrete 

arge >.e ec ance Concrete Radiance-Dump Reflectance (Meas) Reflectance 

(tleas) 


Results obtained are presented in Table II, Radiance plots as well as 
normalized reflectance are presented in Figure 2. 3, 5, 6,1. Study of this 
data reveals the following: 

2. 3. 5. 6.1 Interpretab ility of the four band spectra is greatly improved by the 
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Figure 2. 3. 5. 5 Radiance and Normalized Reflectance Spectra - Soil/ Grass (Crystal Springs) 


table 3 . 3 . 5 . 2 - ERTS-GCT RABIAHGE DATA 

GROUP STATISTICS ARD KOSII&LIZED EEFLEGTAH^^^^ 


ID 


ID 


ID 


ID 


1075-18183 CRYSTAL SPRIMGS 6 October 1972 



4 

■ 5., ■■■ 

6 

■-..7 

Mean 

25.40 

20.40 

20.20 

9.60 

Sta. Dev. 

0.84 

1.06 

1.55 

1.26 

Goef. of ?ar. 

0.03 

0,05 

0.08 

0.13 

Reflectance 

5.11 

6.26 

8.19 

13.97 

Horm. Refl, 

1.00 

1.23 

1.58 

2.73 


1183-18175 

CRYSTAL 

SPRINGS 

22 January 1973 


4 

5 

6 

7 

Mean 

16.80 

U.50 

17.10 

9.10 

Std* Dev, 

0.63 

1.35 

1.73 

0.99 

Coef. of Var* 

0.04 

0.12 

0.10 

0.11 

Reflectance 

3.01 

4.43 

13.65 

20.80 

Norm. Refl. 

1.00 

1.47 

4.53 

6.90 


1201-18182 

CRYSTAL 

SPRINGS 

10 May 1973 


4 

5 

6 

7 

Mean 

29.70 

26.20 

37.90 

21.10 

Std. Dev. 

1.95 

1.93 

4.86 

2.56 

Coef. of Var,. 

0.07 

0.07 

0.13 

0.12 

Reflectance 

4.85 

7.14 

14 . 45 

27.00 

Norm. Refl. 

1.00 

1,47 

2.98 

5.57 


1345-18180 

CRYSTAL 

SPRINGS 

3 July 1973 


4 

5 

6 7 

Mean 

37.50 

37.70 

44.20 22.60 

Std. Dev. 

1.18 

3.37 

2.25 1.35 

Coef. of Var. 

0.03 

0.09 

0^05 0.06 

Reflectance 

8.63 

11.19 

20.86 25.60 

Norm*. Refl, 

1.00 

1.30 

2.42 2.97 


TABLE 2 . 3 . 5 . 2 GOHTINUED 


ID 1345-18180 CRYSTAL SPRINGS 3 July 1973 

(Burn Area) 



4 

5 

6 

7 

Uean 

31.00 

25,40 

25.60 

10.60 

Std. Dey- 

2.45 

4.22 

6.54 

3.6 

Coef. of Var. 

0.08 

0.17 

0.26 

0.34 

Reflectance 

5.60 

6.03 

7.66 

11.38 

Norm. Refl. 

1.00 

1.08 

1.37 

2.03 

1363-18173 

CRYSTAL 

SPRINGS 

21 July 1973 



4 

5 

6 

7 

Mean 

35.10 

34.60 

41.70 

21.00 

Std* Dev* 

2.57 

4.40 

2.11 

1.49 

Coef *of Var* 

0.07 

0.13 

0.05 

0.07 

Reflectance 

7.29 

10.89 

16.80 

23.66 

Norm* Refl. 

1.00 

1.49 

2.32 

3.25 


1363-18173 

CRYSTAL 

Springs 




(Bum Area) 
4 5 

6 

7 

Mean 

Std, Dev. 
Coef. of Var. 

31.00 

3.06 

0.10 

25.50 

3.75 

0.14 

29.20 

9.08 

0.32 

13.50 

5.78 

0.43 

Reflectance 
Nona. Refl. 

5.36 

1.00 

5.98 

1.12 

10.74 

2.01 

14.82 

2.76 


1399-18170 

CRYSTAL 

SPRINGS 

26 August 


4 

5 

6 

7 

Mean 

31.50 

32.40 

35.40 

18,60 

Std. Dey. 

1.27 

1.65 

1.71 

0.52 

Coef. of Var. 

0.04 

0.05 

0.05 

0.03 

Reflectance 

5.96 

9.90 

18.30 

21.66 

Norm. Refl. 

1.00 

1.66 

3.10 

3.63 
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TABLE 2. 3. 5. 2 COHTIHUED 


1399-18170 

CRYSTAL 

SPRINGS 

26 August 1973 


(Burn 

Area) 




4 

5 

6 

7 

Mean 

28.10 

25.10 

24.60 

12.00 

Std. Dev. 

1.73 

3.00 

4.25 

3.13 

Coef. of Var. 

0.06 

0.10 

0.07 

0.26 

Reflectance 

4.20 

6.52 

7.93 

13.54 

Norm. Refl. 

1.00 

1.55 

1.89 

3.22 


ID 1075-18173 FARM HILL ROAD 6 October 1972 


4 

Mean 28 . 80 

Std, Dev. 1.62 

Coef, of Var. 0.06 

Ref iec tance 7 . 28 

Norm. Refl. 1.00 


ID 1075-18183 MIDWAY 


4 

Mean 32,60 

Std, Dev. 3,48 

Coef. of Var. 0,11 

Reflectance 9.49 

Norm. Refl, 1,00 


ID 1165-18175 MIDWAY 


4 

Mean 17,73 

Std. Dev. 2.15 

Coef. of Var. 0.12 

Reflectance 3.37 

Norm. Refl. 1.00 


5 

6 

7 

25.50 

28.30 

13.40 

1.35 

1.49 

0.97 

0.05 

0.02 

0.07 

9.12 

12.9 

19.5 

1.25 

1.72 

2.68 


6 October 1972 


5 

6 

7 

31.73 

30.07 

13.33 

4.77 

3.47 

2.06 

0.15 

0,12 

0.15 

12.58 

13.89 

19.39 

1.33 

1.46 

2.04 


4 January 1973 


5 

6 

7 

12.93 

22.20 

12.60 

3.31 

5.10 

4.36 

0.26 

0.23 

0.29 

4.93 

21.25 

28.80 

1.46 

6.28 

8,55 
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TABLE 

i, \5.2 

CONTINUED 


V/ 

ID 

1291-18182 

MIDWAY 


10 May 1973 




4 

5 

6 

7 



Mean 

34.80 

41.13 

52.73 

29.93 



Std. Dev. 

3.36 

6.53 

8.22 

5.54 

V 


Coef. of Var. 

0.10 

0.16 

0.16 

0.19 



Reflectance 

7.40 

14.15 

22.00 

37.03 



Norm. Refl. 

1.00 

1.91 

2.97 

5.00 


ID 

1309-18181 

MIDWAY 


28 May 1973 




4 

5 

6 

7 



Mean 

39.60 

51.47 

59.13 

31.33 



Std. Dev, 

2.87 

4.82 

4.55 

2.38 

o 


Coef. of Var. 

0.07 

0.09 

0.08 

0.08 



Reflectance 

10.30 

18.37 

26.07 

37.13 



Norm. Refl. 

1.00 

1,78 

2.53 

3.60 


ID 

1309-18181 

MIDWAY 


28 May 1973 



(Bum Area) 





4 

5 

6 

7 


Mean 

29.40 

27.53 

26.33 

11.27 

4 ■ ; 

Std. Dev. 

2.47 

'4.78 

6.10 

3.08 

W 

Coef . of Var . 

0.08 

0.17 

0.23 

0.27 


Reflectance 

5.20 

7.76 

9.66 

13.36 


Norm. Refl. 

1.00 

1-49 

1.86 

2,57 

o 

ID 

1345-18180 

MIDWAY 


3 

July 1973 



4 

5 

6 

7 


Mean 

41.40 

53.6 

58.53 

30.67 


Std. Dev. 

4.12 

7.55 

8.86 

5.01 

U 

Coef. of Var. 

0.10 

0.14 

0.15 

0.16 


Reflectance 

10.45 

17.84 

22.86 

35.16 


Norm. Refl. 

1.00 

1.71 

2.19 

3.36 



TABLE 2. 3. 5. 2 COHTIJIUED 


ID 1345-18180 


Mean 

Std. Dev. 
Coef. of Var. 

Reflectance 
Norm. Refl. 


ID 1363-18173 


Mean 

Std. Dev. 
Coefr of Var. 

Reflectance 
Worm. Refl. 


ID 1363-18173 


Mean 

Std. Dev. 
Goef. of Var. 

Reflectance 
Worm. Refl. 


ID 1399-18170 


Mean 

Std, Dev. 
Coef. of Var. 


MIDWAY 


3 July 1973 

(Burn Area) 



4 

5 

6 

7 

33.07 

33.00 

31.60 

14.93 

3.84 

6.46 

6.99 

4.20 

0.12 

C.20 

0.22 

0.28 

6.57 

9.22 

10.43 

16,51 

1.00 

1.40 

1.59 

2.51 


MIDWAY 


21 

July 1973 

4 

5 

6 

7 

46.93 

60.40 

64.27 

31.60 

6.P8 

11.84 

12.07 

5.93 

0 15 

0.20 

0.19 

0.19 

13.00 

21.43 

27.23 

36.27 

1.00 

1.65 

2.09 

2.79 


MIDWAY 


21 

July 1973 

(Burn Area) 




5 

6 

7 

33.47 

33.07 

31.33 

14.07 

3.23 

5.51 

6.35 

3.17 

0.10 

0.17 

0.20 

0,23 

6.50 

9.33 

11.29 

15.49 

1.00 

1.44 

1.74 

2.38 


MIDWAY 


26 

August 1973 

4 

5 

6 

•7 

43.13 

56.93 

63.07 

jO.87 

3.81 

5.43 

7.01 

3.96 

0.81 

0.10 

0.11 

0.13 

11.99 

21.25 

27.49 

36.7'6 

1.00 

1.77 

2.29 

3.07 



Reflectance 
Worm, Refl. 


TABLE 2. 3.5.2 COHTINUED 


(99-18170 

MIDWAY 


26 August 


(Burn Aroa) 




4 

5 

6 

7 

Meaii 

35.93 

37.27 

37.67 

16.60 

Std. Dev. 

4.61 

5,74 

7.54 

3.40 

Coef. of Var. 

0.13 

0.15 

0.20 

0.20 

Reflectance 

8.26 

12.16 

14.58 

19.20 

Norm. Refl. 

1.00 

1.47 

1.77 

2.32 






application of the atmospheric corrections and normalization of the data 
‘ to band 4. 

2. 3. 5- 6. 2 The normalized reflectance of the soil/grass is at a maximum 
(particularly channels 6 and 7) at the height of the rainy season, 22 January 

^ 1973; roughly twice as high as that during 6 October 1972, near the end of 

the dry season. 

2. 3. 5. 6. 3 The normalized reflectance of the soil/grass gradually diminishes 

S with the end of the rainy season and the entry into the’ summer dry-out 

period. 

2. 3. 5. 6. 4 In the bum area, on 3 July 1973, the normalized reflectance spectra 

O drops to a minimum value, possibly as a result of both the devegetation and 

the residue of the carbonized-grass remains. 

2. 3. 5- 6. 5 By 21 July 1973, the normalized reflectance spectra values in the 

0 - area have increased, probably as a result of the dispersion of the carbonized 

ash by the wind. They now approximate the values at 6 October, 

2. 3, 5; 6. 6 A slight increase in the normalized reflectance is noted at 26 
Ci August 1973, probably due to some revitalization of the grass in the burn area. 

2. 3. 5. 6. 7 The 6 October 1972 reflectances are very close to those of the bum 
area at 3 July and 26 August 1973. 

It would appear that based on the above, a strong likelihood exists 
that the reflectance spectra of 6 October is that of the serpentine soil with 
little or no reflectance introduced by the dead grass. 

A 10 pixel block within the Farm Hill Road test site> (6 October 1972, 
serpentine soil/grass) was also selected and the normalized reflectance 
^ spectra plotted in Fig 235.7. A strong correlation is found with the 
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reflectance spectra at Area I at the same time of year. Based on these 
results, it was decided to broaden the scope of the study somewhat to see if 
the same trends were obtained at a third site in which the terrain soil 
was of a different type and at which a similar burn situation existed. 

This test site (Area III) was located at the Lawrence/Livermore Radiation 
Laboratories Field Test Site 300 at which controlled burns were used to 
reduce the likelihood of uncontrolled fires as a result of explosive tests. 

It is on the east side of the Coast Range, approximately 15 miles east 
of Livermore. The terrain type had been mapped as marine sediments. In 
addition to another soil type, a comprehensive verification program of ground 
reflectance measurements of bare soil was to be accomplished at all three 
test sites. 

The results of the seasonal ERTS-CCT normalised reflectance spectra study at 
Area III (Midway) is presented in Table 2. 3.5.2 and Fig. 2.3.5S. Because of the 
size of this site, approximately 960 acres, a grid pattern, at intervals 
of 10 pixels, was utilized across the burn area to obtain the reflectance 
spectra data. It can readily he seen that the same trends exist, as follows: 

2. 3. 5. 6. 8 The normalized reflectance is at a maximum in the winter, 4 January 
1973 and gradually diminishes with the end of the rainy season and entry 
into the summer dry-out period. 

2. 3. 5. 6. 9 The 6 October 1972 reflectance spectra and that of the burn area 
are again comparable. 

A comparison of the reflectance spectra for Area III (Marine sediments) 
and Areas I and II (serpentine soils) indicate that they are considerably 
different, particularly in bands 6 and 7, thus leading to the conclusion that 
computerized classification could be applied. 
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The statistics relative to the ground reflectance measii Tements taken 
at the test sites are presented in Table III and Fig. 2, 3. 5. 7. An Exotech ERTS 
Radiometer and scaled down satellite geometry were utilized to obtain this 
data. Fig* 2. 3.5.7 compares the ERTS-CCT spectra with that obtained above. 
Study of this data reveals the following ; 

2.3.5.6.10 The ERTS-CCT normalized reflectance spectra for serpentine soils, 

at 6 October 1972 are almost identical at Crystal Springs and Farm Hill Road. 

2.3.5.6.11 The ground reflectance spectra obtained at Midway correlate very 
well with that derived from the 6 October 1972 ERTS-CCT data. 

2.3.5.6.12 The ground reflectance spectra for serpentine soil at Farm Hill 
Road compares very favorably with the CCT data. The ground spectra for 
the roadcut- and outcrop-serpentines, while comparable to each other are 
substantially different than that of the soil. It is believed that because 
of the small areal extent of the outcrops as compared to the soil and the 
limiting resolution of the ERTS system, the outcrops have little integrated 
effect and essentially only the serpentine soi,\ is detectable on the CCT data. 

2.3.5.6.13 The correlation of the ground spectra obtained at Crystal Springs 
with the ERTS-CCT data is not quite as good. It is believed that this is 
due to the infiltration of materials from adjacent soils derived from nearby 
Franciscan sandstone exposures. 

2.3.5.6.14 A significant difference is seen in both the ERTS-CCT spectra 

and the ground spectra for the serpentine soils at Crystal Springs and Farm 
Hil'> Road and the sediments at Midway. 

In general, from the season/ spectral study made and the ground measurements 
obtained, it can be concluded that the four band ERTS spectra for serpentine 
soils and sedimentary soils are sufficiently different from each other at the 
end of the dry or dieback season, to be distinguished from each other and the 
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TABLE 2. 3.5. 3 - CONTINUED 






FARM HILL ROAD (SOIL) 



4 

5 

6 

7 

Mean 

7.43 

11.33 

14.96 

16.97 

Std. Dev. 

0.45 

0.65 

0.58 

0.59 

Coef. of Var. 

0.06 

0.06 

0.04 

0.03 

Norm. Eefl, 

1.00 

1.52 

2.01 

2.28 


MIDV7AY (SOIL-GROUP 1) 



4 

5 

6 

7 

Mean 

15.09 

18.83 

22.95 

25.33 

Std. Dev. 

0.77 

0.98 

1-79 

1.68 

Coef . of Var . 

0.55 

0.05 

0,08 

0.07 

Norm. Refl. 

1.00 

1.25 

1.52 

1.68 


MIDWAY (SOIL-GROUP 2) 



4 

5 

6 

7 

Mean 

10.09 

12.03 

14.22 

16.70 

Std. Dev. 

0.50 

0.62 

0.93 

0.59 

Coef, of Var. 

0.05 

0.05 

0.07 

0.04 

Norm. Refl. 

1.00 

1.19 

1.41 

1.66 


Q 


MIDWAY (SOIL-GROUP 3) 


Mean 

15.79 

18.81 

22.02 

23.97 

Std. Dev. 

1.66 

1.43 

1.78 

3.38 

Coef .of Var, 

o.ii 

0.08 

0.08 

0.14 

Norm. Refl. 

1.00 

1.19 

1.39 

1.52 


O 
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background by computerised clilsterlng^^. 


>.3.5.7 


CLASSIFICATIOlg lECHMIQOE 


The classification procedure utilized to demonstrate the uniqueness 
of the soils spectra studied is an interactive program package called 
STAMSORT, developed at the Stanford Remote Sensing Laboratories (Honey 
et al, 1974). This system provides an extremely rapid, flexible and 
low cost tool for scene classifcationn. It is a non-statist ical^ 
uttsupervised classification technique in which the data is split into its 
distinguishable groups with no prior knowledge of the groups. The 
primary classification procedure, utilizes a search, with variable gate 
widths, for similarities in the normalized or un-normalized digitized 
spectra. Results obtained from the application of STANSORT are presented 
in Fig. 2.3.5. 9 and 2,3.5.10, the Farm Hill Road and Midway test sites 
respectively. Training was accomplished on the serpentine soil spectra 
of Farm Hill Road which is indicated by the clustering in Figure 2. 3. 5. 9. 
Continued search for the serpentine spectra at Midway (the sedimentary 
test site)using identical classifiers, as expected, reveals its complete 
absence except for the single pixel indicated in Figure 2.3.5.10. 
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Fig* 2# 3- 5* 9 Midway - Results of Classif icatioa Program 
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2. 3.5. 8 

As a result of the foregoing study in the San Francisco Bay and adjacent 
Coast Range grassland areas the following may be concluded: 


2. 3. 5.8.1 ERTS soil/grass four band spectra are in fact essentially soil spectra 
at the end of the dry or grass dieback season. 

2. 3. 5. 8. 2 The ERTS font band spectra obtained is a function of the interaction 
of the soil and degree and type of grass cover which in turn is a 
function of the season. 


2. 3. 5. 8. 3 A strong correlation exists between ground measured reflectance spectra 
and ERTS four band spectra for both serpentine and sedimentary deriven 
soils. 


2. 3. 5. 8. 4 The ERTS four band spectra for serpentine and sedimentary deriven 

soils are sufficiently different from each other and their background 
to be classified by application of Stansort the SRSL interactive, 
unsupervised classification program. 
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2,3.6 CORRELATION BETWEEN GRODIjtD METAL ANALYSIS, VEGETATIOM REFLECTAMGE, 

AND BETS BRIGHTHERS OVER A MOLYBDENUM SKARH DEPOSIT, PINE NUT MOUNTAINS, 
WESTERN NEVADA 

2. 3. 6.1 abstract 


In a cooperative study with TJSGS personnel, it has been possible to detect 
a 1.5 by 1 mile anomaly on ERTS CCT data directly, in the pine-covered 
mountains of western Nevada. This anomalous area is about 3-5 times larger 
than that of the known geobotanical anomaly which lies centrally within 
the area. The site has been studied on the ground and bi-directional 
reflectances (relative to BASO 4 obtained for 40 trees, using both in-vivo 
techniques (similar to cherry picker operations) and field determinations of 
cut branches. The anomaly can be seen best by color transparencies made 
from 5/4, 6/4, 7/4 ratioed digital data, the 3 ratios each being coded by one 
of 3 colors (blue, green, and red) . 

Field reflectance measurements of three modes were made, using EXOTECH ERTS- 
type radiometers — cut branches, and viewing the trees both from vertically 
above, and horizontally. Each tree, either a Pinon pine or Juniper, was one 
previously marked by the US6S, who provided the molybdenum analyses of stems, 
twigs and needles (leaves) . In addition sagebrush and bitterbrush shrubs were 
measured together with their background soils and rocks. 

The correlation between Mo and Juniper leaf reflectance was positive, and 
significant at the 99% level (Channel 7 brightness) agreeing with the visual 
observation that even at values in excess of 500 ppm Mo in leaf asb, she 
junipers were healthy. With Pinon Pine however the correlation with leaf 
(needle) reflectance was negative but significant at the 97% level. The pines 
showed significant morphological change (needle lc.;s, profusion of twiggy stems, 
and brittleness of branches) correlateable with mineral uptake of Mo. 

Collectively all 62 tree samples showed a negative correlation with molybdenum 
significant at the 97% level. 

Using unsupervised clustering techniques on CCT taped data (STANSORT program) 
ERTS spectra could be extracted for the total anomaly area, which were used to 
locate similar areas to the south, near Double Springs Flat. Field checking 
located weak gossan mineralization in the bleached andesites there. 

Continuing field studies are aimed at specifically identifying the cause of 
the ERTS anomaly — is it tree vigor, tree species, tree spacing, or sagebrush/ 
soil ratio which can be observed from space over this sicarh zone. 


2. 3.6.2 Table 2.3.6. 1 which follows summarizes the results we obtained from 
several sets of field measurements in the Alpine-Cherokee Miner-, area. 
Pine Nuts Mountains, Nevada 

The detailed descriptions are reproduced in Appendix D. 
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COEEELATION BETWEEN LOG Mo (ppm) , Mo (ppm) WITH VARIOUS 
REFLECTANCE PARAMETERS 


Log Mo (ppm) 07') 


All 

Branches Pine Juniper 
N=62 N=29 N=33 


13 (-.40) 


Mp (ppm) (#4) 


All 



. 

' Branches 

Pine 

Juniper 

N=62 

N=29 

N=33 


:o 

-.13 

.6 

-.20 


-.16 

)4 

1 

O 


.13 

(.42)** 




■■■■■■ 



Signlficanlj 
& at l% t 
(M at 5%| 

.325 .45 

^ A5 \ 

,33 .45 

r ,45 i 

1 

.25 .35 

’ . 35 j 

.25 .35 

r .35 1 

1 




* 

t 


SUMMARY 

Pines ; Negative, Juniper : Positive Reflectance Correlation with Mo 

Significant at 


1% level 
3% level 


/Juniper BP7 vs. Mo (ppm) 0.42 
^All Branches BP 7 vs. Mo (ppm) 0.35 

Pino R4 vs. Log Mo (ppm) -0.40 


/■All Branches BP7 vs. Log Mo (ppm) ,30 
5% level JAll Branches R7 vs. Log Mo .25 

I All Branches R7 vs, no .26 














2.3*7 APPLICATION OF DIGITAL SNOW MAPPING WITH ERTS-*! DATA, USING THE 
STANSORT, IMAGE PROCESSING SYSTEM 


2*3. 7.1 

The STANSORT image processing system (Honey, je^. al* 1974) was 
used in a digital snow mapping experiment v^ith test . Ites in the 
Windriver Mountains (Wyoming) , Using specific density slicing and 
cluster techniques it was possible to separate the following ground 
cover classes: dry siiow,**metamorphi.c'^-(partly"melted) snow, snow covered 

forest, snow/non-snow transition zone, and bare forest/other vegetation .Especially 
in heavily forested areas, where problems in manual image interpretation 
arise, the "snoxjline" - or better "transition” zone ^ could clearly 
be detected. The STANSORT system proved to be a very useful, easy 
to operate, lowcost interactive tool. 


2. 3. 7. 2 INTRODUCTION 

With the first images from TIROS the potential for large-area 
snow mapping was recognized. With the ESSA-^APT System many nations had the 
possibility of instantaneous satellite data for metfeorlogical purposes, 
and in some of the countries these images (with resolutions in the 
order of kilometers) led to ii.teresting snow mapping studies. 

Basically the following main techniques were developed and applied: 

For flatlands - the US/NOAA method of computing 5-day composite minimum 
brightness charts were very successful. In mountainous terrain however 
up to now, digital methods were of less success. Therefore manual 
photointerpretation techniques were used. In the United States, snowline 
drawings v7ere transferred onto topographic maps to locate the altitude 
of snowline. The technique kno^m as SWISS uses transparent isolevel 
contour maps. By overlaying those onto the satellite imagery, a best 
fit method was used to estimate the snowline altitude. Very encouraging 
results have been achieved using this method. 

With tfie advent of RRTS we were able to work with high resolution 
multiband imagery. The same image interpretation techniques were applied, 
but soon two main problems showed, 

2. 3- 7.2.1) It is vei:y easy to map the snowline as long as it is above the 
treeline. Detecting the limit of snow cover in forested areas 
however is very difficult. 

2- 3. 7* 2.2) Secondly, the improved resolution resulted in much more detail 
in the images - details no one had to deal with in satellite 
imagery previously - making the job of snc^^ 7 mapping very time 
consuming, almost impossible for larger areas. 

In the Runoff Prediction model being developed at the Hydrology Branch 
of Goddard Space Flight Center, precise input data concerning the areal 
extent of the snow cover is essential. Therefore the question arose, 
whether existing digital image processing techniqes can be applied to an 
automatic detection and classification of snow covered area. 

In coordination with the above mentioned branch, the Windriver Mountains 
in Wyoming were selected as testsites. Image interpretation studies are 
conducted in the same area, and good "ground truth" in the form of high 
altitude U-2 photography is available (RC~8 )- 
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2. 3. 7-3 


THE STANSORT SYSTEM 
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The STANSORT system allows to read directly the standard NDFF 
ERTS 9- track computer compatible tapes, and to apply the following 
functional handling procedures to the pixel data: Smoothing, 

ratioing, edge detection, normalized (as well as unnormalized) 
clustering, removing of atmospheric effects, calibration, shadeprinting, 
extraction of data values, histogramming* etc* (see Figure 1). The 
interactive access to the PDP“10 computer is achieved through a keyboard 
and instant control over the extremely fast operating process is 
possible via a black and white CRT display. A very fast primary evaluation 
of test results is thereby possible, and parameter changes may be applied 
very quickly. The system operation can be learned very quickly (l^-B hours) 
making it a very valuable tool for the discipline-oriented user-investigator 
who often is not a specialist in digital image processing. 


^^Since the time of writing this report (November 1974) "smoothing" has 
been taken out of STANSORT and replaced by "debanding" (6-line repeat 
bands) and "deconvolution" added. R.L. 


2. 3. 7 .3. 1 The Procedure for the Snow Mapping Job 


Step 1: 


rtep 2 : 


Step 3; 


Step 4; 


Step 5: 


Step 6 : 


Test area coordinates are calculated from the 1:1,000,000 
scale image and directly introduced in the system, the 

region was controlled by viewing the sliadeprint on the CRT 
display, and histograms were produced. 

Optimal shadeprints of the 68X60 pixel area (roughly 
26 square km* or 10 square miles) were created in all 
bands. Careful comp irison of the shadeprints and the 
histograms led to the conclusion to drop Channel 6 from 
further investigations due to the high level of banding and 
noise. This option is available in the "un-normalized" 
clustering" step. 

Establishment of test sites for each expected ground-cover 
class, e.g. , dry snow/ iiietamorphic snov7/forest with snow/bare 
forest 4- other vegetation. 

Special shadeprints of ERTS band 7 were made with density 
slicing of the saturation level 63, by which the dry snox^ 
areas could be obtained. 

Un-no rmalized cluster with Bands 4, 5, and 7. For bands 
4 and 5 the upper values of ranges were set to 126 to exclude 
all snow from the clustering. With this step, and 

combining with the special shadeprint, the metamorphic snow 
area could he determined. 

By changing the tolerance levels (gates) in the STANSORT 
clustering algorithm, different printouts were obtained. 

Careful comparison of these with the test field ground data 
led to the final application of the appropriate tolerance 
level for this specific application. 

Application of the same density slicing and clustering parameters 
to neighboring 68X60 pixel areas until the project area was 
covered* 
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Discussion of the Procedure 

2. 3. 7. 3. 2 Steps 1-3 do not offer problems, when either the coordinates 

(ERTS image coordinates) are known of tlie project area, or enough detail 
in the shadeprinted scene allows for easy location. 

The density slicing in Step 4 requires an experience factor 
knowing that melting snow has a lower brightness in Band 7 than the 
dry snow which reaches saturation. Previous work with other image 
enhancing- and slicing-equipment showed the usefulness of this step. 

In Step 5 the un-normalized cluster proved to be advantageous 
relative to the normalized cluster* approach because of too much loss 
of data variability in normalizing. 

*A11 channels are normalized to Channel 4, which tends to remove 
sunlit-versus sun-shaded-slope problem. 


The setting of the tolerance gate was the most difficult step in 
the whole procedure. Each intermediate result with a specific tolerance/ 
gate setting, had to be carefully compared with a big number of tests! tes 
before rejecting it, or finally approving and apply it to the full project 
area. Good ground data is very essential In this evaluation step. Once 
the parameters were knot«i, they could be saved in the computer memory 
and used for adjacent 68X60 pixel areas. 


|; Evaluation of the STANSORT System= relative to others . 

In the comparison with results obtained by the LARSYS Ver.3 
classification, the STANSORT cluster/slicing showed almost identical 
i2 distribution of classes. The amount of time and effort hov;ever was 

considerably smaller than with LARSYS. On the other hand, the control 
and final accuracy achieved with LARSYS 'supervised technique cannot be 
reached. There is a tradeoff between how fast and how accurate your 
final result should be. Figure 2 shows a direct comparison of hand 
colored samples of the results obtained with LARSYS and with STANSORT 

Positive Points : 

- easy to learn operation of the system 

- very fast cluster algorithm, makes the system comparatively 
inexpensive. 

- B/W CRT display for fast interaction by visual checkout of 
results. 

- easy to operate special features, like edge detection, simple 
ratioing, removing of atmospheric effects etc. 

Negative Points ; 

“ 68X60 pixel area limitations on one run 

- limited statistics 

- no training/ classify functions, except by using a selected 
tolerance set to classify adjoining areas. 


S output. 

2. 3. 7. 3. 3.1 


2 . 3 . 7 . 3 . 3. 2 
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2. 3. 7.3.4 Sufifiestions for Further Develppments 

'2. 3. 7. 3. 4.1 Increment feature : The location of a feature in a scene is 

sometimes difficult when looking at single pixel resolution for a 68X60 area. 
In a first step, to get a better overview, every 2nd or even 10th or 
higher number of pixel sample interval could be used, thereby producing 
a smaller scale. Secondly, the 1^ pixel resolution is not needed in some 
cases, esp. in large area mapping studies, thus computer time could be 
saved. , 

2. 3. 7. 3. 4. 2 Group cluster feature : To be asked by the computer/ system when 

all cluster parameters are established, and an interpretation of the differeiit 
clusters is achieved. In the following example 14 clusters had been built 
up, and nov7 we wish to merge several into three groups. 

CRT Question , Answer 

- "do you want to group clusters?" , . • • *Y " 

- "how many groups (classes) do you wish to create?". .... *3" 

- "assign symbol to group 1" . 

- "assign clusters to be grouped in group 1 *A,B,E,K" 

- "assign symbol to group 2" 

- "assign clusters to be grouped in group 2" *C,H,I 

- "assign symbol to group 3" 

- "assign clusters to be grouped in group 3" *D,F,G, 

“ "do you wish to shadeprint grouped cluster result?" *Y 

and so on 

This feature would allow one to make a sort of a classification, and would 
be a great help in simplifying the display and in the final evaluation 
of the result. 

2. 3. 7. 2. 4. 3 Statistics feature ; A coincident spectral plot showing the 
location (spectrally) of the clusters in the different wavelength bands, 
would help the interpretation and grouping steps considerably. Furthermore 
some statistics would be desirable on number of pixels grouped in group 1, 

group 2 etc as well as a table showing the distribution of the 

clusters within the groups (percentage of pixels of original cluster A 

in group 1, cluster B in group 1 etc..). The tabulated output of the 
spectral values (digital counts) assigned to each symbol is useful but 
additional statistics would be welcome. 


A.-. 
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2. 3. 7. 2. 4. A ' I^rge area feature ; In addition to the increment feature 

(par. 1) , sometimes big areas have to be looked at or searched for 
a specific signature, and the stepwise procedure by adding 68X60 
pixel areas until’ the whole area is covered is too time consuming. 

Therefore a special shadeprint feature - be it for raw ERTS data, 
for cluster/grouping results, for edge detection or anything else - 
could be desirable (see example). 

CRT Question Answer 

- "do you wish to apply the same parameters 

to another area? *Y 

- "type in title for output" Stanford Land Use Map 

“ "after next type in you may have to wait a while 

how many scanlines down to your area?". , . „ . . . *1000 

- "how many lines are in your area?" *1300 

- "what line increment do you want?" *5 

- "how many pixels in to your area from tape edge?" ... *200 

“ "how many pixels wide is your area"" *1500 

- "what column increment do you want?" *5 

and so on ... . 

..... .followed by; 

A confirmation of area coordinates, and a calculation of the 
time involved in waiting for a larger output result, which will 
for the lineprinter be divided in stripes, 68 columns wide, as 
long as stated. Mounting of those stripes would be considerably 
faster than mosaicing the 68X60 rectangles. 

2, 3. 7. 2. 4. 5 Special histogram approach for shadeprinting ; In some cases an 

optimal shadeprint is wanted with a non jLinear scale according to the 
distribution in probability space. As the histogram feature is already 
implemented in the system, the histogram-controlled scaling of the 
shadeprint steps would he easy. 
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Figure 2. 3. 7. 2 Comparison of interpretation samples made 

by using STANSORT (left) and LARSYS Ver. 3 
(right), strip size roughly 4 km X 1.14 km. 
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2 , 3 . 7 . 3 . 5 Further Studies Using Different Processing Systems 

As mentioned In the title, the use of the STANSORT system in the 
- application of ERTS data to digital snow mapping is a part of a larger 

® project, involving the use of also LARSYS Ver. 3 and the General 

Electric Iraage-100. Tlie results achieved and the efforts that had 
to be made when using these different kinds of systems will be compared 
and carefully investigated in a forthcoming report issued from the 
Department of Geography, University of Zurich. 

.•'I.--. 
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2.4.1 HEW TECHNOLOGY 

G ■ 

2.4,1. 1 Hardware 

2. 4. 1.1.1 Development of a projection- type densitometer for radiance 
measurements from ERTS ICS 70 oint transparencies (see Project Al)» 

2. 4. 1.2 Software 

Q 2.4.1i2.1 STAMSORT . a fully interactive program for a PDP-lO con^uter 

with CRT displays ; for processing ERTS CCT Tapes (see Project B5) 

2. 4. 1.2. 2 IMAGE , a program for a PDP-lO computer, to produce tapes of 
250 X 300 pixel areas for use on a DICOMED imaging system, with 
geometric corrections of ERTS distortions (see Project b4). 
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2.4.2 SIGNIFIGAm: RESTILTS 

2. 4.2.1 SCIENTIFIC 

2. 4. 2. 1.1 A far greater appreciation of the seasonal variability of vegetation 
together with the technological abilities to measure these effects 
either stationary (on the ground) or mobile with a truck mounted 
bi-directional reflectance system (see Projects C2.1, C2*B, C2.5> 
C2.4 reports). 

,1.2 Correlation of leaf reflectances for Juniper and Pine, growing in 
Mo-bearing mineralized soils, with ground ERTS- type spectral data, 
significant at the 1^ and 5^ levels, respectively. (See Project 
C2.U, report section.) 

2.4. 2.1.3 ERTS image (Aug, 2h, 1975) covering W. Nevada and the 

copper mines of Yerington was made I5 days after a SKYI4AB BB5T 
underflight obtained high-altitude photography (at 20 km) of the 
same area. A color IR. (false color) image of km made by our 
IMAGE program (Ch. 5 6 and 7) shows a very high degree of 

spatial agreement with the airphoto. Both the airphoto (made at 
0830 local time) and the ERTS image (made at O903 local time) 
show the same topographic relief effects due to the sun angle. 
Significantly^ a color-ratio-composite image made from h ^/k + 
(Blue) + Gh 6/h (Red) + Ch j/h (Green) shows greatly lessened 
topographic effects (from ratioing)^ but the tone patterns now 
closely resemble the published geological maps, specifically, 
precisely localizing three rock types, and differentiating 
'^alteration zones” within three others. The best single ratio for 


this purpose if Ch 6/lt- (see Project Bl|-; and C1.2). This high oorrela 
tion of ERTS, EB57 sir photograph and geological mapping is most 
important and is being further studied with the geological personnel 
of the copper mine. 


I 

2.4.3 

& 

2.4. 3.1 

Successful transfer of the technology of using ERTS tapes, to the 

g 

classroom, with both undergraduate and graduate geology students using 
the STANSORT interactive program in their regular studies of environ- 
mental geoscience monitoring the San Francisco Bay lands (AES 153)- The 
program has now been used also for two successive years in the Remote 
Sensing and Exploration classes in the Department of Applied Earth 
Sciences at Stanford (AES 295)* cases the geographic, map- 

“■ like data displays proved very easy for students to comprehend and use. 

The interactive mode of data processing, based upon their matching the 
computed output (cluster classifications), with their maps as photos 
■ - and field data, has shown itself to be very simple for them to operate. 

2. 4. 3. 2 IHTERHATIOMAL 

" Following visits to Stanford of several foreign geographers and 

geologists (students and professionals), the STAHSORT programs have now 
been installed (with transliterations to other styles of conq>uter) in 
Zurich (Switzerland), Trondheim (Norway), Hanover (W. Germany), Madrid 
(Spain), Sydney (NSW, Australia), and Perth (W.A., Australia). Con- 
tinuing research programs (of varying extents) exist with all of the 
^ above groups, further showing the technological transfer. (See Pro- 

jects A2, B3> Cl.O, C2.2, C 3 .O Summaries.) 

□ 


TECHNOLOGY TRANSFER 

EDUOATIONAL - STANFORD UNDERGRADUATE AND GRADUATE GLASSES USING ERTS 
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Abstract 


Two targets of measured, constant reflectivity in 
the area of San Francisco, California are studied. 
The first standard, a waste (carbon black) treat- 
ment pond at an oil refinery near Suisan Bay, 
having an area of approximately 0.3 square miles, 
(or 215 pixels), and bandpass reflectances of 
<0*5% in all four bands, is assumed to have a 
zero contribution to the radiance recorded by 
ERTS* The radiance observed then arises entirely 
from atmospheric scattering* The variation in 
these radiance values as a function of solar zen- 
ith angle has been analyzed. 

A second target, a concrete parking apron for 
aircraft at Moffett Field, California, assuming 
that it’ remains dry during the period of study 
has constant reflectances of 27.8, 31*0, 30.0, 
and 32.3 percent bandpass reflectances in four MSS 
equivalent chmmcls. Using these values, the 
radiance observed by ERTS may be corrected for the 
atmospheric contribution, and thus values for the 
radiance from the target may be calculated. These 
values were studied as a function of solar zenith 
angle. 

1. Introduction 

In order to be able to compare results from 
ERTS MSS data over a series of tapes, the per- 
turbing effects of a variable contribution due to 
radiation scattered by the atmosphere Into the 
detector £ii*lcj of view, and of the variation in 
the Irradiance on a target with solar zenith 
angle, must be eliminated. These two effects may 
be compensated for, or entirely removed, by 
studying selected targets In a scene, one (or 
more of low (zero) reflectancri, one (or more) of 
high, known reflectance. In some scenes, however, 
suitable reflectance targets may not be obtained. 
When this occurs, atmospheric modelling must be 
employed to arrive at some values for the atmo- 
spheric scattering contribution, and for the 
irradiance on the scene. 

The technique of using standard targets 
within a scene was applied to a specific scene 
which contains an area of measured ref lect:*vity*- 
the Stanford Grassland testsite. 

To increase the contrast between targets and 
background in a UVNDSAT-1 image scene it is 
necessary to attempt to remove the effects of 
atmospheric scattering. Several methods have been 
used generally based upon the scheme of ”dark-area 


subtraction’^ -some area of very low brightness is 
observed on the LANDSAT image. It should be black 
but it Is not and has finite values of brightness. 
These resulting radiance values are presumed to 
have arisen solely from the atmospheric column. 
Cloud shadows were used by Goetz v 2) and by 
Vincent C4)^ hut cloud shadows do not necessarily 
represent the darkest areas in any frame, parti- 
cularly in MSS Channel 7. 

We searched for suitable dark areas near 
Stanford in the San Francisco frames, specifically 
using Channel 5 for our search. We selected Ch 5 
because all water over 3 cm depth Is black in Ch 
7, but in Che same locality may have values as 
high as 30 in Ch 5, if the water is shallow with 
a sandy bottom, or is muddy. A spectrally black 
target should appear dark grey in all channels* 

The area we located, which was black on all 
channels, occurred on the south shore of the 
Sacramento River just east of the Carquinez 
Straight Bridge about 47 km (25 miles) NNE of San 
Francisco. Inspection of the local topographic 
maps (Port Chicago, sheet) Indicated that it was 
a part of an oil refinery complex (Phillips 
Petroleum-Port Chic ago/ Avon Plant) and with a^ 
telephone call we ascertained that the 0.8 km 
area (0.3 sq. mi.) was a dam filled with carbon 
black waste product. Figures la and lb are 
1:250,000 scale enlargements of this area and 
Figure 2 Is the maximum possible enlargement of 
Channel 5 (about 1:32,000 scale). 

Dl-dircctional reflectance measurements made 
on February 1, 1974 using two EXOTECIl-100 ERTS 
radiometers (15“ FOV, relative to BASO. powder) 
showed that the reflectance of the carbon product 
was bciow 0.5% in all four LANDSAT channel band- 
passes. 

We then proceeded to find a high refloctmict 
target, so th.nt we could calculate reflectance ot 
terrain directly from LANDSAT data, using a zero 
reflectance (carbon black dump) and the high re- 
flectance targets, and linearly interpolating 
between them. 

Initially we tried to use the large dumps of 
salt harvested by evaporation on the shoreline of 
San Francisco Bay, near Redwood City. These dumps 
are several pixels in extent, but their bright- 
ness exceeds the 127 count number (7-bit) allowed 
in the dyjifimlc range of the LANDSAT scanner tele- 
metry. Often in the summer periods all four MSS 
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data channels were saturated making a target as 
bright as salt of little use. (Similarly many 
cloud-tops and snow saturate the system making 
their mutual discrimination very difficult). 



Lg. 1. Enlargements at 1:250,000 scale of LANDSAT 
iMges for 1075-18173. Fig. la is Channel 5; lb 
1b Channel 7. White arrows noinl at the Carbon 
Black dump at the Thillips Tetroleum Avon Plant 
(Refinery) here used by us as a "zero- reflect- 
ance" standard. Letters, B^Benecla Bridge; C» 
Cnrqulnez Bridge; S-motliballed ships in Suisan 
Bay; and V-Vallcjo City. 


We decided to use the large expanse of 
concrete used as a parking apron for U.S. Navy 
aircraft at Moffett Field, Mountain View, only 10 
km (6 miles) ESE of Stanford. Similar bi-direct- 
ional reflectance measurements with the same units 
gave the following reflectance values, relative 
to BoSO . . 

4 | 

Concrete Apron, Moffett Field (N«15 measureme nts) 


C hann el ^ 

Ave.Rofl.(x) 27.88 
1 a 2.71 

Cov (o/x) (0.10) 


Both the carbon black 
sea level. 


1 6 7 • 

31.0 30.0 32.1 

3.42 3.53 3.14 

( 0 . 11 ) ( 0 . 12 ) ( 0 . 10 ) 

id concrete target are at 


Using the 14 LANDSAT CCT tapes of the 
Stanford area (San Francisco Bay) we extracted the 
average brightness (at the satellite) for eacli of 
the two standard localities - carbon black and 
concrete. Great care was taken not to use edge 
pixels, and just to take the central portion with- 
in each area to avoid effects of terrain contam- 
ination, and also forward scattering by other 
adjacent terrain. Table I lists the tapes used, 
together with tlieir dates, solar elevation and 
zenith angles for the 24-month period of tills 
study. 



Fig. 2. Maximum enlargements (scale =* 1:52,000) 
of the same Carbon Black area on Channel 5 image, 
as in Fig. la, showing the large area contained 
in an earthen dam near a small river. 
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SUMMER 


©CONCRETE 

•CARBON 

APUMICE 


;s. 3~6. Four plots of normalized LANDSAT (Channels 4»5»6 ind 7) bandpass data for 17 overpasses. Tor 
each plot the ordinate is brightness (in DN) , over brightness for the Winter scene (DN,min); the 
abscissa is cosine zenith angle of each pass, divided by the cosine of the maximum zenith angle 
(winter). For full adherence to the cosine-law all points should lie on the thin dashed line (1:1). 
Open circles arc for the concrete at Moffett Field; solid circles for the Carbon Black; triangles for 
the Mom* l.nke pumice snnd. Thick solid lines are least-squares fitted to the carbon and concrete 
data, d.isii-dot for the Mono sand. Krror bars denote tiie l-sigraa limits, where N-50 (carbon) and N**15 
(concrete). Channel A and 7 plots sliow departures from the 1:1 relat ioiiship. 






Table IX contains all tile measured LANDS AT 
data, arranged in increaaiag zenith angle (rouglily 
Winter to Summer). Mean brlglitnenses and standard 
deviations are calculated for eacli channel and 
date. Values are given in digital numbers (DN) 
directly off the bulk CCT tapes. This table also 
contains data for two sites at Mono Lake, Caiifor- 
nla where one of us (Gary Ballew) has similarly 
measured a dark and a light target on the ground, 
and extracted their digital values from the 3 Mono 
Lake tapes (referred to in Table I). 


Tiiu carbon black target (solid circles) shows 
!ar less relation with the 1:1 cosine- law line, 
particularly with Channel A data, wliich are some- 
what constant r nardless of die zenith angle. 
Tliere are Iiiglier values In the other channels 
(particularly Channel 6) but the data points are 
mucli more spread. In Channel 7 the^baslc DN 
values are so low that a change of -1 in the A-D 
decision (at the spacecraft) creates a^high 
standard deviation (up to AOK cov., c/sc). Errors 
are high. 


Study of the Table II listings suggest a 
pattern of relationship with zenith angle, so we 
calculated their cosine Z values. To make the 
plots more clear we normalized the brightness 
data to that of the maximum zenith angle (winter 
data) and compared these ratioed values with a 
similar ratio of cosine zenith angle/cosine zenith 
angle maximum. From these normalizations Table 
III was produced, and the values were also graphed 
at Figures 3,4,5, and 6 for Cliannel 4 5 6 and 7 
data respectively, as seen above. 


2. Analysis 

From the Tables and the graphs It is quite 
clear that for the bright target- concrete (open 
circle), in all channels, there Is a high correla“ 
cion between cosine zenith angle function and that 
of their LANDSAT brightnesse,'-', iiuJicating that 
the change In brightness is directly related to 
the elevation of the sun. For carbon black there 
is a different relationship (witli much more spread) 
although some cosine-law aspects are preserved. 
One-sigma bars are shown on some of the points. 


LEAST-SQUARE S FIT NITH CORJtElAlTON AND ELOPE 



Ch 4 

Ch 5 

Ch 6 

Ch 7 

Concrete 





R 

0.987 

0.991 

0.993 

0.988 

Slope 

.86 

.99 

.98 

1.00 

Mono Pumice 




R 

0.981 

0.980 

0.992 

0.986 

Slope 

.85 

.96 

.87 

.99 

Carbon Black 




R 

0.936 

0.890 

0.883 

0.796 

Slope 

.39 

.45 

.64 

.82 


In detail, the IrJ straight line rel*jtionship 
for the bfight target is most clear In Channels 
5,6 and 7 data, but with almost ail Channel 7 
values lying slightly above tl»e 1:1 line, whereas 
Channel 6 closely approximates the cosine zenith 
angle line. Channel 3 data lie slightly below 
the line. 

Ciiannel 4 data are tiiucli more spread, parti- 
cularly with the higher sun (summertime) dates 
which all lie below the line. The concrete is less 
bright than one would expect directly from the 
sun elevation. Tills lowering (and probably also 
the spread) of Channel 4 data no doubt are related 
to the Increased smog and alr-pullution over the 
San Francisco Day in the suitimur months. Air 
pollution levels are quite low in the winter 
(when also a much greater horizontal vlsabiiity 
can be documented). 


Some yearly variation in the values can be 
tioted, particularly in the 1974 tapes (ID: 1687 
and 1669) which show wildly high values for 
Channels 5,6 and 7. This may be drift in the 
satellite system, or a calibration change in 
ground data-processing. It may also be a change 
In the surface conditions of the dump Itself 
(albedo change, etc.), although the reflectance 
measurements we made at the refinery were made 
on February 1, 1974, 

3, Mono Lake Standard Terrains 

Ballew^^^ has measured two natural soil 
materials at Mono Lake, California, about 290 km 
(180 miles) duo east of San Francisco, and calcul- 
ated their bi-directional reflectances for a wide 
series of sun elevations. His data for a dark 
target (basalt lapilli) and a liglit target (pum- 
ice sand) have the following values, relative to 
Fiberfrax; 


Average Bl-di rectionai Reflectance ( ) 


Dnrk Target 

Cli 4 

Ch 5 

Ch 6 

Ch 7 

CN=7> 

7.6 

8.4 

8,4 

8.0 

Light Target 

(N=6) 

21,6 

23.1 

24.0 

23.5 


From five LANDSAT tapes over Mono ],ake (Table 1), 
similar brigJitness ratios were prepared for the 
light target to compare with the concrete 
(Miiffett Field), witii allowance made for their 
different reflectance values. These are tabulated 
in Table II IC and plotted cm the graphs of 
Figures 3, 4, 5 and 6 as open triangles connected 
by a dash-dot line. These lines are rub-parallel 
to the 1:1 cosine-law line, but art*, lower in all 
channels, probably due to slight differences in 
the absolute values of the reflectance of the pum- 
ice and the concrete, the ratio of which was used 
to relate the two data sets. (IZ error in either 
curve would change the Mono Lake data by 006 units 
at 180, Figs. 3-6). 

4. Application 

Our main rese.irrli task was relating the 
LANDSAT (lata to ground measured rcfJertanccs of 
vegetated (grass) targets — four main rock cypes-- 
on tJie 1850 hectare (5000 acres) of the Stanford 
Grassland Site (145 meters (475 feet) above sea 
level). We used the two calibration targets 
(carbon black and concrete) to reduce all our 
LANDSAT tape coverages over the site, to reflect- 
ance, .ind then compared them (whenever Lapes and 
field d.ita coincided!) with ground bidirectional 
reflectances. A reasonable agreement was found 
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and three of the areas have been selected from 
the 4A studied, and listed in Table IV. Channel 
7 invariably shows a calculated reflectance much 
higher (5-15%) than the values for the ground 
measurement, a feature also noted by Levine^^' in 
other similar measurements elsewhere in the Bay 
Area. 


Apart from this un-exp lained increase, the 
data in this calibrated ''reflectance^' form are 
c{uite useful for evaluation in clustering and 
other decl 5 lan*<making algorithms. 


5* Conclusions 

a. A very- low reflectance target has been 
located near Port Chicago, California and used 
for dark-area subtraction for atmospheric correct- 
ions. Ground measurements confirm the carbon- 
black material has a low reflectance (<0.5%). 

LANDSAT brightnesses have been extracted for 14 
overpasses in 2 years. Some cosine-law effects 
appear, 

b. A light- target (concrete, with about 30% 
reflectance has been used to extract similar 
brightnesses, which obey the cosine-law closely. 

c. Use of both targets (which occur on the 
same tape and scene) enable one to calculate 
"reflectances" from LAKDSAT tapes for any other 

set of pixels. jj 


d. These reflectances compare favorably with 
ground-measured data, despite the many locational 
and sampling problems. A positive bias in reflect- 
ance for LANUSAT Channel 7 data gives high re- 
flectance values. This may be due to problems 
with tlie low values of all Channel 7 data, to 
look“Up table biases in decompressing tiie Channel 
7 data, or problems with our basic concrete 
standard. 
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Dates 

and Sun Angles LAND5AT Tapes 

Used 

tA**IWAT 


Gun 

Zont 


in 

lintc 

elevation 

urwLl* 


lX3-im75 

07/2&/7Z 


3T.3* 


lOS*.-lU05i 

n'j/lf,/72 

47 

43 

h*na Lltfc 

mt-tmi 

10/72/7? 

36 

54 

tU/nO lll'C 

iiLii^ia75 

Ol/04/n 

24.2 

65, B 

StOlklODi 


01/72/71 

26.3 

63.7 

Stanford 


04/04/73 

49,4 

60.6 

StnnioTd 

l273-inifl3 

06/72/73 

5G.2 

14. a 

Stnnipr.J 

13Q7-1B0G6 

^ii/7r./71 

61 


ISnilO in’.rV' 

13a9'lB181 

ii^/7n/n 

61. 0 

29.0 

Stnnf^td 

U4S-1E174 

07/03/7J 

61. b 

2D, 4 

St.iniord 


<17/1*1/71 

59 

31 
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07/71/71 

59.0 
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Stanford 


O0/;4/73 

53 

37 
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J399-mi70 
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52.<l 

3B.0 
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Ik/24/71 

7/.0 

61.0 

St.trfPii’ 


U/V^/73 

23. 0 

67,0 


LLG^IQILI 

0'»/7J/74 

r.U.O 

30.0 

Staiifont 

lGU7-iaiD4 

Ob/ 10/ 74 

61.0 

29.0 

Stanford 


LANDS AT Radiances from Standard Targets 

Mean no. Standard <ievlatl<in 


_JD 

7* 

C114 

dl5 

016 

CM7 

P.4 ... 

as 

v6 

o7 

Moffett 

Field— Ught 

biirf;ec 

(cDiicretc) 

N - 14 






1365 

20,4 

75.9 

79.4 

67.9 

26.4 

3.97 

5.7fl 

5.16 

2.15 

1309 

29.0 

75.2 

77.4 

65.5 

26.6 

1.42 

2.61 

2.79 

1.0^ 

lfrH7 

20.0 

74. G 

77.7 

67.9 

26.1 

2,51 

2.;: 

4 .46 

i,il 

l(i69 

30.0 

70.1 

73.6 

63.7 

24.5 

2.62 

3.55 

2.55 

1.93 

13^3 

31.0 

72.9 

75.9 

64.6 

26.0 

2,50 

3.32 

3.66 

1.3f 

1003 

3U3 

67.7 

74.7 

66.9 

V7.3 

1,19 

2.72 

3.73 

la? 

1273 

34.8 

71.2 

76.8 

64.6 

26.0 

2.50 

3.32 

3.86 

1.3t 

1399 

30.0 

71.7 

73.8 

63.7 

25.6 

2.H3 

4,56 

3.00 

i.tr- 

1255 

40.6 

62.4 

G3.A 

56.7 

23.7 

2.66 

3.67 

3.26 

1.5* 

1075 

48.4 

61,6 

61.6 

52.7 

20.4 

1.29 

2.46 

2.06 

1.31 

1AH9 

63.0 

42.9 

47.2 

35.8 

14.5 

1.71 

2.31 

2.63 

1,39 

11B3 

63,7 

38.3 

38. 7 

33.2 

13.3 

1.64 

1.04 

1.40 

0.;: 

1165 

65.8 

35.9 

34.1 

30.3 

12.4 

1,31 

1.57 

1.18 

0.73 

1525 

67.0 

36.2 

35.3 

30.6 

11.8 

1.30 

2.35 

1.59 

0.B8 

Bcnccia Ref tnery—'d^rk target (carbon 

black) 

H " 50 





1365 

28.4 

21.0 

12.3 

10.0 

1.9 

0.96 

0,76 

0.79 

o.so 

13G9 

29.0 

19.7 

11.0 

B.I 

1.5 

0.72 

0.63 

0.62 

n.6i 

1657 

29 .0 

23.5 

13.9 

13.5 

2.6 

0.99 

0.68 

0.97 

0.5f 

1669 

30. D 

23.4 

15.0 

14.0 

3.7 

1.12 

0,82 

1.32 

0.53 

3291 

30.4 

22.1 

12.9 

10.7 

2.1 

1.4H 

0.88 

1.16 

0-59 

1363 

31. 0 

21.5 

12.8 

9.2 

2.0 

1.13 

0,77 

0.54 

0.73 
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31.3 

18.8 

10.9 

10.5 

1,9 

0.91 

0,50 

0.71 

0,56 

1273 

34.0 

20.3 

n.9 

11.0 

2.6 
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0.64 

0.85 
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1399 

38.0 

21.1 

12,6 

10.5 
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0.B6 

f),‘16 

0.75 

1255 

40.6 

18.6 

10.5 

9.3 

2.0 

0.7B 
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0.57 

1075 

48.4 

18.6 

9.8 

7.9 

1.6 

0.96 

0.40 

0-85 

0.95 

1689 

63.0 

n.6 

6.6 

5.3 

0.5 

o.eo 

0.97 

0.81 

0.61 

1IB3 

63,7 


area clouded over 
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0.53 

0.89 

0.64 

1525 

67.0 

14.2 

8.2 

5.8 

1.2 

0.79 

0.44 

0.S2 

0.68 

Mono Lokc-^llchC tar|;ct (puiLlec aiiinl) 

11 - 6 

(nfi<*r lUillfUf 

1975) 
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29 

59,0 

60.0 

50.2 

19.0 
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55.2 
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47.0 

17.5 
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50.5 

51.3 

44,7 
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43 

48.3 

47.a 

41.8 

15.0 





1091 

54 

40.0 

39.7 

33.2 

12,2 
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27.5 

24.0 

17.7 

4.6 

1351 

31 

27.0 

22.9 

19.6 

6.1 

1397 

37 

23.6 

19.7 

Ifl.l 

S.O 

1055 

63 

24.1 . 

20.7 

15.6 

4*9 

1091 

54 

20.3 

16.3 

12.3 

3.1 
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SXAMSORT: STANFORD REMOTE SENSING LABORATORY PATTERti 

RECOGNITION AND CLASS IFICATIOM SYSTEM 


F.R. Honey, Research Associate 
A. Frclat , Research Associate 
R*J*P* Lyon, Principal Investigator 

School of Eacth Sciences 
Stanford U^versity 
Stanford, California 94305 

ABSTRACT 

The principal barrier to routine use of the ERTS roultispeccral scanner computer compatible 
tapes, rather than photointerpretacion examination of the images, has been the high 
computing costs involved due to the large quantity of information (4 X 10^ bytes) contained 
in a scene* SXANSORT, the interactive program package developed at Stanford Remote Sensing 
Laboratories alleviates this problem, providing an extremely rapid, flexible and low cost 
tool for data reduction, scene classification, species searches and edge detection* Tlie 
primary classification procedure, utilising a search, with variable gate widths, for 
similarities in the normalized, digitized spectra Is described in section 2, with associated 
procedures for data refinement and extraction of information- The more rigorous statistical 
classification procedures are described in suction 3 * The programs have been developed 
on an interactive computer (PDP-10) with the non-specialist operator in mind, requiring 
very little computing experience for their operation* 


1* INTRODUCTION 

When confronted with the overwhelming quantity of data available on magnetic tapes from the 
ERTS-1 multispectral scanner system, it may appear to an investigator that reduction, analysis and 
presentation of significant interpretations of the taped data using a digital computer would be an 
expensive and time consuming approach* In comparison, visual examination of the standard imagery 
product generated by NASA from the original data (or color-combinations of the data) Is less expensive, 
though probably more time consuming. However, this photo-geologic approach can not be correlated 
readily with field-or laboratory-measured data (here referred to as "ground data")* 

Crouping data for a scene into distinguishable classes, for comparison with known ( oi suspected) 
geologic, geobotanical, crop or urban features, may be accomplished with either a statistical or u 
non-statistical classification procedure. These procedures may be divided further into supervised 
procedures (requiring training groups consisting of either areas which are known to be uniform, or 
of digitized "ground data"), and unsupervised (seif-training) procedures in which the data is split 
into its distinguishable groups with no prior knowledge of the number, or species, of these groups- 
The statistical classification procedures are reviewed and evaluated in section 3 . 


^Thls research report is based upon work performed under NASA Contract NAS 5-21884, 
the receipt of which is gratefully acknowledged* 
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2, KQN-STATISTICAL PATTEHN RECOGNITION, CLUSTERING PROCEDURE 
2,1 UNSUPERVISED CLUSTERING PROCEDURE 

The tjon-statistlcal gating procedure described below developed as a result of manual examination 
of digitized spectra plotted for a area using the four NSS bands* It was realized that* for an area 
of reasonable size (e*g. 30 X 30 pixels* where a pixel* 57 meters X 79 meters, covers approximately 
1 acre) only a finite number of patterns appeared. Figure la illustrates the patterns appearing when 
traversing across a row of pixels for a scene in Mono Lake, California (Figure 2)* Although overall 
lavels vary it may be ofaservetl that similar pattern s appear across this traverse* This variation in 
level is due mainly to a topographic effect: slopes facing the sun appear brighter than slopes facing 
away from the suni and partially due to the texture of the surface within the pixel: smoother surfaces 
(generally composed of smaller, closely-packed particles) appear brighter chan coarser*rough surfaces 
of the same material. This change In brightness level may be reduced considerably by normalizing to 
one of the channels. The effect of normalizing to channel 4 the patterns plotted in Figure la is 
illustrated in Figure lb. 

These normalized patterns could be grouped Into classes now, just by their shape. This Is a 
tedious approach when performed manually* yet the concept provides a very simple, rapid and economical 
technique when performed by computer. The computer is required to perform the m^nimum of operations, 
all arithmetical , merely comparing values within st.te preset range, to discriminate different classes* 

2.1.1 CONVERSION OF ERTS DIGITAL VOLTAGES TO REFLECTANCE VALUES 

To compare the satellite results with ground data it is necessary to convert the ERTS digital 
values for each channel to some more absolute measurement which w-il be virtually independent of sun 
elevation and atmospheric effects. For this conversion, two (or more) "standard** targets are required 
in any ERTS scene. One of these should have as ^qw a reflectance as 'ossible (prcferrably zero 
percent), so chat It may be assumed that energy impinging on the dc • from the direction of the 
low reflectance target arises only from Che radiation back scattered i to the sensor field of view 
from Che atmosphere. The four ERTS channel values for this low reflectance target may then be 
subtracted from the corresponding values for all other pixels within the scene to give a measure of 
the radiation impinging on the detectors which arises specifically from the radiation reflected from 
each pixel. Obviously* not all scenes have zero-reflectance targets within them-in this case several 
targets having low reflectances must be chosen, and a linear extrapolation performed to give reasonable 
values for a zero-reflectance target. 

To convert these corrected ERTS voltages to reflectance , a standard high reflectance target (or 
targets) within a scene must be chosen. By ratioing the corrected channel voltages for an unknown 
pixel to the corresponding corrected channel voltage for the high reflectance target, and then 
multiplying by the respective band pass reflectance (known from ground data) of the standard target (s> 
yields the band pass reflectance of the unknown pixels within the scene. This procedure may be 
represented explicitly by 

yu . rl - - . Vz>l X (1) 

Vx.i - Va.i 

where i is the channel number. Pu,i is the bandpass reflectance of the pixel being examined. Ps,X 
is the band pass reflectance of the standard targeu(s). Vu,i is the voltage in the ith channel of 
the unknown pixel. Vz,i is the voltage of the zero reflectance target Vs,l is the voltage of the 
high reflectance standard target. 

The low and high reflectance standard targets must be chosen to cover an area at least three 
pixels square, preferrably larger, so that the center pixel (or pixels) of the standard target give 
voltages arising entirely from the standard targets, not affected by bordering species, particularly 
if the area is to be repetitively sampled by other ERTS overpasses. 

Again* as would be expected, only a small, finite number of patterns appear when these reflectances 
are plotted* Levels of the reflectance patterns vary, due to the topographic and texture effects, 
but these variations may be removed by normalizing to one of the four channels (c.f. the ratio 
technique described by Vincent (1972)). 

The technique described above has been developed into a rapid, inexpensive clustering program 
for an interactive computer system. With the man-machine interaction the investigator can rapidly 
choose hla scene* display shadeprlnts as maps (for location) and optimize the gate used in the 
clustering to suit his particular requirement and the size or complexity of the area being examined. 

The claasificacion procedure searches through the array for the first unclassified pixel and a 
descriptor (alphabetic) assigned for this pattern. The remainder of the unclassified pixels ate 
then compared with this "standard" pattern, if the pattern of an unclassified pixel agrees with 
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Che cur teat ” standard” pattern within the gate widths it is given the same descriptor as the current 
"standard** « The program recycles until all pixels have been classified, or until Che number of 
classes exceeds twenty six (arbitrary)* The set of "standard” patterns generated during the search 
are stored in an array to be u ed, 1£ required, to classify another scene in the same area. In this 
manner very large areas may have a cluster analysis performed on them*. 

2,2 SUPERVISED CLASSIFICATION PROCEDURE 

tfhllst the unsupervised clustering technique described above is useful for examination of unknown 
scenes, separating them into their spectrally distinguishable species, a similar but more powerful 
technique may be used Co search the band pass reflectance patterns for known types, using the ground 
spectral data* The gate generally employed in this method is two or three times the largest standard 
deviation of the normalized bandpass reflectances of the measured ground target. Obviously for this 
approach the standard targets in a scene must be chosen carefully, according to statistical sampling 
technique, and their reflectances measured* 

2-3 NOISE AND SMOOTHING 

When radiance data for a large ^ uniform scene (eg. water) is examined, noise with a six row 
periodicity may be observed, resulting from the basic detector design in the* MSS. This noise is in 
phase for channels A, 5 and 7, but out of phase by two scan lines for channel 6* The origin of this 
noise is not clear-it may be due to slight differences in detector responses or to a misalignment of 
the detector array. This noise must be removed .is completely as possible for a reliable cluster 
analysis to be performed, not to do so makes each sixth line a separate class. Since it is in the 
farm of one or two unit spikes every sixth line, it is ideally suited to treatment by a digital 
smoothing technique described by Savitsky and Golay (1964). As the convolution blurs the image 
slightly, it must be performed in two dimensions in order that the change in contrast will be uniform. 

The result of this smoothing is evident in figures 3a and 3b* Figure 3a is the result of a 
cluster analysis without smoothing on an island Lti Mono Lake. Figure 3b shows the analysis of the 
same scene (with the same tolerance) af ter smoothing of the data* The water surrounding the island 
appears non-uniform in the unsinuothed cluster analysis result, but becomes uniform after smoothing* 

2,4 BOUNDARY SEARCHES- EDGE DETECTION 

After smoothing it is possible to search for abrupt changes in contrast, such as occur at sharp 
boundaries, deep valleys or borders between materials with large differences in reflectivity. In this 
manner, a search for linear, curvilinear or elliptical features may be pursued and, hopefully, some 
correlation between the presence of these features, their intersections and changes in classification 
using tlie clustering algorithm may be observed. Figure 4 illustrates the result of "edge detection" 
for the same scene clustered in Figure 3b. 

^•5 COST IK COMPUTATION 

The clustering technique as outlined above has proven to be very rapid. No direct comparison 
with a statistical clustering program is available yet, although the B^^D07^l atepvjise discriminant 
program hai! been employed and found to require approximately ten times the computing time, even when 
using training groups initially generated by our imsupervised clustering program. Obviously the time 
required for classification is a function of the number of pixels in the scene, and also the width 
of the gate. Figure 5 illustrated the times required as a function of Che area of the scene ln*t!ie 
vicinity of the island in Mono Lake for different gate widths* 

. At present the program is being extended to "defocus" the scene, so that large areas may be 
examined by averaging four, nine or sixteen pixels (in a square), clustering them to look for broad 
patterns, then examining sub-scenes of the large scene in 1 X 1 pixel detail. Statistical procedures 
are being inserted to provide means , standard deviations and histograms of areas classified by the 
clustering algorithm, or of areas selected by the operator. 

The program has been developed with the non-specialist in mind. It Is completely Interactive 
and seif explanatory so that a person with no computing experience is able to examine ERTS tapes* 

It is designed specifically for use with a limited budget, with fast turn atound time. 

3. SRSL NUMERICAL CLASSIFICATION TECHNIQUES 

This section outlines the theoretical basis of numerical classification technqiues used in 
conjunction with the procedure described above* Together they constitute the software system for 
analysis of ERTS multispectrai data in operation at the Stanford Remote Sensing Laboratory. Four 
numerical classification procedures are discussed, two of which are supervised and two of which 
are unsupervised. 


3* 1 SUPERVISED CLASS IF ICAIIONS 

A classificacion is supervised when data points of unknown origin are assigned into a priori 
defined classes. 

3.1.1 NEAREST NEIGHBOR 

Most of the classification techniques depend upon the assumption that samples have been drawn 
from a normal population. The Nearest Neighbor method makes inferences without any assumptions as to 
the form of distribution in the population. Such procedures are said to be non-paratnetric or 
distribution-free. The technique consists of classifying unknown data points into known categories 
through comparison with knoi;n data. Each unknown sampLu is allocated to that group to which it is 
nearest in terms of the generalized-distance si.itistic* Thus> the degree of similarity between 
two samples is provided by Che distance that separates them; the shorter the distance the greater 
the degree of similarity, and vice-versa. 

The nature of non-parametric statistical infc'cences usually requires testing with large amount s 
of data to achieve a respectable degree of accura.v (Swtizer el-al * , 1968) . 

3.1.2 MULTIVARIATE DISCRIMINANT ANALYSIS 

Multiple discriuinant analysis is a statistical method of assigning samples in a probabilistic- 
manner to previously defined populations on the basis of a number of variables considered simultaneously. 

The use of the discriminant function may be considered in terms of a population consisting of 
X variables, which forms a cluster of points in X-dimensional space. A second population^ described 
by the same X variables, consists of a second clusster of points. The linear combination of variables > 
that defined a multi-dimensional plane efficient ly separating the two clusters of points is the 
discriminant function . The degree of distinctness of the two clusters can be analyzed by measuring 
the ‘^distance*’ between theiv multivariate means. Once this distinctness has been established and 
the separating plane computed, additional unknown samples can be assigned *n one or the other of tiie 
groups depending on which side oE the discriminant plan*' they fall.. 

The basic assumptions about the data are; (t) the observations in each group are randomly 
chosen; (ii) the probability of an unknown observation helonging to either group is equal; (iii) 
the frequency discributions of the groups are each multivariate Gaussian distributions with a common 
covariance luatrlx. This means that the distributions have identical bell-shapes and differ only in 
^ha ^ they are centered at different points . 

The BMD07M is a stepwise discriminant analysis program, and is part of a series of bio-medical 
statistical analysis programs compiled by the UCLA Health Services Computing Facility. The stepwise 
discriminant analysis indicates that tlie computation of discriminant functions is not simply based 
on the original variables considered as a whole, hut rather that the variable- arc entered separately 
and consecutively by order of discriminatory power. The advantage of thi.s procedure is to recognize 
the relative importance of each variable in classifying the samples into t[w. different groups. Rank- 
ing the variables by predictive power permits a concent ration of efforts on those factors which are 
important for classifying groups, and this ran represent a highly effective means of reducing costs 
of data collection and processing. 

The computational procedure of the stepwise discriminant technique is described in the user^s 
guide of the BMD07M program (Dixon, W, J . ,ed . ,1972) . 

3.2 UNSUPERVISED CIASSIFICATIONS 

Classification is unsupervised when similar data points are placed into an unknown number of 
distinct classes in which the data points of each class have a closet similarity to each other than 
to the data points i ail other classes. 

3,2.1 CLUSTER ANALYSIS BASED ON DISTANCE-SIMIbVRJTY MATRIX 

A distance-similarity matrix is obtained to determine the relationship uf the data. ITie use 
of distance is based on the concept that a quantitative measure of the degree of similarity between 
two variables or t\-m samples is provided by the distance that separates them in a rectangular coordinate 
system. As the distance becomes shorter the degree of similarity increases and vice versa. The 
sample points are grouped or clustered in a hierarchical dendritic network Cdondogram) in which thfir 
interrelationships, as contained in the distance-similarity matrix, are shown with greatest simplicity. 

In a two-dimensional case, two samples ace plotted according to the values of the two variables, 

X and Y. The distance between these two points is, by simple geometry, ths square root of the sum 
of the sqiiared differences between X and Y values of the two points; as in a right triangle the 
square of the hypotenuse is equal to the sum of the squares of the two sides of the triangle. 

This calculation of the distance assumes that the input variables (or the axes from whicl* they 
are measured) are uncorrelated, that is, orthogonal or at right angles to each other- However, most 
raw variables are correlated to different degrees so that the coordinate axes would not be at right 


Angles And the simple Euclidean distance formula would be inaccurate « To overcame this difficulty, 
the raw variables are transformed to a new set of uncorrelated orthogonal variables by a series 
of linear transformations (for details see Sebestyen, 1962). In calculating the distance coefficients 
for the slcdlariCy matrix it is convenient to limit its value to the range 0.0 to 1.0. To satisfy 
this requirement I the original data Is transformed, so chat all the measurements are positive and 
range from zero to one. 

Finally e a cluster analysis is performed to measure the degree of similarity between samples 
on the basis of the distance-similarity matrix* Distances close to 0.0 represent icaximuni similarity, 
distances close to 1*0 represent minimum similarity. A cluster diagram is printed out with the value 
of the distance coefficients* Groups of similar samples can be selected at any desired level of 
aimllarlty , and each group can be plotted on a geometric matrix or map. 

The present procedure accomplishes clustering by computing a matrix to measure all pairwise 
similarities between data points on the basis of the measurements corresponding to the channels of 
the scanner* The procedures cannot be used when large data sets are to be analyzed because the 
size of the distance-similarity matrix becomes too large for the core storage requirement of the 
computing equipment* 

^».1.2 ISOMIX: AN ITERATIVE CLUSTERING PROGRAM 

Similar cluster programs have been developed by Stanford Research Institute (Ball and Hall, 
“ISODATA**, 1965), Purdue University (Wacker and Landgrebe, 1971) and Lockheed Electronic Company 
(Kan, Holley and Parker, **ISQCLS*', 1973). ISOMIX (Stanford) essentially follows the iterative 
clustering procedure of ISOCLS; however, new statistical techniques have been added to help the 
analyst in the interptetation and evaluation of the final data points grouped into clusters. The 
following Is an outline of the main steps; The program first computes the initial cluster centers 
and assigns them to regions of high sample-point density. Then the samples are sorted, one by one, on 
the basis of their distance from a set of initial cluster centers which create a cluster of data 
points or vectors X and ¥ is defined as 

d(X.Y) - I iX.-Yj (2) 

i«l® ^ 


After the samples have been sorted the mean and standard deviation for each subset in each 
dimension (variable) is computedi 


Those clusters which contain only a few sample points are discarded* Splitting of the clnscers 
takes place if the standard deviation in any dimension is greater than a suitable threshold specified 
by the analyst. When the cluster is divided two new centers are formed. These centers are (u ,U , 

.....a-kt. ... U ) and (w ,y , .,p. -O. , .li ) where(M ,U ) and (a ,c ,...,0 ) are the* 

mean and standard uevlation^foi the dimensions in the ori^ln^l cluster, and lA tSe ICth ^dimension 
the original cluster contains the largest stanuard deviation* 


The degree of distinctness of the clusters is measured by the similarity of ’’cluster centers" 
attached to regions of high density of data points. The distance or measure of similarity between 
two clusters and C^, where is characterized by (n^ ^ and standard deviation 

and by and (Kan, Holley and Parker, 1973), is defined as 




n 

I 

i«l 




2 n 


1/2 


(3) 


If the distance between two clusters is less than a specified threshold, the two cluster 

centers are merged into one at a weighted mean of the two original clusters. 

The progress is cycled repeatedly until the standard deviation In every channel of the generated 
clusters is less than the specified threshold, or Che maximum number of clusters desired by the 
analyst is reached. 

ISOCLS* a chaining algorithm is used to link those subclusters which are close to at least one 
other subcluster. - This linking process determines Che subpopulations, the union of which constitutes 
the parent populacloru 

In the last step the overall areal proportions various clusters are obtained* For example, 
if p* is the areal proportion of a specified cluster j, is the number of sample points counted 
of tne specified cluster j, and N the total number of sample points, then the usual estimator of 
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the areal proportion p. is p * n./N, Finally In the last step» the pattern complexity which gives 
the spatial scale of variation is^also obtained* A pattern that has a cluster A with its samples: 

In a contiguous body is less complex than another with the same proportion of cluster A distributed 
in many scattered smaller units* One index to express the pattern complexity (Switzer » 1973) .is 


X « total length of boundaries between different sample points 
(area of region) 1/2 

The value of X is invariant to the choice of the measurement unit. As the X value increases the 
pattern grows in complexity. 

The output gives the statistics for each cluster and Includes a map showing the final cluster 
assignments of all the points in the area analyzed. These maps are geographic matrices preserving 
the original position of the data points. 
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Figure 1. Digitized spectra resulting from 
plots of 


(a) Raw ERTS digital voltages 
plotted against channel number 


(b) ERTS digital voltages normalized 
to channel 4, plotted against channel 
number , 


for section of a scan line across 
Negit Island and portion of Mono 
Lake, California. 


Figure 2. Negit Island, Mono Lake, California 
Dark areas of island are basaltic 
lava flows and cones of varying 
texture, white "beaches" composed 
of calcareous tuffs. 
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Figure 4 4 Result of ^'edgs detection" on Negit 

Island scene. Only those borders with 
a large change in contrast are emphasized. 
Lower steps in the scanning technique 
would result In more detail. 


Figure 3. Cluster analysis result 

of the Mono Lake scene for 

(a) unsmoothed data 

(b) smoothed data 

The water surrounding tl«t- 
islatkd appears more uiiifonik 
for the smoothed data , 
However » some structural 
detail of the Island Is 
removed due to the 
"defocusing" effect of 
the convolution. 
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Figure 5. Computation Limes for 
unsupervised '.‘lassif-^ 
ication as a function 
of the number of 
pixels In the scone 
and of the gate widths 
In a first approach, 
gatewidths of 
approximately 0. L 
normalized units are 
used* 
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APPENDIX C 


CQHRELATIOM Q? ERT3 SPECTRA VJlTH 
RaCK/SOIL TYPES i:: CAIiIFQRtJlAN 
GRASSI,A:tD AREAS 


Saul Levine 

(Tech Report 7 ^“10) 

3tani"ord Remote Sensing Laboratory 
Department of Applied Earth Sciences 
Stanford University 
Stanford, California 9^305 

SUMMARY 

A seasonal study of ER,T3-CCT data, accomplished 
by means of four band spectra plots of normalized 
reflectance, indicates that in the San Francisco 
Bay and adjacent Coast Range grassland areas, 
soils mapping or classification by computer tech- 
niques is possible at the end of the dry or grass 
dieback season- Excellent correlation is shown 
between ground reflectance measurements and CCT data 
at three test sites and two different soil types: 
serpentine and sedimentary. The uniqueness of 
their spectra is then demonstrated by the success- 
ful application of STAUSORT, a computerised 
classification technique developed by the Stanford 
Remote Sensing Laboratory. 


I. IJtTRODUCTIOU 

Tne primary purpose of this investigation v/as to determine if the serpentine 
exposures and soils on the San Francisco Peninsula could be detected uniquely 
ty means of ERTS imagery and/or the related CCT data. In doing so it was also 
r*;ped to evolve a methodology vjhieh would be useful in conducting similar studies 
ir. the future. As envisioned, the imagery (individual bands and color composites) 
■.-.■ere to be studied first to determine if these serpentine areas could be detected 
visually and then a study made to determine if any iiniqueness existed in their 
four band spectra. This property, if existent, could then be utilized as a basis 
:''or the development of a computerised classification program to automize the 
i^vectlon and mapping procedure. 

In the course of this investigation, it became evident that the seasonal 
response of the vegetative rover could be most important in obliterating or 
rrhanclng the information - cttlng to the serpentine soil. Therefore, a careful 
3/.-::eratic, spectral study the CCT dat-^i for the yearly cycle v/as undertaken 
of four band radiance and reflectance plots of the test areas. Off 
correlations of serpentine soil spectra vs. serpentine soil/grass spectra 
■-ere also made possible by means of a fortuitous grass fire in the study area 
■■‘"Ich had exposed a large area of bare soil. These correlations ultimately led 
:: the conclusion that the soil/grass spectra were in fact essentially soil spectra 
the end of the dry or dieback season. 

After an extensive ground measurement program had substantiated the unique 
:JV3racter of the serpentine soil spectra the study was expanded to Include a 
’TTlnentary area on the east side of the Coast Range upon v/hich a yearly controlled 
- -r- occurred. The sane seasonal trends v/ere evident and a strong correlation 
ERTS reflectance spectra and ground measured spectra was again found. 
2-ddition, the spectra of the sedimentary soil was found to be distinguishable 
fro". the background as well as the serpentine soils studied on the San Francisco 




Peninsula. (See Figure 1 for general location oT the test sites-) The clustering 
program STA’^SORT developed by the Stanford Remote Sensing Laboratory was then 
applied to the study areas with significant success. 

rC. AREAS STUDIED 

Two major exposures of serpentine rocks and soils mapped by the USCS, on the 
San Francisco Peninsula, were selected for study- The first. Area I, consists 
of exposures of highly weathered blue-gray serpentine, only a small percentage 
of which is outcrop, the remainder decomposed fragments, grading to a serpentine 
soil. These exposures are to the east of, and adjacent to, the Crystal Springs 
Reservoir segment of the San Andreas Fault Zone* They are surrounded by and 
occasionally penetrate the various rocks of the Franciscan assemblage- To a 
large degree the northern section of Area I is obliterated by housing developments 
and roads. Therefore, the study was focused on the southern section which is 
largely within the Crystal Springs watershed and Is public land- The vegetation 
of this site, composed of annual broad leaved herbs and annual grasses, is readily 
distinguishable from the surrounding grassland on nonserpentine soil- It is 
marked by a different species composition, smaller size (height less than one 
foot), sparser cover and earlier onset of senescence and drying. The dominant 
broad leaved herbs are Layia olatyglossa (tidy tips), Orthocarous sp . (ov/l’s 
clover) and Plant ago erecta (California plantain) and the dominant grasses are 
Bromus mollis (s^ft chsssT and Lolium mulfciflorum (ryegrass). 

Area II is approximately ^ miles south of the Crystal Springs Reservoir, 
again on the east side and adjacent to the San Andreas Fault Zone. Because of 
the housing developments and roads an open field area of roughly 60 acres at 
the south end was selected for study. The serpentine is heavily weathered and 
blue gray in color vrith only a small percentage of outcrop] the remainder decom- 
posed fragments and serpentine soil. Interstate 280 transverses the south end 
of the area exposing large amounts of fresh serpentine in the roadcuts. The 
serpentine vegetation of the Farm Hill Road site Is clearly differentiated from 
the’ surrounding nonserpentine vegetation by the same features that distinguish 
the Crystal Springs Road serpentine vegetation, i.e*, a different species 
composition, smaller size, sparser cover and earlier onset of senescence. It 
Is made up of annual broad leaved herbs and grasses and shares several species 
in common with the Crystal Springs site. The dominant plants are a grass, 

Festuca so . (fescue), and the broad leaved herbs, Layia platvslossa (tldv tips) 
and Hemizonla sp . (tarweed). 

The third area studied consists of a sedimentary area located In San Joaquin 
County, California, at the base of the eastern foothills of the Coast Range. 

This area lies within the Lawrence/Livermore Radiation Laboratories Field Teat 
Site 300 and is largely a yearly controlled burn area. The study site is typical 
of the rolling grassy eastern foothills of the Coast Range. It is roughly 600 
acres in extent, crossing elevations varying from 250 to 400 meters* The sedi- 
ments are semi-consolldated sandstones with the outcrops again a minor percentage 
compared to the soils derlven from the sandstone. The vegetation Is that typi- 
cally described as a California valley grassland community, dominated by annual 
species of the grasses Bromus (bromegrass) , Festuca (fescue), A vena (cat) and 
others. 

III. VISUAL STUDY 

A visual study of available ERTS imagery, both the individual bands and color 
composites Goverino the San Francisco Peninsula was accomplished- Also included 
was U-2 Imagery taken during the EETS Simulation Program* It was noted, in the 
EP.TS frame date 6 October 1972, that a distinct dark gray pattern existed which 
seemed to coincide generally with the Area I serpentine east of Crystal Springs 
Reservoir. Study of the imagery before and after this date Indicated that the 
pattern persisted with diminishing intensity back to 26 July 1972 after which it 
could not be seen. The pattern was not ffvident again until 26 August 1973 » at 
which time it was faintly discernible* The appearance and 'disappearance of the 
observed pattern correlated with the dioback and growth cycle of the grass in 
this area* 

Review of the ERTS color composites substantiated the above, with the pattern 
readily discernible at the dates noted, as a dark purplish tone- In addition, 
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nini-iar tones were also e*/'ident within Area ±1, south of Farm Hill Roadj coinciding!; 
with the open field mentioned previously, 

IV, RADIANCE 3FSCTRA 


?o study the possible uniqueness of the tones associated with the serpentine 
areas, the radiance values of ERTS-CCT pixels traversing these and adjacent areas 
were obtained and their spectra plotted. These pixel spectra traverses, across 
"he Crystal Springs Reservoir and Farm Hill Road areas, are shown in Figure 2* 

The mean radiance values, standard deviations and coefficients of variation 
of the various terrain types, across these traverses were then computed^ Radiance 
throughout this report is presented as digital or v/ord count levels. Should 
absolute value of radiance be desired conversion factors must be applied. At 
this point no atmospheric corrections viere made. The location of specific fea^ 
tures^was accomplished by means of a shewed pixel overlay of the proper scale and 
ar. ortho-photomap (1:24000), as well as aerial photographs of the areas. 

It can be seen from examination of the spectra plots that the serpentine and 
grass areas as well as Interstate 28o, water and the forested areas appear to 
have distinctive spectra. It is interesting to note that while the pixel spectra 
across the forested area in traverse AA are generally the same shape, peak values 
are evident at four points. The aerial photographs indicate that these coincide 
with the hilly terrain across which the traverse v;as made. Apparently, this effect 
is caused by the variation in sun angle due to hill slope. The repetitiveness 
of the individual water spectra is also very striking. 

In traverse BB, the constancy of the pixel spectra of the serpentine soil/ 
grass area on the east side of Interstate 2S0 as contrasted to the west can be 
correlated bo the tones evident in the ERT3 imagery. The slight variability 
evident on the west side is attributed to variability in the soil, grass cover or 
br-";h. It seems apparent that the spectra obtained is a function of she inter- 
action of the soil and degree and type of grass cover which In turn is a function 
of the season of the year. 

V. SEASOflAL REFLECTANCE SPECTRAL STUDY 


Based on the results of the visual and radiance spectra study, it was evident 
that the serpentine soil/grass signature vras unique at the 6 October 1972 date. 

It also seemed possible that because the grass dieback in this part of California, 
was complete by this date, that the recorded spectra was essentially that of the 
soil. To substantiate the above, a systematic study of the soil plus grass 
interaction at the Crystal Springs area through the yearly cycle was instituted. 
Due to a fortuitous 15 acre grass fire, within Area I which occurred 1 July 1973 
and easily seen in the EHTS Imagery, it was also possible to Include spectra of 
the devegetated burn area in the study for comparative purposes. This study was 
accomplished by means of normalized four band spectra plots of the mean ground 
radiance values of the selected test areas. In so far as was possible, identical 
ten pixel areas within the ERTS frames covering an 11 month time interval were 
utilized. 


In order to compare results from the ERTS multispectral scanner data over 
this time period, corrections were made for the perturbing effects of radiation 
scattered by the atmosphere and the variation in irradiance on the scene with 
solar zenith angle. These effects were removed by studying selected targets 
of low (zero) reflectance and high known reflectance {Honey and Lyon, 1974). 

In the scene studied, a waste products treatment pond at an oil refinery near 
Suisun Bay, with bandpass reflectance of <0.5^ in all four bands was utilized 
as the zero reflectance standard. A concrete parking apron for aircraft at Moffet 
Field NAS California with reflectances of 27-8, 31-0, 30*0 and 32.3 percent’ 
bandpass in the four ERTS channels was used for the high reflectance standard. 

The factors derived were applied as follows: 


Target 


Reflectance =* 


Target Radianoe-Dump Reflectance (Meas) 
Concrete Radiance-* Dump Reflectance TPIeas } 


X Concrete Reflectance (Meas) 
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Radiance plors as well as normalized reflectance are presented in Figure 3* 
Study of this data revealed the following: 

a. Interpret ability of the four band spectra la greatly Improved by the 
application of the atmospheric corrections and normalization of the data to band 

b. The normalised reriectance of the soil/grass is at a maximum (particu- 
larly channels 6 and 7) at the height of the rainy season, 22 January 1973; 
roughly twice as high as that during 6 October 1972, near the end of the dry 
season* 

c* The normalized reflectance of the soil/grass gradually diminishes with 
the end of the rainy season and the entry into the summer dry-out period. 

d. In the burn area, on 3 July 1973s the normalized reflectance spectra 
drops to a minimum value, possibly as a result of both the devegetation and the 
residue of the carbonised- grass remains. 

e. By 21 July 1973 j the normalized reflectance spectra values in the area 
have increased, probably as a result of the dispersion of the carbonized ash by 
the wind. They now approximate the values at 6 October* 

f. A slight increase in the normalized reflectance is noted at 26 August 
1973 » probably due to some revitalization of the grass in the burn area. 

g. The 6 October 1972 reflectances are very close to those of the burn area 
at 3 July and 26 August 1973- 

Based on the above, a strong likelihood exists that the reflectance spectra 
of 6 October is that of the serpentine soil with little or no reflectance Intro- 
duced by the dead grass. 

A 10 pixel block within the Farm Hill Road test site, (6 October 1972, 
serpentine soil/grass) was also selected and the normalised reflectance spectra 
plotted in Figure A strong correlation was found with the reflectance spectra 
at Area I at the same time of year. Based on these results. It was decided to 
broaden the scope of the study somev/hat to see if the same trends were obtained 
at a third site in which the terrain soil was of a different type and at which a 
similar burn situation existed. This test site (Area III) was located at the 
Lai'/rence/Livermore Radiation Laboratories Field Test Site 300 at which controlled 
burns vjere used to reduce the likelihood of xmcontrolled fires as a result of 
explosives testing. It is on the east side of the Coast Range, approximately 15 
miles east of Livermore. The terrain type had been mapped as marine sediments. 

In addition to the study of another soil type, a comprehensive verification program 
of ground reflectance measurements of bare soil was to be accomplished at ail three 
test sites. Because of the size of the Midway site, approximately 960 acres, a 
grid pattern, at intervals of 10 pixels, was utilized across the bum area to 
obtain the reflectance spectra data. Study of this data again indicated that: 

a* The normalized reflectance is at a maximum in the winter and gradually 
diminishes with the end of the rainy season and entry Into the summer dry-out 
period. 


b. The 6 October 1972 reflectance spectra and that of the burn area are 
comparable. 

A comparison of the reflectance spectra for Area III (marine sediments) and 
Areas I and II (serpentine soils) (see Figure 4) indicates that they are consid- 
erably different, particularly in bands 6 and 7> thus leading to the conclusion 
that computerized classification could be applied. 

Data relative to the ground reflectance measurements taken at the teat sites 
are plotted in Figure 4. An Exotech ERTS Radiometer and scaled down satellite 
geometry were utilized to obtain this data. Figure 4 compares the ERTS-CCT 
derlven spectra with that obtained above. Study of this data reveals the 
following: 
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a* The ERTS-CCT normaliised reriectance spectra for serpentine soils, at 
6 October 1972 are almost Identical at Crystal Springs and Farm Hill Road. 

b. The ground reflectance spectra obtained at Midway correlate very well 
with that derived from the 6 October 1972 ERTS-CCT data. 

c. The ground reflectance spectra for serpentine soil at Farm Hill Road 
compare very favorably with the CC? data. The ground spectra for the roadcut- 

and outcrop-serpentines, while comparable to each other are substantially* different 
than that of the soil. It Is believed that because of the small areal extent of 
the outcrops as compared to the soil and the limiting resolution of the ERTS 
system, the outcrops have little integrated effect and essentially only the 
serpentine soil is detectable on the CCT data* 

d. The correlation of the ground spectra obtained at Crystal 'Springs with 
the ERTS-CCT data ir. not quite as good. It is believed that this is due to the 
infiltration of materials from adjacent soils derived from nearby Franciscan 
sandstone exposures. 

e. A significant difference is seen in both the ERTS-CCT spectx.. and the 
gr id spectra for the serpentine soils at Crystal Springs and Parra Hi 'JL Road 
an^ the sediments at Midway. 

In general, from the seasonal spectral study made and the ground measurements 
obtained, it can be concluded that the four band ERTS spectra for serpentine soils 
and sedimentary soils are sufficiently different from each other at the end of the 
dry or dleback season, to be distinguished from each other and the background by 
computerized clustering. 

VI. CLASSIFICATION TECHNIQUE 

The classification procedure utilized to investigate the uniqueness of the 
soils spectra studied is an Interactive program package called STANSORT, developed 
at the Stanford Remote Sensing Laboratories (Honey et al., 197^*)* This system 
provides an extremely rapid, flexible and low cost tool for scene classification. 

It Is a non-sfcatistical, unsupervised classification technique in which the data 
are split into distinguishable groups with no prior knowledge of the groups. 

The primary classification rirocedure utilizes a search, with variable gate widths, 
for similarities in the normalized or un-normalized digitized spectra. Training 
vias accomplished on the serpentine soil spectra of Farm Hill Road and the clus- 
tering resulted in coverage approximating the serpentine soils mapped. Continued 
search for the serpentine spectra at Midway (the sedimentary test site) using 
Identical classifiers revealed virtually complete absence of serpentine spectra. 

VII. CONCLUSIONS 

As a result of the foregoing study in the San Francisco Bay and adjacent 
Coast Range grassland areas the following may be concluded: 

1. ERTS soil/grass four band spectra are in fact essentially soil spectra 
at the end of the dry or grass dieback season. 

2. The ERTS four band spectra obtained is a function of the interaction of 
the soil and degree and type of grass cover v/hlch in turn is a function of the 
season. 

3. A strong correlation exists between ground measured reflectance spectra 
and ERTS four band spectra for both serpentine and sedimentary derlven soils* 

. The ERTS four band spectra for serpentine and sedimentary derlven soils 
are sufficiently different from each other and their background to be classified 
by application of STANSORT the SRSL interactive, unsupervised classification 
program. 


239 


REFEREMCES 




P, R* Honey, and J. P. Lyon, 197^, A Comparison of Observed and Model- Predicted 
Atmospheric Perturbations on target Radiances Measured by ERTS, Stanford 
Remote Sensing Laboratory, Stanford University. 

P. R. Honey, A. frelat, and R. J. P. Lyon, 1974, STANSORT: Stanford Remote 

Sensing Laboratory Pattern Recognition and Classification System, Stanford 
Remote Sensing Laboratory, Stanford University. 


240 







GENERAL LOCATION OF TEST. SITES 


pj«5PUCmi^ 
TrtfNAL PA®* * 




FIGUnE 2 ; RADIAHCS SPECTRA TRAVERSES (BY PIXEL) 
. ERT3 FRAi'tE ID 1075-18173 




72 22Jan73 19May73 3Jul73 3Jul73 2 UuI 73 ; 2lJul73 2EAug73 2EA\(073 

(Burn) (Burn) (Burn) 





CRYSTAL SPteS 
(Area I) 


FARM HILL RD 
(Area II) 


MIDWAY 
(Area 111) 


Group I 


Sp. outcrop 


Sp. roadcut 


COMPARISON OF NORMALIZED ERTS/CCT AND GROUND 
MEASURED REFLECTANCE 




APPENDIX D 




... 1 . 


r:UPHliLATIOIJ RETWKEN '^JROU:rD METAL ANALYSIS, VEGETATION REFLECTANCE, 

AND ERTS BRIGHTNESS OVER A MOLYBDENUM SKARN DEPOSIT, PINE NUT MOUNTAINS, 

WESTERN NEVADA 


R. J. P. Lyon 

Stanford Romot** Seriating Lab. Tech. Publ. 75-16 
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1 . ABSTRACT 

In a cooperative study vjlth U.S.G.S. and American Metals Climax (AMAX) per- 
sonnel, it has been possible to detect a «?.0 by 1 mile anomaly on ERTS CCT data 
nJrectly, in the pine- and .jun i f> M*-covered mountains of western Nevada. This 
arioiiialous area is. about 3-5 times larger than that of the known geobotanical 
anomaly which lies centrally wiUiin the area. The site has been studied on the 
grruncl and bi-directional r*f'f le«-tances (relative to BaSO^ ) obtained for 40 trees, 
'»slng both in-vivo techniques (similar* to cherry picker operations) and field 
determinations (d* cut brandies , The anomaly can be seen best by color trans- 
parencies made fi*'>m 5/^i , 6/H , 7/^1 ratloed ERTS data, the 3 ratios each being coded 
by one of 3 colors (blue, green, and red) (Figure 3B). Ratio-image 7/^1 is tiie 
best single view of these three »*atlos (Figure ?B). 


BACKGROUND 

Molybdftium m Ineral i satlon is. pt*cs'^‘nt as a skai'n deposit in Upper Triasslc 
limestones and mafic volcanics in .'lections 23y ?5 and ?6; T 12 N; R 21 E, 
approximately 2D mils GE of dar:-)n iMty, Nevada. DoLong (1971), of AMAX, In a 
r-.'iper pr*esented orally at an AIMK meeting at Vail, Color*ado , described the deposits 
as tactis.ed Limestones from whi'-h Uiere has been sporadic tungsten production 
pclncipally during World War* II from t,he Alpine (Mill), Divide and Cherokee mines. 
Work by AMAX resulted in defining a "sizeable molybdenum deposit" related to a 
quarts monziuiit* (Aipino Mill stock) in Section 25, although no grade or tonnage 
figures have been I’eleased. Geo(!hemical sampling of the soils was hampered by 
the lack of mobilit.y of Mo in calcareous environments. Vegetation sampling of 
tir^ deeper-rooted pines and junipers was used to define tiie distribution of the 
heavy metal by analyzing leaf-ash (needles, and second-year twigs). Junipers 
wf-r^' more usrfui than the pines because of ''heir highei* metal contents (DeLong, 
1971). Sui'ser|ucnt work by the U.3.G..S. has further documented this pattern 
(W'O :;on, 197/1). 

Mi> I yhd*.‘fiufn t.iuids lu be mohil*.-, liowever, in tLie volcanic soils and moves 
«b‘v;n i fie hydraulic grade out i rd •) Uie recent gravels rnd alluvials of Pine Nut 
t. * t.he riorUiear.t. In U:i.; area pines and junipers also shov/ high Mo 
coritt'nts . 

I’lius tn mi n-'ral i ::ed ii*ea h.a.* a well-documented geebotanleal anomaly witfi Mo 
as i-'ir*ing in I h. jurii{>**r riee<llps up to >500 ppm and in the pine needles to ^200 
: pm. dm*. .]*i uisr i . iy for our purposes no analyses of trie ubiquitous sagebi'usi. 

sfis'usis inadt'. 

R. Vi'it e )ii arui F. t’anney of t.he U.S.G..R. discussed this botanical anomaly wltfi 
the aur/f.L' * April 197^1, hut at Uie Ime w*irc not specific as to the locat ^ -n 
In tlie elic Nut Mountains, bo maintain company confidentiality. We suggesf.t’d 
using our .t EHT.S CCT tapes and the 3TANL0RT interactive program (Honey, 

et . al . , 197/i) to condU'!t a search for the section locations based upon ERT.3 
ref lect.ancM data alone. This proved successful in indicating an N-S elongated 
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elliptical area of about 2 square milea^ times larger than the known vegetation 
anomaiyj ahqi centered on the wiestern edge of Secliion 25 (I^yon and Honey t 1974) . 

Our results were sent to \fatsoh In Ha^ and field checked July 24-26* 1974 with hini^ 

This paper deals w^th subsequent att^ to document by field reflectance 
measurements that we had defined fz^om ERTS 


3- VEGETATIVE COVER 

The central area (Alpine Mlll^Dlvlde and Cherokee Mines) Is on the east facing 
slope at an average altitude of 650ti-^7Qflp teet about 1500 feet of relief. 

The hill Is bounded on the west by Buffalo Canyon. To the east across Pine Hut 
creek the main range of the Pine Nut rises to 9^ with Ht. Siegel at 9450 feet. 

The principal source of moisture are the winter snows which dust the lower slopes but 
leave Mt. Siegel snow^-covered for several months a year. The area lies in the 
^'precipitation shadow" of the Sierra Nevadas. Minden* Nevada* 11 nllies to the north- 
west and at 4600 feet averages 22 cm (8.5 Inches) per year* water equivalent. 
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The vegetation is mainly pinon pine (Plnus mbnophylla) and juniper (Junlperus 
utahensis) (Billings* 1950). The trees average 3-4 meters In height and are 
moderately bushy In shape, A morphological change of the pliie to a excesslvely- 
twiggy* sparsely-needled* brittle form was noted as the molybdenum content increased. 
The Junipers appeared more healthy with metal contents* 

Significantly for remote sensing* these trees (at 4Q per acre in unaltered 
volcanic s) covered only about of the ground surface* when viewed from above 

using air photos at a scale 1:40*000. The tree-frequency^count varied somewhat in 
the Immediate Alpine Mill area* showing about 30-35 per "acre, but they are more 
closely bunched* leaving more open space* which appeared white on the airphotos. 
in the unaltered areas the space between the trees is a light gray on the photo 
indicating a more-even cover of smaller shrubs like sagebrush. Elsewhere In the 
carbonates about 0.5 mile to the west of the Mill, the tree count was as high as 50 
per acre, still on the E-facing slope. No bare patches could be seen. 

By far the greatest cover is from sagebrush (Artemesia trldentata) * but because 
of its pale color it is difficult to estimate a percentage from the available black, 
and white photography. Visual estimablon would range froi». 10-40?. but as the shrub 
has an open* sparsely -leafed form, mostly soil and plant debris would be seen from 
above. This becomes important as it appears that a greater tree spacing* with an 
Q increased bleaching of the soll*determinea the "ERTS-average spectrum" of the anomaly. 


4* PREVIOUS WORK 

At the central area the earliest efforts in World War II were mainly connected 
with tungsten mlherallzatlon* American Metals Climax (AMX) spent considerable time 
mapping in 1967/68 and in studying the Mo geobotanioal anomaly, culminating in 
drilling the prospect (DeliOng, 197D • No logs were made available . 
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R. U, Watson and T. Hessin of the U.S.G,S. conducted measurements of the anomaly 
with the Praunhoffer Line Discriminator (FLD)* a device to measure secondary fluor- 
escence from chlorophyll in leaves as seen In a black solar line (Praunhoffer Line). 
They were in the field with us July 24-25, 1974 and used the Identical branches cut 
from the marked (analyzed) trees as we did for our EXGTECH measurements. Watson 
(1975 a-b.) has indicated a high degree of correlation with PLD-fluorescence 
anomalies and Mo content of leaves, and significantly, was able to repeat the pattern 
of anomalies in a second program of helicopter-FLD measurements in 1975 (Watson* 

1975, pers. comm. ) 


5. ERTS DIGITAL TAPE ANALYSES 

ERT3 GCT digital tapes* for scenes 1289-18063 (5/8/73) and 1397-18051 (8/24/73) 
were already available to us at Stanford and these we processed to enhanced-lmages 
using our STANSORT program (Honey* , 1974). We quickly found that all four 

bands gave somewhat similar images * principally showing the topography of the central 
ridges (Alpine-Divide and Mt. Siegel). By dividing any of the three bands by the 
fourth (usually Channel 4) this "sunlit- and shadowed-enhancement" of topography 
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Gould be removed from the digital imagery giving a "flat 'earth" presentation, the 
ratio- Imagery retalrilrif^ essentially only the spectral (colored) components. If one 
then combined any three of these black and white ratio images, using colored filters 
(red, green and blue) onto color film, a "color-ratio" image, enhancing all the 
spectral Informaitlori may be prepared (Levine, 1975)* These color-ratio images hold 
the key to Uie information content of the ERTS system, the image form being very 
significant to the geologist for photolnteroretatlon, and anomaly location in the 
field. 

Figure 1 shows enhanced ERTS images for Ch Gh 6 and Ch ?■ 

Figure 2 A-C show B/W ratio-images (Ch 1/Gh 4) and Figure 3B, a B/W copy of an 
orlginally-Golored ratio image (Ch 5/4 R; Ch 6/M - G; Gh 7/k - B). Figure 3A is 
a B/W copy of a false color image (Ch 4 + Ch 6 + Ch 7) , 

The Important anomaly covers about 2 square miles and Is centrally located in 
these figures, forming an elliptical to pear-shaped mass elongated in a NNE 
direction across the NW trending Alpine-Divide ridge. The eastern edge is formed 
along Pine Nut Creek, the southern edge by the low pass between Pine Nut and 
Buffalo Creeks, The northeastern limb is parallel with a strong Jointing pattern 
across the NW strike of the limestones, and may be also the location of two NNE 
fan Its 4 


The feature is subtle and varies In strength with differing combinations of 
ratios although remaining fixed in geographical position, centered Just north of 
the mutual Junction of gectlons .?5, 26, 35 and 3^- 

Tiieae images cover about 300 square miles at an original scale of 1:210,000. 

An important feature of their areal coverage is that other anomalous areas may be 
LdentU’led simultaneously , either with similar or differing "color-tones". 

Another ].arge anomaly, centered on Double Springs Fiat along Highway 395 is also 
enhanced on the same color-ratio image. This has a circular shape with a central, 
alluvlal-filled depression and several radial and concentraic drainage anomalies. 

We propose this as a caldera (andesitic) centered in section 23, T 11 N, R 21 E, 
and is particularly well shown on the Ch 6/Ch 4 image (Figure 2B) , and on Ch 7 
image (Figure 1C, arrow points outlining the shape). 

V/lthin a circular anomaly to the south a smaller one was located on Section 2/11 
by using cur search pattern In STANSORT/C LUSTER, trained on the Alpine-Divide area. 

A field traverse in this area located a new patch of weak gossan with many of the 
general appearances of the Alpine area. 


6 . FIELD MEASUREMENTS 

We have taken extensive sets of ERTS-band reflectance measurement, using a 
pair of EXOTECH-100 ground radiometers on 5 field days in the area. These were 
concentrated mainly on the 60 analyzed (and tagged) Juniper and pine trees along a 
drllllng-road down the N-S axis of our color anomaly. A map of the Mo analysis 
results was provided by R. Watson with data on Mo contents of needles, twigs and 
stems of these marked trees. 

Our Joint field measurements with Watson on July 24, 1975, ensured that the same 
trees were used for both sets of experiments. 

The reflectance measurements have been made with several sets of geometries 

1) Cut branches, observed vertically (front and back) at about 1 m with an EXQTECH 
15^ FOV unit, giving concurrent bi-directional reflectance relative to BaSOj^ 
powder, 

2 ) Living (in-place) trees, observed vertiaally down with a boom-mounted radiometer 
5 m above the gT'ound (concurrent bi-directional reflectance relative to BaSOq 
powder. 

3) Living (in-place) trees, observed horlgontally from 3^30 m distance, with a 
150 FOV unit. Bl-dlrectlonal reflectance. With intermittent (10-20 minutes) 
observation of a BaSO;^ standard. 
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^rraverses* with ineasurefnents every 30 ni, vertically viewing soils ^ rocks, sage-bruati and 
trees (horizontally), with a 15^ FOV unit, Interaiittently viewln^^ a standard. 

Traverse A; Alpine Mill ridge, to E to V/, over 0.6 inlles length. Ti’averse Bj Double Spriru-^^ 
anoitialy (Sec. 11), with 0.5 tniles in lenfdh- (Sept. 25, 197^)* 


7. DATA ANALYSIS 

7.1 Correlation of Reflectance with F^folyUicnuni Contents ’ 

The four EFTS bandpass bri^tness vabies for each pass target and BaSOi^ standard views 
were reduced usln^^ our EBTSRATS progrruu (Salem, 1975) < "Vbis yields several types of 
output; 

'H 

7.1.1 Target and Reference Standard separat« tly 

a. Bandpass radiant enilttance (w/cm^) 

b. Spectral radiarit end.ttance (vj/cm^/um) 

c. Bandpass ratios (7/6, 7/5, 7/^U 6/5, 6/4, 5/4) (Vincent, 1973) 

d. Specti‘^ ratios (7/6, 7/'>, 7/4, 6/5, 6/4, 5/4) (Vincent, 1973) 

e. Fseudo-C.I.E. cooi^inates (usiru^ Ch. 4, 5 and 7, instead of visible-region data) - 
this is a running check on the sky-color for error detection}. 

7.1.2 Using both Target and Reference Data 

a* Bidirectional reflectance (If botli have- ^aire PQV and ^«:eoiieU'y, or else direct ioiial 
i-enGctanre If hemispherical icv'/rllance, 21, is useil inst-eaa of a reference stai^ard. 

b. R eflectance ratios iCh^A; OH^/!-; CHj_/^ GHjA'; 1 = 4 ....7). Six non-redundEuit 
ratios Mtiy be pt^ared fr«3ir roui'-diannel data. 

Using the "cut-branch^* set of data following arjalyses been coiiipleted, to date 
(Table I). 

A cross-correlation (i-ir»P‘.v2D) analyst-- of Flo (ppm) and log versus 

a. Pseudo CIE cooi*^llnates RASMX arui ilASliiY (tarf^et). 
b- Bandpass brl[^tness BP4, UP5, BPn, RP7 (target). 

c. Bl-dii*eGtional reflectance R4, R5i Rb, H7* 

d. Reflectance ratios R76, R75, R74, Rii5, li64, R54. 

The results, tabulated in Table I sh'.w tlia? the cfjri»elal ions arc sh;?il ncniitly t> ordt lvr 
v;ith I'lo content for* Junipet^ and ^/g with log Mo ^-ontcnt foi‘ pines. 


Slpjnlficant at 

vjjth 

Parameters 

"r" 

1% level 

Juniper 

vs. 1-V> 

+0.1)2 


All branches 

b£‘T vs. i‘lo 

+0.35 

355 level 

Fine 

R4 vs .lor' Mo 

-O.llO 

555 level 

All bifinohe-’ 

R?7 vs. log M? 

+0.30 


n 

R7 vs.U'r III 

+0.25 


u 

R7 vi>. I-V. 

+0.2P 
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CoT’relation of Reflectance I tati r:- wli^h Anomalous Molybdenujn Content In 
Leaves 

Six Ratios of the four refi^nMan.^e values RT5> R7^l, R65> 

H5^) were calculated for pt« iU. ion of the t ree samples relative to the 
Mo anomaly boundai‘’les^ (!.'?. lir:lde or outside the anomaly). For pines, 
needle assay values above ^?0 ppir and for juniper values above 75 Ppm were 
considered anomalous. The .hree different viewing geometries . (days 1, 2, 
and 3) were kept separate tn avoid obscuring special aspect-dependant 
effects. Table II lists the data for the cut-branches. Table III for the 
vertical data and Table IV ^‘or the horizontally-viewed data. Ratio means 
(x 1000) and ocefficient of variability (COV) , where COV - o/x are shown. 

The following points may be madei 

a. All three aspects have very similar ratios, with R75, R7^* R65 showing 
the highest values. 

b. Variability within a group ratio (COV) is above 15 to 25^, higher in 
the R75 and r 65 ratios. Greatest variability occurs in the vertically 
-viewed trees where location of the radiometers relative to the needle 
masses is most difficult. Conversely, lowest COV values appear in the 
horizontally-viewed data whei'e location in the field of view of an 
instrument is most easy. 

c. The coefficient of variability of the mean s of each group may be used 

as a "criterion of separability" for each ratio tOoetz, et ^ 1975) 
l.e. in Table U, whlob t^onuaUis some soils and 16S sagebrush In 
addition to the pines and Junipers, a COV value of 0.39 for R75 makes 
H the most-useable ratio (0.07 for R5i» Indicates the least use- 

able one) . 

d. The three most-usable ratios (R75, H?4 and R65) show higher mean 
values Inside the anomaly for all the pine spectra, and also for the 
cut-branch and horizontally-viewed Juniper. Again variability exists 
for the vertically-viewed data, especially for the juniper. The 
horizon tally -viewed set (Table ^0 shows drops in mean ratio values of 
20 to 37? outside the anomaly, for both pine and Juniper. 

* . As R75, R7*^ and R65 are corrc^lated above 8c? with each other, any one 
may be used, in preferance ratio R75 because of its oonsistantly 
higher COV. Color ratios should be made with R75 or R7^» R6^ and 
to maximise their Inforniatloi! content, R5^ being In some ways a 
negative of R6*^. 


8. ;^I1MMARY 

Field reflectance measurements uf three modes were made, using EXOTEOH ERl'G- 
typ^ radiometers — cut branches, and vlevjing the trees both from vertically 
above, and horizontally. Each tree, eith^'r a plnon pine or juniper, was one 
previously marked by the U.S.H.S,, who provided the molybdenum analyses of stems, 
tv;igs and needles (leaves). In addition sagebrMish and bitterbrush shrubs were 
measured together vritli their background soils and rocks . 

The correlatii.'n b<=tween Mo and .’unip* r cut-branch reflectance was positive, 
and significant at. the 99? level (Chr-nnel 7 brightness) agreeing with the vlsua. 
nhservfit Ion that even at values in •• cess of 500 ppm Mo in leaf ash, the Junipers 
vjoro healthy. V/ith pinon pine however, the correlation with cut-branch (needle* 
t‘<H’ier’iaxico vjr’is norrai. ivc- but significant at the 97? level. The pines showed 
;■ I’l -unt nifjrf>hQ toKical change (needle loss, profusion of twit’:»''y stems, and 
f;rl tl Lr:ri‘.-‘sr* of branches) correiuteable with mineral uptake of Mo. Viithln each ret. 
the iuc Titian uf the anomaly can be found statistically by the R7«j, R65 ci* Rb^l 
l^-’vels, belni'; higher within the Mo-rich areas. 
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Using unsupervised clustering techniques on CCT taped data (STANSORT 
program) ER^ spectra could be extracted for the total anomaly area, which were 
used to locate similar areas to the south, near Double Springs Plat, Field 
checking located weak gossan mineralisation in the bleached andesites there. 

Continuing field studies are aimed at specifically identifying the cause of 
the ERTS anomaly — is it tree vigor, tree species, tree spacing, or sagebrush/ 
soil ratio which can be observed from spare over this skarn zone. 



LOCALJTy MAP of ICRTD-delected an'mily (heavy dashed lines) and tfiolybderuim 
anomaly (h- avy yt-ipplr^d area at»ovr 7‘> Ppm Mo in Juniper needle ash), 
riertiona given are In TliJN; H.'UK, Ml. ^Siegel quad, Douglas Co., L'allf* 

JIW- trending ltru*urs arc t!*eridi' o\ limn; tone beds; NE-t rends are slronr;ly 
developed Joints (fault;;?)* (^uat’lp; im^nzonite outcrop is solid black. 
Scale about 1:30, 
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i i. .1 . .■ 


S.-. iJs'satSi.i'^kSvS.'iL'C’. 
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TA‘if.E 1 

CORRELAfiOJJ BETWKT-JM hOa Mo (fipra). Mo 'ppm) WITH VARIOU:?: HKKr.KCTAMCE PAHAMETERr- 


RASNX 

RASNY 


Log Mo 
All 

Branches 


Junip^ r 


Branches 


Mo (ppm) 


Juniper 

M»33 


BP^t 


-.06 

-.28 

-.22 

.Oil 

-.16 

.28 

BP5 


-.Oil 

-.05 

.07 

.Oil 

-.20 

.31 

BP6 


.17 

-.03 

.29 

.13 

.07 

.13 

BP7 


(.30)» 

.19 


C.35)*« 

.07 

(.43] 



.13 


.23 

-.07 

-,2i| 

.19 

R5 


*11 

-.18 

.07 

-.02 

-.26 

.25 

R6 


.08 

-.15 

.27 

.00 

-.02 

.01 

R7 


C.25)* 

.11 

.32 

(*26)» 

.00 

.32 

R76 


.11 

.17 

.0(5 

.18 

.00 

.24 

R75 


.?n 

.21 

. i; 

.09 

.16 

-.05 

R7'« 



.28 

.It; 

.19 

.11 

-13 

R(>‘; 


.10 

.1)8 

.06 

-.Oil 

.li| 

-.22 

R6'J 


.13 

.19 

.or 

.03 

.15 

-.16 

RS'J 


-.Oif 

.06 

-.13 

.02 

-.09 

.12 

Signif- 

icant 








(»») at 

lit 

.3^5 

.'<7 

.i»5 

.33 

.47 

.45 

(*) at 

5i! 

.25 

.37 

• 35 

.25 

.37 

.35 





... 


(Snedecor^ 


^gD®^*S,l*WO? 






o 
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TABI.E 2 




32 ^ 
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TABLE 3 

VARIATIOMS IM RATIOS OP USO'JMC SPEOTRA-VBETICftL VIEWiaq 


DAY 2 DATA, July 25, 197** (IN-VIVO) 
ll7< spftctra of vegetanlon; 3* soils (cov « o/x, in parenthesis) 

VERTICAL VIENINU H RTd HT5 ^ R55 R51T 


Fines ; 

Inside anonaly >50 ooia Mo 22 1572 **803 5630 3056 35**1 i: 

’ " 1.13) (.26) (.28) (.24) (.20) ( 

Outside anOTTialy 30 I045 **310 5236 2ol5 3173 

(.13) (.26) (.22) ■ (.22) (.16) 


Juniper ; 

Inside anomaly >50 ppm Mo 31 -397 **^7^9 5530 30^8 34^17 

(.11) (.16) C.^1) (.18: (.15) 

Outside anomaly 15^^^ SUHS 5532 315? 3^W 

(.11) (.30) (23) (.30) (.17) 


Sagebrush 6 1377 1812 235? 133= 1796 ) 

(.13) (.11) (.09) (.14) (.11) ( 


Sells and Rooks 4 1174 1599 193^ 1360 1649 

(.04) (.1 ") ^.U'; ml;- (.12) 


tdtal data set 


127 


i:;82 47'-,2 4^73 3475 

(.13' (.42) (.40) (.35) 


2874 

(.31) 



255 


1 


<• I 

\ 


;jT 


VARIATI'-.i:-^ Zli 

R.H?:.;i OP 

iR:;nriD otPOTRA - 4 oh:i:::tai 'izv/::]*; 

Pine ;I'^t 

’•’.c j^nz 3iUrz , 

A’jvada, Pay 5, August 15, 157- 

c2% 

:3p‘^?’:ra oV 

^vegetat i on ; ZdZ of soils 


; ?OV s ?' 

X, In paT'^nthesia } 


Horizontal VI -v; 


rjt 

R75 

^74 

PITS' 

Rc-4 

R59 


Fines i 

Ins l ie anoziaiy 

15 

1Ac7(.26 : 

-211 (.26) 

551S(.03) 

33 •9u03) 

3735( .20 • 

1Q67(.09) 


Inside ancmaly 


1J<79(.05) 

5015(.15) 

5395,.,,, 

34li4 

• - ll *- 

3d60 
. 12 ) 

1076 

(.06) 


Outside 

13 

1399( .05) 

3793! .03) 

^209 

(.11) 

2"08 

(.06) 

3000 

(.05) 

1109 

(.08) 


Juniper: 

Inside anomaly 
(analyzed) 

25 

1993 

(.09) 

5031 
( .15) 

5627 

(.12) 

3368 

(.19) 

3762 

(.11) 

1123 

(.10) 


Outside 

20 

1921 

(.05) 

9096 

(.21) 

9396 

(.19) 

2615 

(. 18 ) 

3061 
( .15) 

1095 

(.05) 


Sagebrush 

- 

1192 

(.27) 

2090 
( .15) 

2290 

(20) 

1616 
( .12) 

1807 

(12) 

1122 

(. 08 ) 


Soils and Rocks 

25 

1077 

(.07) 

1153 
( .12) 

1505 

(.09) 

1069 
; .09) 

1396 

(.06) 

1307 

(.05) 


Total Data Set 

i^I 

1369 

(.1) 

38£’Jj 

(.39) 

9217 

(.35' 

2699 

2969 

(.30) 

1123 

( .07) 
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Figure? 3. LEKT(A): Falise Color Imago, prepared directly frcin ERT3 trjpe.s, 

adding CHll+CHf>+CI17 , closely resembling a normal alt* 
photo. White Is snow, lake in lower left Is Double 
springs Flat, white line Is Highway 395. 

RUJHT(B): Color-Ratio Image, prepared from ratlolng ERTS br Ightnes.? 

values. Snow area appears unusual due to saturation 
on the original data tapes. Anomaly is Indicated by 
ari’ow points. 


